{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "58ecd0f8-90fc-429c-bb36-1d5002ed3972",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import os\n",
    "import matplotlib as mpl\n",
    "from scipy.stats import lognorm\n",
    "from scipy.stats import beta \n",
    "import scipy as scp\n",
    "import matplotlib as mpl\n",
    "\n",
    "mpl.rc('font',family = 'Times New Roman')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "edfad454-d769-4777-a452-2946e2ed7c7c",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "d = np.zeros(30)\n",
    "\n",
    "for i in range(len(d)):\n",
    "    disp = np.loadtxt('./Results/dispI_GM'+str(i+1)+'.out')\n",
    "    d[i] = max(abs(disp[:,1]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "3c79e4d7-9bfa-47c3-9ffb-320086573219",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "d = np.sort(d)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "29206d0a-973b-4b29-bddd-0c993b20a1e9",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.0173556"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.median(d)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "0e3c9007-42c8-41d6-bfd3-94d0c4fedfca",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "pc = np.linspace(0,1,len(d))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "a872d1c9-2b2b-4007-bcf5-a8eef699b595",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([0.        , 0.03448276, 0.06896552, 0.10344828, 0.13793103,\n",
       "       0.17241379, 0.20689655, 0.24137931, 0.27586207, 0.31034483,\n",
       "       0.34482759, 0.37931034, 0.4137931 , 0.44827586, 0.48275862,\n",
       "       0.51724138, 0.55172414, 0.5862069 , 0.62068966, 0.65517241,\n",
       "       0.68965517, 0.72413793, 0.75862069, 0.79310345, 0.82758621,\n",
       "       0.86206897, 0.89655172, 0.93103448, 0.96551724, 1.        ])"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pc"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "e5a76820-8495-44cb-9834-52679a8b2500",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkQAAAG9CAYAAAD0lWkWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAACG6UlEQVR4nO3dd1gU1/s28HvpHQtYUCk2xF6QaGJBLAG7BjH22LuisWu+9qjR2Fts2GLB3rECURNFwS6KBUTALr3DnvcPf8zLSpeylPtzXXvpzpyZeWZnl3125sx5ZEIIASIiIqISTEXZARAREREpGxMiIiIiKvGYEBEREVGJx4SIiIiISjwmRERERFTiMSEiIiKiEo8JEREREZV4TIiIiIioxGNCRERERCUeEyIiIiIq8QpNQnTu3DnY2NggICAgwzYnT57ExIkTMWrUKFy6dKnggiMiIqJiTU3ZAQDA+/fvkZSUhFu3bmXY5smTJ1i0aBFu3rwJIQSsra1x6tQpVKpUqQAjJSIiouKoUJwhKleuHDp16pRpm9WrV8Pe3h4ymQwqKipo3rw5Nm3aVEAREhERUXFWKM4QAYCKSua52ZUrVzB9+nTpeY0aNXDkyJEM28fHxyM+Pl56LpfL8fnzZ5QtWxYymSz3ARMREVG+E0IgMjISJiYmWeYKuVFoEqKsBAcHo0yZMtJzfX19hISEZNh+yZIlmD9/fkGERkRERPns9evXqFy5cr6tv8gkRDKZDFpaWtLzhIQEqKurZ9h+5syZmDx5svQ8PDwcpqameP36NQwMDPI1ViIiouLo1atX2LdvX5bt+vbtCzMzszzZZkREBKpUqQJ9ff08WV9GikxCZGJigvDwcOl5yumzjGhqakJTUzPNdAMDAyZERERE36BOnTooV64cIiIiMmxjYGCAOnXq5Pnlrfzu7lIoOlVnR9u2beHn5yc9f/78Odq0aaPEiIiIiEoWFRUV2NvbZ9rG3t4+X/v65JdCE7EQQuFfAFixYgWePn0KAApjDyUlJcHLywvDhw8v+ECJiIhKMCsrKzg5OaW52mJgYAAnJydYWVkpKbLcKRSXzKKiorBnzx4AwK5duzBu3DgYGRlh//79MDc3h6WlJRo0aIDBgwdjypQpSEhIwKpVq1ChQgUlR05ERFTyWFlZwdLSEoGBgYiMjIS+vj5MTU2L5JmhFDKR+pRMMRYREQFDQ0OEh4dn2ocoOTkZiYmJBRgZUfGgrq4OVVVVZYdBRMVMdr+/c6tQnCEqDIQQePv2LcLCwpQdClGRVapUKVSoUIFjfRFRkcOE6P+kJEPlypWDjo4O/6AT5YAQAjExMXj//j0AoGLFikqOiKjwkMvlxerSUnHFhAhfLpOlJENly5ZVdjhERZK2tjaAL7UJy5Urx8tnRAB8fX3h5uamcJu6gYEB7O3ti2zn4+KKCREg9RnS0dFRciRERVvKZygxMZEJEZV4vr6+cHV1TTM9IiICrq6uRfqOrPyWlJSEoKAg+Pv7w9fXt0C2yYQoFV4mI8odfoaIvpDL5XBzc8u0jZubGywtLUvk5TMhBN68eQN/f3/pERAQIP3/9evXSE5OLtCYmBARERHlscDAwExHcwa+nCkKDAyEubl5wQRVwIQQ+PjxI/z8/PDs2TM8e/ZM+v/z588RHR2d6fIaGhowMzNDlSpVcOXKlXyPlwkRfZNTp07B29sb+/btw5IlS/DTTz8pO6Ri6f79+1i4cCHq16+P3377TdnhEFE2RUZG5mm7wiwmJgZPnjzBkydP0iQ/qUtufU1VVRVVqlSBhYUFLCwsYG5urvBvxYoVoaKiIt12n9+YEOWhgryTwMPDA87OzvD398ePP/6IKlWq4OnTpyhbtiymTZuGOnXq5Mt2gS8Vh7dt24YTJ06gZs2aOHbsWL4nRIcOHcLChQsRFhaGu3fvokyZMgCA6OhoHDhwAFOmTMHGjRthY2ODdevWYc2aNRgwYAAmTJgAa2traT1aWlro2LEjdHR0cPz4cdSrVw8WFhbw8fFBs2bNMHDgQEyePBn+/v7YsGEDevXqhbi4OOzfvx/jxo1Dq1atsHTpUtjY2GQY66hRo7Bq1Sqpk3FuVK9eHeHh4QV+6piIcie7hUjzu2BpXoqKisKTJ0/w+PFjPHr0SPo3ICAAGQ1pKJPJUKVKFdSsWRM1atRAjRo1pP9bWFhkWqS9oDEhyiMFfSeBra0tOnfuDA8PD6nTnhACS5cuhbW1NTw8PPDdd99luo7379/D3d0dvXv3ztG23dzcpGy9b9++6Nu377ftRA706tULderUQZ06dTB48GCcOHECAKCrq4uhQ4fi4sWL6NOnDwBg3LhxWLNmDWbOnJnmtV++fDnGjx8PADA3N8fQoUMxbNgwxMfHY8OGDbCzs0OXLl1w5coV9O/fH8CXQsGjRo3CwoUL0b9//0yToQ8fPsDFxQXff/89Bg4cmOv91tHR4S3sREWQqakpDAwMsiyCampqWoBRZU9CQgKePHmCu3fv4uHDh1LyExAQkOEyRkZGsLKyQs2aNRWSn2rVquXJj8OCwIQoDyjrTgI1NcXDJ5PJMHPmTFy/fh3Dhw/H/fv3M1w2ISEBgwYNynEyBHwZs0kZnQB1dHTw448/4ty5c1izZg0mTpwozdPS0pL+n/K6pPfLIyVp+pqmpqaUwKiqqqbbOVgmk2W53y4uLmjcuDG2bt2aJwlRynaJqGhJKYKa3ndDisJQBDUsLAz37t3D3bt3pcejR48yrNhQrlw51KlTB7Vr10bt2rWl/xsbGxdw5HmPCVEuFcY7CYYOHYqePXviwYMHsLKywqRJk2BsbIybN2/C1tYWU6dOxalTp/Do0SMcOnQIYWFhcHZ2xv/+9z8IIfD06VOYmJhg9erVada9Y8cOnDt3DhEREZgyZQqGDx+O33//HZqamoiNjcWdO3dw//59PH78GH/99Rf09PRw48YNLF++HNbW1nB1dcUff/yBOXPm4ODBgzh37hzmzZuH7777DmPGjMGbN29w+vRphctcqTVr1gw//vgjpk2bhhYtWqBJkyY5em2MjIy+aV52CCFw584drF69Gs2aNcOTJ09Qq1YtAEB8fDw2btyIP//8EydOnMCwYcPw+fNnXL58GdWrV0dwcDBmz56NGjVq4Ny5c5g7dy7at2+vsP7w8HD8/PPPOH/+PC5dugQ7Ozu8fPkSHTp0gIuLCxo2bIjly5dDR0cHmzZtwl9//QV7e3u8ffsWGzZskP7wubi4oFq1arhy5Qpu376NwMBAHDt2DMHBwbnafyJSlFIEtbCMQxQSEoJbt27hzp07uHv3Lu7du5fhWR9DQ0M0bNgQ9erVU0iAcvt3sjBjQpRLhfFOgpQv4WfPniEgIACPHz/G5cuXcfv2bbRp0wZTp07FTz/9hHXr1qFXr1745Zdf8ODBA/z999948eIFPn78CGNjY/z2229pBqocMmQIAgMDERAQgBUrViAxMRHa2trw9PTE0aNHcePGDcTHx8PJyQn//fcfDA0NcejQIXTq1AkvXrxAp06dMHz4cPj4+GDHjh04fvw4Ro4ciV27dsHHxweTJk3CypUrsW/fvgz3b9KkSbh69Sp+/vln+Pj45Ms1eH9/fzg7OytMy6xzIABcvnwZdnZ2+O6771CvXj1s27YNK1asAPDlbglbW1tMnjwZz549g7e3N7p164Zt27Zh6dKlWLlyJapVq4bZs2dDU1MTa9euTZMQGRoaYufOnTA3N5cKGxsbG6NLly5o2bIlNmzYgLp168LJyQktWrRAVFQUAGDChAnYtm0bDAwMMHHiRIwaNQoXL17EokWLcOHCBaipqaFUqVJ588IRkQJlFUENCwvD7du34eXlhVu3bsHLywshISHptjU3N0eDBg3QsGFD6WFmZlbizk4zIcqlwngnQcqbODk5Ge3atYOpqSmioqLwzz//SF+SX7O0tMS+ffuQkJAg3d4YFRWV5cjd6urqKF++PGxsbGBlZQUrKyscOXIEOjo6Uj+jnj17YsSIETh+/Dj69++PUqVKwc7ODtra2rCxsUFkZCR69OgBAKhfvz727t2b5T66uLigSZMmGDVqFP7+++9svzbZZWFhkeYM2eHDhzNdZt++fVi7di0AYNiwYVi4cCF+//13aGhoQCaTSUnHzz//DABo0KCB9Adq9OjR0NXVRVBQEB48eJDhcSpfvjz69euHDRs2YMOGDXB1dZXWZ2BggClTpkBDQwPdunXDx48fpV+EO3bsAPDlcmBKcUSZTIaffvoJq1atUrj8SER5S0VFJV9/ECckJODOnTu4efOmlPz4+fmlG0edOnXQuHFjNGrUCA0bNkT9+vVRunTpfIutKGFClEuF8U6ClA9C7dq1YWBggH/++Qdubm6ws7PLcBkNDQ28ePECZ86cwbBhwwAgw7sGviaTyRR+STx79kzh+rOqqirMzMykSzKp2349mrGKigrkcnmW20w589S8eXO0bds2W3Hmp3fv3uHRo0dYtGgRgC93v4WHh+P48eNwcnICkLYvkJqamrSvFStWxKJFi2BpaYnvv/8er169ynBbzs7OaNasGRYuXIibN29i6NChAID+/fvjxYsX+Pnnn2FtbY1Dhw7h9evXUFdXT3O2CwB27tyJX375BbVq1cKMGTOwYMGCvHgpiCifhYWF4b///sO1a9dw/fp13Lx5E3FxcWnaVa1aFU2bNoWNjQ2aNm2Kxo0bQ1dXVwkRFw1MiHKpMN5JsGvXLlhbW6NOnTrYsWMHbt26hd27d2d6h8CVK1fw119/wdPTM9fbNzU1hb+/PxISEqChoQHgS3KVcikvrzRq1Ahr1qzB+PHjs5UUJSYm4sGDB2jcuHGexgF86Vu1detW1K9fX5r24cMHbN26VUqIMjN69GjY2tril19+wc6dOzNtW7duXTRv3hxTpkxR6IPw6tUrzJs3DwMHDkTv3r0xffp0LFy4EM+fP8f9+/el2Ly8vGBjY4Pk5GRcvnwZBw8exODBg2Fra5tp0kxU0imjSKsQAoGBgVLyc+3aNTx8+DDND9ayZcuiWbNmsLGxgY2NDaytrYt1f5/8wIQol5R5J0FSUlKaaevXr4e7uzs8PDwAAHfu3MHnz5+RlJSEy5cvAwCeP38OU1NTaGho4PPnz3jy5Anu3LmDiIgIxMfH4+rVqwC+nPUwMDCQxvxJvd2v70BIPU5O9+7d8euvv+LQoUPo168fwsLCkJycDHt7ewBfPuCZnX3KaF5ycnKafR45ciSuXr2KoKCgNK/L1zH+/vvvGDt2rEK79F7DjKYnJCSkOx5QUlISvL29MXPmTIXpgwcPhoODA54+fQpLS0vpbJAQQjpblLKvd+7cQd26dRETE4Pr168jNjYW/v7+sLCwSPf1cnZ2Rrdu3fD69Wtp2sGDB9GvXz9UrVoV//vf/+Dq6gozMzM0b94c3bt3x/Lly6GqqoqHDx/CxsYGK1aswPr169G7d2+4urpm+4wgUUlUkEOrBAYGwt3dHVeuXIG7u7vC5zxF9erV0aJFC/zwww9o0aIFLC0tS1yfn7zGhCgPKONOAnd3d5w8eRLBwcEYM2YMDA0N4e/vDz09PXh7e8PCwgIAMGjQIHTp0gVNmzbFypUrUalSJezZswfz589H7969MWPGDBgYGMDR0RGbN29GvXr1sGrVKtSvXx9btmzBtm3bFLbr7e2NkydPIiwsDGfOnEH16tVx8eJF6RdMixYtoKOjg1OnTsHZ2Rl37txBXFwcDh06BE1NTZw5cwZv3rzB0aNHYWlpid27dwMAtm7dCgcHB5w4cQKPHz/Gv//+i++//17abmBgIDZs2AAPDw/8+OOPaNmypTTvr7/+wtSpUwEAAQEBWLdunbTvNWvWRHJyMl68eAEVFRXMnTsXUVFROHDgAN68eYNDhw6hSZMmaNq0KQDA09MTZ86cgZ+fH3bu3Ik+ffogLi4Of//9Nz58+ID9+/ejfv36Uvvw8HBMnjwZDx8+xPXr1/HDDz8A+JIk3bx5E0IIDB48GGvXrsXp06cBANu3b0fbtm1x9epVhISE4NGjR5gwYQKmTJmC69evY9CgQTh58iQePXqE+Ph43LhxA8+fP8eAAQNQvXp1AICDgwMcHR0VxiiKjY1F69atMXToUHz48EG6fLd//34MGTIEgwcPxo8//ojt27cDAE6fPo3Xr1+jZcuWqFu3bqG49EhUGOX30Crv3r2TEqArV67gxYsXCvPV1NTQuHFjKfn54YcfUL58+W/eHqVPJkrIz8KUob/Dw8OlTqUp4uLipF/jqcezySllnE6lksnT0xPv379Hr169lB2Kgrz6LBEVFnK5HGvWrMmyW8TEiROz/fc+NjYWnp6ecHNzw8WLF/H48WOF+aqqqrC2toadnR3atGmD77//vkT3/cns+zsv8QxRHsrvOwmIjh49iri4OJw8eRJ79uxRdjhExV5eDK0ihICfnx/OnTsHNzc3eHp6KnSClslkaNCgAezs7GBnZ4eWLVvm6xc/pY8JEVERcvjwYdy8eRP79+8vVDWAiIqrbx1aJTo6GpcuXYKbmxvc3NzS3NRSuXJlODg44Mcff4StrW2WQ5xQ/mNCRFSEZDZgJRHlvZwMrfLu3TucOnUKJ06cwKVLlxTOAmloaKBVq1awt7eHg4MDrKys2Am6kGFCRERElIGshlb58OEDAgICcP78edy4cUPhbk1zc3N06tQJDg4OsLW1LdH9gIoCJkREREQZ+HpoFSEE3rx5g0ePHuHJkyf49OmTQntra2t069YN3bp1Q926dXkWqAhhQkRERJSJWrVqoX79+ti4cSN8fHzw+fNnaZ6amhrs7OzQrVs3dO3aFZUrV1ZipJQbTIiIiIjS8fjxYxw8eBAHDx7E06dPpelaWlpo06YNBgwYgI4dO0p1G6loY0JERET0f169eoW9e/di//79ePTokTRdU1MTDg4O6N27Nzp37gw9PT0lRkn5gQkREREVa1kNmhsZGYkjR45g9+7dcHd3l6arq6vD3t4eTk5O6Nq1K8cGKuY4jHIx1rBhQ1y/fj3P1uft7Y2OHTti165daebdv38fvXr1wsKFC/Nse8oWFRWFhQsXSmU6MhMSEoIxY8Zg5MiR0rS8eP1jYmJgbm4Of39/PH36FH379s2T1/j69eto2LBhrtdDVNj5+vpizZo12LVrF44ePYpdu3ZhzZo1ePjwIS5cuID+/fujfPnyGDx4MNzd3SGTyWBnZ4cdO3bg/fv3OHnyJPr3789kqATgGaIi6sqVK5g8eTI+fvyIn3/+GQDw+fNn/P3333j37h1KlSqFWbNmoXbt2nm2TWNjYzx9+jTdIqDVq1dHeHh4usVPiyoNDQ2YmZnhw4cPWbZVV1dHWFgYNDQ0pGl58frr6Ohg7ty5MDExwefPnxEcHIyaNWvmap3Al/p7s2bNyvV6iAqz9GqQhYaG4sqVK5g7d67CYIo1a9bEoEGD0L9/f5iamhZ0qFQIMCEqouzs7NClSxe4u7tjxYoV0vTURU+dnJwUllm3bh3Gjx+f4fOsmJqaolKlSunO09HRUSg0WhxoaGigSpUq2WprbGyM6tWrIzAwUJr29eufnvv37yM0NBStW7fOsM3gwYMBABUrVoSZmVm24klP6uNdpkyZbMVHVFTJ5XK4ubkBAJKTk/H06VP4+Pjg+fPnUhsdHR0MGjQIv/zyC5o2bcpb5Es4JkRFmKqqappp/fv3T3f67t27cfToUekLccGCBXj58mWOt5lZ8cLi+MckJ/uU00K+ERERGDRoEFatWpVv20jxrcebqKgKDAxEQEAAfHx8cOfOHURHR0vzqlatiiZNmsDS0hJDhw5lDUoCwD5E6RJCIDo6WimP9C5HZdeWLVugrq4OFRUVuLq6wtraGh4eHvD394erqytevHiBGTNm4Pr167hw4QK8vb0xY8YMxMTE4O3bt/jtt98wfvx4tGrVCi9evADw5VfWzJkzMWfOHAwdOjRNPZ7MHDlyBNOmTcOYMWPQo0cPhUtP27Ztw+bNm/HLL79AJpOhQ4cOOH/+PFatWoXKlSvD29sbjRo1gpmZmcIvus2bN2POnDkYMGAAfvnlF8TExCAsLAyzZ8+GtbU1Nm3ahPLly+PatWvo3r07FixYgN9++w0VKlRA79698fbtW3Tq1AmGhobYsmWLtN7//e9/+O233+Dk5ARnZ+ds7+PevXsxYcIEODs74/Tp09L01K8/8OVM0KJFizBjxgzo6OggLi4Orq6uePnyJbZs2YJdu3Zh06ZNqFWrFg4ePAgTExPs2bMHmzZtgqWlpcLrHhcXBycnJxgaGmLChAmQy+W4desWLC0tMW/ePABfisBqamrCw8MDd+7cUTjeISEhmD17Npo0aSKtMygoCBMmTMCcOXNgZ2cn/bL29PSEnZ0dNm3ahEGDBsHQ0BDr1q3L9utDVNDkcjnOnj2Lfv36Ye3atbh27Rqio6Ohq6uLFi1aYMKECRg4cCDq1KkDNTW1bNcqoxJAlBDh4eECgAgPD08zLzY2Vjx+/FjExsYKIYSIiooSAJTyiIqKyvY+zZ07V1SoUEGMHDlS9O3bV5QrV06aFxUVJQwMDIS7u7sQQggXFxfRunVrhWUHDRokPe/Vq5f02kyYMEG0a9dOCCHE6tWrxezZs6XXqUyZMsLFxSXdeAYNGiTmzp0rhBDi3r17olWrVtK8MWPGCAcHByGEEI8fPxb169cXQggRHx8vypQpI/766y8hl8uFj4+PACD2798vkpOTRefOncX06dOFEEKcOXNG9O3bV1pnx44dxejRo0ViYqL4+++/hb6+vnB3dxe7du0S79+/F/b29qJr164iODhYBAYGCjU1NbFgwQIRGxsrjh07JipWrCiEEOL+/fuiatWqQgghPnz4IACIjx8/CiGEcHd3F2ZmZunur7e3t2jZsqX03MHBQXpNv379HR0dxbt374QQQixcuFB6r5mZmUltXr16JQCI1atXiwsXLoj79++LwMBAAUD4+/tLr7GdnZ14/PixcHNzE7q6umL37t1pXn8hhDA1NZXWnfp4JyUliX379insl42NjfDz8xNCCOHl5SW0tLTEq1evRFJSkmjQoIFwcnISnz9/FocPHxbly5dP9/UQIu1niaighIWFidWrV4vq1asr/E2tVq2acHJyEr/99puYN29emkfKZ4sKr8y+v/MSL5kVcdWqVcPmzZsBAKtXr5am6+rqolSpUtlaR0hICG7duoUdO3YA+HIpLuWOij/++AOHDx8G8GUwsnr16mVrnVu3blW4O2vo0KFo0qQJgoKCcP/+falSu4aGBqpWrYrY2FjIZDIp5pSO4g0aNEBISAgAYNOmTfjxxx8V1tm3b1+sW7cOFStWROnSpWFrayvNL1euHCwsLGBiYgIAqFChAlq2bAktLS3Ur18fb9++BQBYWlpi3759SEhIwJUrVwB8ucMsq+rTf/75JxwcHKTnNjY20pmcr19/AwMD9OzZExs2bMDo0aMVOl+nSOmv1LVrV1hYWABAumcMW7ZsCSsrK1hZWWHIkCE4d+4cBgwYkKZdRpf7VFVVUaFCBem5t7c3Xrx4gRo1agAAmjZtiho1amD37t2YM2cODA0N4eDggNKlS6NBgwZ4//59pq8LUUF6+vQp1q9fj507dyIqKgoAYGhoiKFDh0JbW1v6W5MeAwMDdqAmCROidOjo6EgfLGVs+1s5OjoqPM9u/5fXr19DXV09zaWiT58+ISQk5JsGIHv27BksLS2l51WrVgUABAcHo3nz5ggKCsK7d+9QunRpfPjwAV26dEk3ZjU1NcjlcmmddnZ2CuuMj4/Hx48fIZPJ0iz79fPUfatUVFSkZENDQwMvXrzAmTNnMGzYMADpJyJfu3//Ppo1a5bh/NTbX7FiBUaOHIlGjRph2LBh2LBhQ5r+QCntUy+X1TGsXbs2jh8/nmWsmcX27NkzJCYmKsyvWrUqgoOD07RVU1PL1WVdorwghIC7uzuWL18uXd4Fvtw9OWHCBPTv3x96enrp3mWWmr29/Tf3y6PihwlROmQyWZGsSvytNXQqVKiA58+f4/79+6hfvz4AwMvLCw0aNICKigp8fX2zfWYohampKZ48eSI9F0JAVVUV1atXR9myZTFv3jxs27YN+vr6OHLkiJQw5XSdZcuWhbGxMXx9fXMUX2pXrlzBX3/9BU9PzxwtZ2BgkO3thoeHw9XVFVeuXEHv3r3RrFkzDBky5FvCVRAWFiYlnioqKt807IGpqSkiIiLw5s0b6U5BIQRq1aqV6/iI8lJycjKOHj2KP/74A7dv3wbw5e91586dMWHCBLRt21YhgbeysoKTkxPc3NwUqtUbGBjA3t4eVlZWBb4PVHgxNS7CkpKSMv0CFEIonAUJDQ1FbGwsAgICoKGhgc+fP+Pjx4/Q09ND8+bN0b17dxw5cgTHjx/HhQsXoKmpiU6dOmHJkiUIDw9HTEwM3r17hw8fPiApKSnT7Y0YMQIeHh7SJSQvLy/07t0bZcuWxatXr3D8+HF06dIF7dq1Q+nSpaWzQCn/pj4LkfL/UaNG4ciRIwgPD5fWOXLkSOkX3tevRep4MnuN7ty5g4iICMTHx+PSpUsAgHfv3uHz58+ZruOnn37Cnj178PDhQwCAv7+/wmuTetnVq1cjKSkJdnZ2GDhwoMJx+fz5s0Kil3o/UtqljiEmJkZqd+7cOYwbNw4AUL58edy8eRNxcXE4duwYIiIi8P79eyQnJysc70+fPinE1rx5czRo0EC6ZJqcnIyXL1+ib9++AL4ck69fA54looIUFxeHv/76C7Vq1YKTkxNu374NbW1tjBs3Ds+ePcPJkyfRrl27dM+oWllZYeLEiRg0aBB69uyJQYMGYeLEiUyGKK187aFUiOSkU3VR4OHhIRo2bCh0dXXFjh07RFxcnML8c+fOCXV1dTF27FgRHh4ugoODhZmZmejWrZuIjY0Vd+/eFcbGxmLEiBFCiC8detu2bSv09fWFo6Oj9Dq9efNGdOjQQVSuXFmMGzdO2NrainHjxong4GCF7fn6+oqaNWuK77//Xjx79kwIIcTOnTtFy5YtxezZs8WkSZNERESEEEKId+/eiapVqwoDAwOhoqIiAAgzMzMRGBgo5s6dKwCIrVu3ipcvXwo7OztRq1Yt8fDhQyGEEEuXLhVt27YVc+bMEbNmzRIJCQkiNDRUDBs2TKioqIjt27cLIYR48OCBqFmzpmjVqpV4/vy5OHnypNDU1BTjxo0TQUFB4rfffhMAxJYtW0RAQICoXr26qFGjhjh9+rSoX7++GDJkiIiKihKjR48W6urq4vjx42mOQUJCghg7dqwwNjYWTk5OYuDAgaJr167i1q1baV7/1q1bi9atW4vly5eLX3/9VcTHxwshhJg3b54wMTERJ06cEOvWrRMAxIQJE6TXf/v27QKAWLBggUhISBBeXl7Czs5ODBs2TEybNk3cvHlTiicoKEhYWFiI6tWri0uXLolGjRqJmTNnitDQUIXjHRkZKcaMGSM0NDTEuXPnhBBCPH/+XLRp00aMGzdOTJgwQVy7dk0IIcSNGzdE2bJlhaOjowgICBBLliwRAKTX+WtF8bNEhVdoaKhYvHixKFeunNRJukyZMuJ///ufeP/+vbLDowJSUJ2qZUKUjJ96ERERMDQ0RHh4eJoh2OPi4uDv7w8LCwtoaWkpKcKS48yZM5DL5VK/oaSkJJw9exbJycno0aOHkqOj3OBnifJCWFgYVq9ejdWrV0tnhE1NTfHrr79i6NChRbJLA327zL6/8xIvmVGBmzVrFt6+fStdWoqNjYWPj4/CHWJEVPKEhoZi7ty5MDMzw/z58xEeHo66deti7969eP78OcaNG4cPHz7gwYMHCAgIkC6xE+UFdqqmArdixQrMnTsX06dPh5GRERo3bowlS5agdOnSyg6NiJQgNDRUOiOU0vm5bt26mDt3Lnr27Cnd3MHO0ZSfmBBRgWvfvj3at2+v7DCISMkiIiLw559/ZpoIAekXaU1Z3tXVFU5OTkyKKNeYEBERUYGKj4/H5s2bsWjRInz8+BFA+okQoFikNSNubm6wtLTkmEKUK3z3pFJC+pcT5Rt+higzcrkce/fuRa1ateDs7IyPHz+iZs2acHV1xb179+Do6JgmqQkMDFS4TJaeiIgIBAYG5mfoVALwDBEgDe0eExMDbW1tJUdDVHSljJGUWbkEKnmEEDh//jxmzJiBe/fuAQAqVqyI+fPnY/DgwVBTy/irKLvFV1mklXKLCRG+lHQoVaqUVKNJR0cn22UviOjLF15MTAzev3+PUqVKKZRJoZLN29sbU6dOhbu7O4AvdcamT5+OiRMnZqtUkb6+fra2k912RBlhQvR/UopdsnAl0bcrVaqUQuFYKrnevn2LWbNmYefOnRBCQENDA+PHj8fMmTOzLJycmqmpKQwMDDK9bMYirZQXmBD9H5lMhooVK6JcuXJpCl0SUdbU1dV5ZogQHx+PNWvWYNGiRdJlrH79+mHx4sUwMzPL8fpUVFRgb2/PIq2U75gQfUVVVZV/1ImIckgIgRMnTmDKlCl48eIFAMDGxgZr1qxBs2bNcrVuFmmlgsCEiIiIcsXPzw/jxo3DxYsXAXzpML106VL0798/z87cWFlZwdLSEoGBgYiMjIS+vj5MTU15ZojyDBMiIiL6JrGxsViyZAmWLVuGhIQEaGpqYsqUKZgxYwb09PTyfHsqKiowNzfP8/USAUyIiIjoG5w7dw7jxo3Dy5cvAXzpx7N+/XpUq1ZNyZERfRsmRERElG1BQUFwdnbGkSNHAACVKlXCmjVr0LNnzzwbrkQul/PSGBU4JkRERJQluVyOzZs3Y/r06YiKioKqqiomTpyIefPm5ekYQCziSsrClJuIiDLl5+cHW1tbjB07FlFRUWjevDl8fHzw559/5nky5OrqmmbMoZQirr6+vnm2LaKvMSEiIqJ0JSUl4Y8//kCDBg1w9epV6OrqYu3atbh27Rrq16+fp9vKbhFXuVyep9slSsFLZkRElMa9e/cwZMgQ+Pj4AADat2+PLVu25NtdXjkp4so7zSg/8AwRERFJkpKSsGjRIlhbW8PHxwelSpWCi4sLzp8/n6+JCIu4krLxDBEREQH40ldo4MCBuHnzJgCgR48e2LhxY4HUp2MRV1I2niEiIirhhBDYuHEjGjZsiJs3b8LQ0BB79uzBkSNHCqxYb0oR18ywiCvlJyZEREQlWHBwMOzt7TF27FjExsaibdu2ePDgAfr3759n4wplR0oR18ywiCvlJ76ziIhKqMOHD6NevXq4cOECtLS0sGbNGly4cAFVqlRRSjwpRVy/PlNkYGAAJycnjkNE+Yp9iIiISpiYmBhMmjQJW7ZsAQBYW1tjz549qFWrlpIjYxFXUp5CkRBFR0dj2rRpKF26NKKiorBs2TJoamqmaff7779DLpdDRUUFMTExWLhwYYGe0iUiKuoePnyIn3/+GY8ePYJMJsP06dOxYMECqKurKzs0CYu4kjIUipR79OjRaNeuHRYtWoTGjRtj5syZadqcPn0a/v7+mDNnDmbNmoWQkBC4uroqIVoioqJHCIHNmzejadOmePToESpUqIALFy5gyZIlhSoZIlIWpSdEISEhOHToEBwcHAAADg4O2Lx5c5qxJnx9fRWmaWtrIzw8vEBjJSIqikJDQ9GrVy+MHj0acXFxcHBwwL1799CuXbtcr1sulyMgIAAPHjxAQEAAR5KmIkvpl8w8PDxgZGQELS0tAICxsTE0NDTg5eWFtm3bSu26du2KuXPn4vDhw2jbti0+fvyIAQMGZLje+Ph4xMfHS8+zGgGViKg48vHxgaOjI/z9/aGuro6lS5fC2dk5T/rksBArFSdKP0MUHByMMmXKKEzT19dHSEiIwjRLS0vs27cPAwcORN++feHi4gJtbe0M17tkyRIYGhpKD2XdNUFEpCzbt2/H999/D39/f1hYWODff//F5MmT8ywZYiFWKk6UnhDJZDLp7FCKhISEdK9px8XF4eTJk3jx4gX69OmDpKSkDNc7c+ZMhIeHS4/Xr1/neexERIVRbGwshg4dimHDhiE+Ph6dO3eGt7c3rK2t82T9LMRKxZHSEyITE5M0fYGioqJgYmKiMM3DwwO3b99Gu3bt4OHhgTt37mDlypUZrldTUxMGBgYKDyKi4u7ly5f44YcfsGPHDqioqGDx4sU4ceIESpcunWfbyEkhVqKiQukJka2tLYKCgpCQkAAA0qUyGxsbhXaHDx9GjRo1AHxJohYtWoR//vmnYIMlIirEzpw5gyZNmuDOnTswMjLC+fPnMWvWrDwfw4eFWKk4UnpCZGJiAnt7e3h6egIALly4gDFjxkBLSwsrVqzA06dPAQANGzbEnTt3pOVkMlmapImIqCQSQmDp0qXo0qULwsLC8N1338HHxydP7iJLDwuxUnGk9LvMAGDz5s2YMWMGbt68ic+fP2Pp0qUAgP3798Pc3ByWlpYYMmQIAgMDsXTpUhgZGeH169eYNWuWkiMnIlKulP5C+/fvBwCMHDkSa9euhYaGRr5tM6UQa2aXzViIlYoamRBCKDuIghAREQFDQ0OEh4ezPxERFQtBQUHo3r07vL29oaamhrVr12L06NEFsu2Uu8wywtpjlFcK6vtb6ZfMiIgo5/777z9YW1vD29sbZcuWxcWLFwssGQJYiJWKn0JxyYyIiLJv165dGDFiBBISElC3bl2cPHkSFhYWBR4HC7FSccKEiIioiBBC4LfffsPixYsBAN27d8fu3buV2nmZhVipuGBCRERUBMTHx2Pw4MFS5+nZs2djwYIFGZ6NkcvlPHNDlANMiIiICrlPnz6hR48euHr1KtTU1LBlyxYMHjw4w/asMUaUc/y5QERUiL148QLff/89rl69CgMDA5w7dy7LZIg1xohyjgkREVEh9d9//6FZs2bw8/ODqakprl+/nulgi6wxRvTtmBARERVCJ06cgJ2dHT5+/IjGjRvjxo0bqFu3bqbLsMYY0bdjQkREVMjs2LEDPXv2RFxcHDp37gxPT09UrFgxy+VYY4zo2zEhIiIqJIQQ+OOPPzB06FDI5XIMHjwYx44dg56eXraWZ40xom/HhIiIqBCQy+WYOnUqpk+fDgCYPn06tm/fDjW17N8MnFJjLDOsMUaUPiZERERKlpiYiMGDB+PPP/8EAKxYsQJLly6FTCbL0XpUVFRgb2+faRt7e3uOR0SUDn4qiIiUKCYmBj169MDu3buhqqqKXbt24ddff/3m9bHGGNG34cCMRERKEhkZic6dO+Off/6BlpYWDh06hM6dO+d6vawxRpRzTIiIiJQgLCwMDg4OuHHjBgwMDHDmzBm0aNEiz9bPGmNEOcOEiIiogH369AkdOnSAj48PSpcujQsXLsDa2lrZYRGVaEyIiIgK0Lt379CuXTs8fPgQxsbGuHTpEurXr5+rdbKQK1HuMSEiIiogwcHBaNu2LZ4+fYqKFSvi8uXLue7kzEKuRHmDPyGIiArAq1ev0KpVKzx9+hRVqlTBP//8kyfJEAu5EuUNJkRERPns1atXsLW1xcuXL1G1alX8888/qF69eq7WyUKuRHmLCRERUT56/fo17OzsEBAQgOrVq8PT0zNP7v5iIVeivMWEiIgonwQHB8POzk46M+Tu7o7KlSvnybpZyJUobzEhIiLKB2/evIGdnR2eP38Oc3PzPE2GABZyJcprTIiIiPLYu3fvYGdnBz8/P5iamsLd3T3PC6qykCtR3mJCRESUh96/fw87Ozs8efIElStXhru7e76MGM1CrkR5i58UIqI8Ehoaivbt2+Px48cwMTGBu7s7qlatmm/bYyFXorzDgRmJiPJAVFQUOnbsiPv376NChQpwd3fP9a312cFCrkR5gwkREVEuxcfHo3v37rhx4wZKly6NixcvombNmgW2fRZyJco9/oQgIsqFpKQk9OnTB5cvX4auri7OnTuHunXrKjssIsohniEiIvpGcrkcw4YNw7Fjx6CpqYmTJ0/iu+++y5P18hIYUcFiQkRE9A2EEJg0aRJ27doFVVVVHDx4EHZ2drleL4u1EikHf3IQEX2D+fPnY+3atQAAFxcXdOvWLdfrZLFWIuVhQkRElEMbN27E/PnzAQDr1q3DgAEDcr1OFmslUi4mREREOXDs2DGMGzcOADBv3jzp/7nFYq1EysWEiIgom65fv46+fftCCIERI0bgf//7X56tm8VaiZSLCRERUTb4+vqiS5cuiIuLQ5cuXbBhwwbIZLI8Wz+LtRIpFxMiIqIshISEwN7eHqGhofjuu+9w4MABqKnl7U26LNZKpFxMiIiIMhEREYGOHTsiMDAQNWrUwOnTp6Gjo5Pn22GxViLl4ieLiCgDCQkJ6NmzJ+7du4fy5cvDzc0NRkZG+bY9FmslUh4OzEhElI6UjtOXL1+Gnp4ezp49m6+V61OwWCuRcjAhIiJKx9KlS6VRqA8dOoTGjRsX2LZZrJWo4PEnBxHRVw4dOoRZs2YB+DLwYlZ9e4io6OMZIiKiVG7evImBAwcCACZOnIjRo0fn2bpZtJWo8GJCRET0f169eoVu3bohLi4OnTp1wp9//pln62bRVqLCjT9NiIjw5fb6zp074927d2jQoAH2798PVVXVPFk3i7YSFX5MiIioxEtKSkLv3r3x8OFDVKxYEadOncqzEaFZtJWoaGBCREQl3q+//go3Nzdoa2vj1KlTqFKlSp6tm0VbiYoGJkREVKLt2LEDa9euBQDs3bsXTZo0ydP1s2grUdHAhIiISqz//vtPuots/vz56NmzZ55vg0VbiYoGJkREVCIFBwejZ8+eUnmOOXPm5Mt2WLSVqGhgQkREJU5cXBx69OiBt2/fom7duti1a1e+jQfEoq1ERQM/gURUoqTUKLt16xbKlCmDEydOQE9PL1+3yaKtRIUfB2YkohJl9erV2LNnD1RVVeHq6logBVsBFm0lKuyYEBFRiXHx4kVMmTIFALBy5Uq0bdu2QLfPoq1EhRcTIiIqNjKrFRYQEIDevXtDLpfjl19+wfjx45UcLREVJkyIiKhYyKxWmIWFBRwdHREaGgpra2ts2rQJMplMidESUWHDhIiIiryUWmFfS6kV9vDhQ3h7e6Ns2bI4cuQItLS0lBAlERVmTIiIqEjLqlaYj48PTp48CZlMhn379nG8HyJKF29vIKIiLbNaYSEhIThz5gwAYNKkSejQoUNBhkZERQgTIiIq0jKqARYTEwNXV1ckJyejZs2aGDhwYAFHRkRFCRMiIirS0qsBJpfLcezYMYSFhaF06dLo0aMHDA0NlRAdERUVTIiIqEhLr1bYP//8g2fPnkFNTQ1OTk4oX748+w4RUaaYEBFRkfZ1rbAXL17Aw8MDANC5c2dUrFiRtcKIKEuF4i6z6OhoTJs2DaVLl0ZUVBSWLVsGTU3NdNvGxMTgr7/+QtmyZVGzZk00a9asgKMlosImpVaYq6srjh49CgBo3LgxWrVqBXt7e9YKI6IsFYqEaPTo0ejRowd69OiB3bt3Y+bMmVi5cmWadp8/f8aAAQOwfv16WFhYKCFSIiqsatasCU9PT0RHR6NGjRr4+++/UbNmTZ4ZIqJskQkhhDIDCAkJQbVq1RAaGgotLS18+PABZmZmePfuXZrOkh06dMCsWbNga2ub4+1ERETA0NAQ4eHhafobEFHRN3/+fMybNw+6urq4ffs2atWqpeyQiCgPFNT3t9J/Onl4eMDIyEgaOdbY2BgaGhrw8vJSaHf69Gk8e/YMXl5e6NixI2bOnInExMQM1xsfH4+IiAiFBxEVT1euXMH8+fMBAJs3b2YyREQ5pvRLZsHBwShTpozCNH19fYSEhChM+/vvv/Hdd9/B2dkZQ4cORaNGjSCEwNKlS9Nd75IlS6Q/kERUNGRWnDUj7969Q79+/SCEwJAhQ9C/f/8CipaIihOlJ0QymSxNXaGEhASoq6srTHv06BFGjBgBDQ0NlC1bFsOHD8emTZsyTIhmzpyJyZMnS88jIiJQpUqVvN8BIsoTmRVnzahTdHJyMvr374+3b9+iTp06WLduXUGFS0TFjNIvmZmYmCA8PFxhWlRUFExMTBSmJSUlITk5WXpev359fP78OcP1ampqwsDAQOFBRIVTSnHWry9tpxRn9fX1TXe5JUuW4NKlS9DR0cGhQ4ego6NTEOESUTGk9ITI1tYWQUFBSEhIAADpUpmNjY1Cu/r16+PZs2fSczU1NfYTICoGsirOCgBubm6Qy+UK065evYq5c+cCADZu3Mhb64koV5SeEJmYmMDe3h6enp4AgAsXLmDMmDHQ0tLCihUr8PTpUwDAxIkTcebMGcTHxwMA/v33X4wfP15pcRNR3sisOGuKiIgIBAYGSs/DwsLQv39/yOVyDBw4EIMGDcrvMImomFN6HyLgy10hM2bMwM2bN/H582epX9D+/fthbm4OS0tLNG/eHAsXLsT48eNhaWkJLS0tDB06VMmRE1FuZVScNaN2QgiMGjUKgYGBqFatGtavX5+f4RFRCVEoEiIjIyNs27YtzXRvb2+F5/379+cdJETFTHrFWTNrt3v3bhw8eBBqamrYt29ftpcnIsqM0i+ZEVHJll5x1q8ZGBjA1NQUz58/x9ixYwF8GYjx676GRETfigkRESnV18VZ02Nvb4/k5GT07dsX0dHRaN26NaZPn15AERJRScCEiIiULqU469dnigwMDODk5AQrKyvMmzcPt27dQunSpbFnzx6oqqoqKVoiKo4KRR8iIiIrKytYWlqmO1K1h4cHlixZAgDYsmULB1klojzHhIiICg0VFRWYm5srTPv8+TMGDBggleZwdHRUTnBEVKzxkhkRFVpCCIwcORJBQUGoUaMG1qxZo+yQiKiY4hkiIsp331K0FQBcXFxw+PBh6RZ7PT29AoiWiEoiJkRElK++pWgrAPj7+2PixIkAgEWLFsHa2jrfYyWikouXzIgo33xr0dbk5GQMGjQIUVFRaNmyJaZMmVIQ4RJRCcaEKJtatWoFmUwGmUwGdXV1WFpaYt++fQpt7t69i59//hkVKlSAhoYGqlWrhnnz5iExMTHDdaV+9OvXL09j3rhxIywsLKClpYUmTZrg6tWr2V52yZIlkMlkcHZ2VpgeGRkJZ2dnmJmZQVtbG99//z1u3bqVp3FT8fCtRVsBYPXq1bh69Sr09PSwc+dO3mJPRPmOCVE2CCFw9+5dLFmyBG/evIGfnx9atGiBQYMGwd/fH8CXvg42NjYoX748Tp8+DV9fX/z2229Yu3YtfvnllzTrWrFiBd68eaPw+Ouvv/Is5oMHD8LZ2RmzZ8/GnTt30LJlSzg4OCgUyMzIrVu3sGXLFtSvXz/NvGHDhuHixYvYs2cPHjx4gA4dOqBdu3YIDg7Os9ipePiWoq0A8OjRI8yaNQsAsGrVKlStWjXfYiQiSiETQghlB1EQIiIiYGhoiPDw8CzLBHzNz88PlpaWuHr1Klq0aAEAePHiBapXr47jx4/D0NAQbdu2xY4dO9JU3d6+fTuGDRsGX19f1KpVS1qXl5cXmjZtmmf797XvvvsOjRs3xqZNm6RpVlZW6N69uzSeS3qioqLQuHFjbNy4EYsWLULDhg2xevVqAEBsbCz09fVx4sQJdOrUSVqmYcOG6Ny5MxYtWpRv+0NFz4MHD3D06NEs2/Xs2RP16tUDACQkJKBZs2a4c+cOOnXqhFOnTkEmk+V3qERUiOXm+zsneIYoG7y9vSGTyRTOmAQFBQEAypcvj4kTJ8LBwSFNMgQAbdq0AQDcu3dPWpeamlq6Z1++9vvvv0NPTy/TR3qXwRISEuDt7Y0OHTooTO/QoQP+/fffTLc5duxYdOrUCe3atUszLykpCcnJydDS0lKYrq2tjWvXrmW5P1Sy5LRoK/Cl8/SdO3dQpkwZbN26lckQERUY3mWWDT4+PrCwsJAy06dPn2LKlClo2LAh1NXVcf/+fSxdujTdZWNjYwEAGhoa0rqSk5NRtmxZhXZ9+vTB1q1bFaaNGjUKTk5OmcZWqVKlNNM+fvyI5ORklC9fXmF6+fLl8fbt2wzXdeDAAfj4+GTYJ0hfXx/NmzfHwoULYWVlhfLly2P//v24efMmatSokWmcVPKkFG3N7LJZStFWAPDy8sLvv/8OANi8eTMqVqxYIHESEQFMiLLF29sbAQEB0NPTQ1JSEmQyGXr37o1ly5bhzJkzAL5cNkqPj4+Pwnxvb2/06tULixcvVmhXunTpNMuWKVMGZcqU+ea4v/51LYTI8Bf369evMXHiRFy4cCHNGaDU9uzZgyFDhqBSpUpQVVVF48aN0bdvX2k/iVKkFG11dXXNsI29vT1UVFQQExODgQMHIjk5GX369EGvXr0KMFIiIl4yy5Y7d+5g6tSpuHv3Ll6+fImYmBjs3LkT5cuXl+4gyyiJ2LBhA1q2bAkLCwtpXS1atED16tUVHl+fMQK+/ZKZkZERVFVV05wNev/+fZqzRim8vb3x/v17NGnSBGpqalBTU4OnpyfWrl0LNTU1JCcnAwCqVasGT09PREVF4fXr1/Dy8kJiYqK0f0SpZadoKwDMnDkTT58+hYmJCdavX6+MUImohMv2GaKePXtCX18f/fv3R/v27fMzpkLl5cuXCAsLQ/v27VG9evU08xs1agQA8PT0RPfu3RXm/fnnn7h79y7+++8/hXWlLJOVb71kpqGhgSZNmuDixYvo0aOHNP3ixYvo1q1buutp27YtHjx4oDBt8ODBqFWrFqZPn57mtmddXV3o6uoiNDQU58+fxx9//JGtfaKSJ7OirQBw+fJlrF27FgCwY8eOXJ0VJSL6VtlOiDw9PXH79m3pTMDJkyeleTo6Oul2wi0OUjpUN2nSJN35NjY2sLe3x9ixY5GYmAhra2u8e/cO27Ztw969e3Hw4EEpAfL29gaQfl+ecuXKpSllkJtLZpMnT8aAAQNgbW2N5s2bY8uWLQgMDMSoUaMAAOvXr8exY8dw+fJlAF/6B9WtW1dhHbq6uihbtqzC9PPnz0MIAUtLSzx//hxTp06FpaUlBg8e/E1xUsmQXtFWAAgPD5feO6NGjcKPP/5YwJEREX2R7YTI1tZW4bJI1apV8f333+PAgQPSnVTFkY+PD6pXr45SpUpl2ObIkSOYO3cufv31V7x9+xZGRkaws7PD3bt3UatWLYV1AUDNmjUVlldXV0dkZCQ0NTXzLO7evXvj06dPWLBgAd68eYO6devi7NmzMDMzA/Cl4/WLFy9yvN7w8HDMnDkTQUFBKFOmDH766ScsXrwY6urqeRY7lRxTp07F69evUa1aNSxfvlzZ4RBRCZbtcYgGDhyI3bt3K0zr3LkzTp8+naZtXFxcph1zlaGgxjEgKom+pXjrpUuXpMvvnp6eaNWqVUGESkRFTEF9f+fqLrOUW8m/tn//fl5CISohvqV4a2RkJIYNGwYAGDduHJMhIlK6bJ8hql+/Pnr06IHUzQ8dOpTm9ti4uDgcOHAgWyUiChLPEBHlvZTirRlJfSdZamPHjpVq7d2/fx96enr5GSYRFWEF9f2d7YQoq9PfCiuVyaTbtAsLJkREeUsul2PNmjVZDrw4ceJEhb8f7u7usLOzA/DlDrOU/xMRpafQle4YMWIEYmNjIZfLM31ERUVh3Lhx+RYwERUO31K8NTo6WrpUNnLkSCZDRFRoZLsP0ZgxY7J1F5SOjg6GDh2aq6DyU3R0dJoxdQBAVVVVoSN4dHR0hutQUVGBtrb2N7WNiYlBRiflZDIZdHR0vqltSrKaEV1d3W9qGxcXl+nZvpy01dHRkUbKjo+PR1JSUp601dbWls5AJCQkSINl5ratlpaW9F7JSdvExEQkJCRk2FZTUxNqamo5bpuUlIT4+PgM22poaEh3++WkbXJyMuLi4jJsq66uLvUXTN323bt3aWJXUVGR4pXL5UhKSsK7d+9gbGwMAJg1axZevnyJKlWqcOwqIipcxDcICwsThw8fFitXrhR//vmnOHbsmAgLC/uWVRWY8PBwASDDR8eOHRXa6+joZNi2devWCm2NjIwybGttba3Q1szMLMO2tWvXVmhbu3btDNuamZkptLW2ts6wrZGRkULb1q1bZ9hWR0dHoW3Hjh0zfd1Sc3R0zLRtVFSU1HbQoEGZtn3//r3UdsyYMZm29ff3l9pOmTIl07YPHz6U2s6dOzfTtl5eXlLbP/74I9O27u7uUtv169dn2vb06dNSWxcXl0zburq6Sm1dXV0zbevi4iK1PX36dKZt169fL7V1d3fPtO0ff/whtfXy8sq0bevWrcW8efPEvHnz0j1uMplMABBubm6CiCg7Ur6/w8PD83U7ObrLLD4+HlOnTsW2bdsQHx8vnb2QyWTQ1taGs7Mz5s2bJ/1CJCJKTQiBAQMGcABGIip0st2pWi6Xo02bNrh69Srq168PBwcHmJmZQSaTITAwEOfPn4ePjw969eqFgwcP5nfcOZbSKSskJCTdTlm8ZJZ+W14y4yWzjC6ZAcCTJ09w9OhR6fnXl8y6du2KWrVqYcaMGVi/fj3Kli2LJ0+ewMjIKMPtERGlVujuMtu4cSOmTp2KHTt2oHfv3um2OXnyJAYOHAgXFxeFGlqFAe8yI8ofWY1DdOPGDfzwww+Qy+U4c+YMOnbsqMRoiaioKXQJkZ2dHYYPH44+ffpk2u7cuXPYvHkzTpw4kScB5hUmRET5J6ORquPj49G4cWM8fvwYAwcOxK5du5QdKhEVMYXutvu4uLgskyEAcHBwyPRyDBEVPynFW+vVqwdzc3PpsuSyZcvw+PFjGBsbo2nTpti5c2emlx6JiJQl272fK1eunO2VVqxY8ZuCIaKiITu1y3x9fbF48WIAwPLly/HLL78AAHr16sViwERU6GQ7IXr//n22V5rVYG1EVHRlp3aZXC7H8OHDkZCQgE6dOuGnn36SEiIiosIo2wnRvXv3MHfu3Azvekrt2rVruQqKiAqnjGqXRUREwNXVVapdtmXLFly/fh16enrYuHGjdMcgEVFhle2EKDw8HAsXLsxWW/7xIyp+5HI53NzcMm3j5uYGPT09TJs2DQDw+++/w9TUNNOhKYiICoNsd6oeP348oqKislXLbMyYMfkZMxEpQXZql4WHh2Po0KGIjIxEs2bN+LeAiIqMbCdEP/30EwICAvDvv//C3d09zfzk5GRcv34dOjo6GDFiRJ4GSUTKFxkZmWUbX19fXLx4Eerq6ti6dWu6dQOJiAqjbCdEZ86cQYMGDXD27Fmp42Rqqqqq0NbWxqZNm1CvXr08DZKIlE9fXz/T+bGxsTh79iwAYMaMGahbt25BhEVElCey3YfI1NQUu3btQt++fTNs07hxY4SHh+Po0aPo2bNnngRIRIWDqakpDAwMMrxsdvHiRURFRcHS0hKzZs1SmKepqSl1xtbU1Mz3WImIcirbZ4h8fHwyTYZStGnTBidPnsxVUERU+KioqMDe3j7deQEBAfDx8QEAbN26VaEuIACoqamhV69e6NWrF4s/E1GhlO2EKCdyMmYRERUdVlZWcHJyUhg+PzExEWfOnAEAjBw5Ei1btlRWeERE3yzbP9Wye9usXC7Hs2fPvjkgIircrKysYGlpKY1UvX37dnz48AEVK1bE0qVL010mKSkJx44dAwD06NGDZ4mIqNDJ9hkiIQTu3LmTZbv9+/fDyMgoV0ERUeGWUrtMQ0MDmzZtAgCsWbMGpUqVSrd9fHw8nJyc4OTkhPj4+AKMlIgoe7KdEE2aNAldunSR7iJJz4EDBzBq1CgMHTo0T4IjosJLCIExY8YgISEB9vb2cHR0VHZIRETfLNvnrZs1a4YpU6agS5cuMDU1RcuWLVGhQgWoq6vj7du38PT0hL+/P7p3745hw4blZ8xEpCSpi7peunQJV65cgZaWFjZs2MAR6omoSMvRhXxnZ2fUrl0b06ZNw969exXmlS5dGgsXLsSMGTPyNEAiKhxSF3WNiYnB+vXrAQCjRo1C1apVlRwdEVHu5LhnY4cOHdChQwc8e/YMT58+RVRUFKpUqYKmTZtCQ0MjP2IkIiX7uqjr5cuXERMTA2NjY+jp6cHX1zfdAVuJiIqKb77Vo0aNGqhRo0ZexkJEhdDXRV0DAwPh7e0NAOjUqRPU1NTg5uYGS0tLqKjky0geRET5jn+9iChTqYu6Jicn4/Tp0wCAhg0bwtzcHAAQERGBwMBAZYVIRJRrHAyEiDKVuqjrjRs38P79e2hra6N9+/YZtvuahoYGXFxcpP8TERU2TIiIKFMpRV3DwsLg4eEB4EtfQl1d3XTbpUddXR2//PJLfoVIRJRrvGRGRJkyNTWFvr4+zp49i8TERJiZmaFhw4YKbQwMDGBqaqqcAImI8gDPEBFRplRUVKCurg4/Pz+oqKigc+fOacYcsre3z7RDdVJSEs6fPw8A+PHHH1m6g4gKHf5VIqJMRUZGYtmyZQAAOzs7GBsbS/MMDAxgb2+f5S338fHx6Ny5MwAgKiqKCRERFTr8q0REmZo7dy6Cg4NRtWpVHD9+HB8+fEBkZCT09fVhamrKW+2JqFhgQkREGXrw4AHWrl0LAFi/fj10dXXTdKYmIioO+NOOiNIlhMC4ceOQnJyMHj16wMHBQdkhERHlG54hIqJ0HThwAP/88w+0tbXx559/IiAggJfKiKjYYkJERGlERkZiypQpAIDhw4fj+PHj0mjVQPY7UxMRFRX8iUdEaSxYsAAhISGoUqUK9PX1FZIh4EupDldXV/j6+iopQiKivMUzRESk4PHjx1i9ejWAL2MGqaurZ9g2u0VdNTQ0sH79eun/RESFDRMiIpIIITBhwgQkJSWhbdu2qFy5cqbtU4q6phR5zYi6ujrGjh2bh5ESEeUtXjIjIsnhw4dx+fJlaGpq4tdff83WMpkVdSUiKioKRUIUHR2NsWPHYs6cOXB2dkZ8fHym7Q8cOABbW9uCCY6ohIiKisLkyZMBANOnT892h+nMirqmSE5OhoeHBzw8PJCcnJyrOImI8kOhSIhGjx6Ndu3aYdGiRWjcuDFmzpyZYduQkBAsXry4AKMjKhl+//13BAUFwdzcHDNmzICpqSkMDAwyXSa7RV3j4uLQpk0btGnTBnFxcXkVMhFRnlF6QhQSEoJDhw5Jg745ODhg8+bNGZ6GX7p0KUaPHl2QIRIVe35+flixYgUAYPXq1dDW1oaKigrs7e0zXS6roq5EREWF0v+SeXh4wMjICFpaWgAAY2NjaGhowMvLK03brVu3on///tDR0clyvfHx8YiIiFB4EFFaKR2pExMTYW9vj65du0rzrKys4OTklOZMkYGBAZycnDgOEREVG0q/yyw4OBhlypRRmKavr4+QkBCFac+fP0dERARsbGzw+PHjLNe7ZMkSzJ8/P09jJSqOjh8/jvPnz0NDQwNr166FTCZTmG9lZQVLS0sEBgZypGoiKraU/hdNJpNJZ4dSJCQkKIx9kpycjI0bN8LZ2Tnb6505cybCw8Olx+vXr/MqZKJiIyYmBpMmTQIATJkyBTVq1Ei3nYqKCszNzVGvXj2Ym5szGSKiYkfpZ4hMTEwQHh6uMC0qKgomJibS83///RebN2/Gjh07AHxJmBISElCqVCmEhYWlu15NTU1oamrmW9xExcHSpUvx6tUrVKlSBbNmzVJ2OERESqP0hMjW1hYjRoxAQkICNDQ0pEtlNjY2UpumTZsqXCY7fPgwDh8+jAMHDhR4vERFnVwuR2BgIB4/foxly5YBAFauXAldXV0lR0ZEpDxKT4hMTExgb28PT09PtG/fHhcuXMCYMWOgpaWFFStWoEuXLrC0tFQYCTelE3ZWo+MSkSJfX1+4ubkhIiIC+/btQ0JCAmrWrInatWvn63bV1dXxxx9/SP8nIipslJ4QAcDmzZsxY8YM3Lx5E58/f8bSpUsBAPv374e5uTksLS2VHCFR0efr6wtXV1cAwNOnT+Hn5wcVFRW0b98ehw4dyte7xjQ0NDB16tR8WTcRUV6QCSGEsoMoCBERETA0NER4eHiWg80RFTdyuRxr1qxBREQEEhMTsXHjRoSGhuL7779Hhw4dAHy5lX7ixInsME1EhUpBfX/zLx9RCRAYGCiNxXX9+nWEhoZCX18frVu3ltqkFGrND8nJybh16xZu3brF0h1EVCgViktmRJS/UkZ+Dw0NxbVr1wAAHTp0SHMnZn4Vao2Li5NulIiKimIHbiIqdHiGiKgESCnAev78eSQlJcHc3Bx169bNsB0RUUnDhIioBDA1NUVwcDCePHkCmUwGBweHNCNSZ7dQKxFRccSEiKgESExMxIULFwAA3333HcqXL5+mDQu1ElFJxj5ERCXAypUr8erVK5QtWxZdu3ZFQkKCNM/AwAD29vYs1EpEJRoTIqJiLjAwEIsWLQIArFq1Cv369WOhViKirzAhIirmfv31V8TExKBFixbo378/ZDIZR3knIvoKEyKiYiClPtnXZ30uXbqEw4cPQ0VFBevXr0/TkbqgqKurY+7cudL/iYgKGyZEREVc6vpkKQwMDGBnZ4fx48cDAMaOHYsGDRooK0RoaGhg3rx5Sts+EVFWmBARFWGp65OlFhERgalTp+LJkycwNjbGggULlBAdEVHRwYSIqIiSy+Vwc3NLd15ERAQ8PT0BAEuWLEGpUqUKMLK05HI5fH19AQBWVlbsxE1EhQ4TIqIiKnV9sq9duHABCQkJqFy5Mtq0aVPAkaUVGxsrjYzN0h1EVBjxZxpREZVR3TF/f388fPgQANCxY0dER0cXZFhEREUSEyKiIiq9umPJyck4e/YsAMDa2homJiasT0ZElA1MiIiKKFNTUxgYGChM8/LywocPH6CtrQ07OzvWJyMiyiYmRERFlIqKCuzt7aXnkZGR8PDwAAC0a9cOOjo6rE9GRJRN/EtJVIRZWVnByckJBgYGuHjxIuLj42FiYoLWrVvDycmJ9cmIiLKJd5kRFXFWVlb48OEDJk+eDJlMhg0bNqBr1648M0RElANMiIiKuKSkJEyYMAEAMHToUHTv3l25AaVDXV0dU6ZMkf5PRFTYyIQQQtlBFISIiAgYGhoiPDw8TUdUoqJs/fr1GD9+PEqXLg0/Pz8YGRkpOyQiojxTUN/fPKdOVIS9f/8ec+bMAQCMHj0aUVFRkMvlSo6KiKjo4SUzoiJs9OjRCA8PR4UKFaCmpoZdu3bBwMAA9vb2hapDtVwuR2BgIIAvwwWwfxMRFTb8q0RURB04cABHjx4F8GVE6pQkIyIiAq6urlLtsMIgNjYWFhYWsLCwQGxsrLLDISJKgwkRURGUmJiIqVOnAgAaNGiQ7uCLbm5uvHxGRJRNTIiIiqDly5cjKCgImpqaaN++fbptIiIipMtURESUOSZEREXMp0+f8McffwAA2rRpAz09vQzbZlQAloiIFDEhIipiZs2ahfDwcJQrVw5NmzbNtC0LuxIRZQ8TIqIi5Pbt29i6dSsAoFevXlBVVc2wLQu7EhFlHxMioiJCLpdj7NixEEKgX79+GDt2bKbtWdiViCj7OA4RURHh4uICLy8v6OnpYfny5ahYsSKcnJzg5uaGiIgIqV1hHIdITU0NY8aMkf5PRFTYsHQHURHw+fNnWFpa4uPHj1ixYgV+/fVXaV7KoIeRkZHQ19fnwIdEVKwU1Pc3f6oRFQFz5szBx48fUadOHamQawoVFRWYm5srJzAiomKCCRFRIeft7Y3NmzcD+FLItShWixdC4OPHjwAAIyMjyGQyJUdERKSICRFRIfL15a/KlStjzJgxEEKgb9++sLW1VXaI3yQmJgblypUDAERFRUFXV1fJERERKWJCRFRI+Pr6pukg/fjxY3h5eUFfXx/Lly9XYnRERMUbEyKiQsDX1xeurq4K02JiYnDq1CkAwJgxY2BiYqKM0IiISgTeikKkZHK5HG5ubmmmX7lyBbGxsShXrhyMjIxYqJWIKB8xISJSssDAQIXLZAAQHByM27dvAwA6duyI6OhoFmolIspHTIiIlOzrAqxyuRxnzpwBANSrV0+6pZ6FWomI8g8TIiIl+7oA6507dxASEgINDQ106NAhw3ZERJR32KmaSMlMTU1hYGCAiIgIxMTE4NKlSwCANm3aSElQUS/UqqamhkGDBkn/JyIqbPiXiUjJVFRUYG9vD1dXV1y+fFnqSG1jYyO1KeqFWjU1NbFz505lh0FElKGi+xeWqBixsrKClZUVfHx8AACdOnWCqqoqDAwM4OTkVKgKtRIRFUc8Q0RUCMjlcqxYsQJCCHTv3h2TJk0qVoVahRCIiYkBAOjo6LB0BxEVOkyIiAqB7du349atWzAwMMCmTZtQoUIFZYeUp2JiYqCnpweApTuIqHAq+j89iYq4T58+YcaMGQCABQsWFLtkiIioKGBCRKREcrkcY8eOxefPn2FpaYnRo0crOyQiohKJCRGRkvj6+sLZ2RkHDx4EAHz//ffYsGEDfH19lRwZEVHJw4SISAl8fX2xf/9+KRlq1KgRzMzMEBERAVdXVyZFREQFjAkRUQFLKeZ648YNvH//Htra2mjfvr1CGzc3NxZzJSIqQEyIiApYYGAgAgMD4eHhAQDo0KEDdHR0FNpERESwmCsRUQHibfdEBSwiIgJnz55FYmIizMzM0LBhw3TbFadirqqqqnB0dJT+T0RU2DAhIipgN27cgJ+fH1RUVNC5c+cMByksTsVctbS0cOjQIWWHQUSUIV4yIypAUVFRWLhwIQDghx9+gLGxcbrtinoxVyKiooYJEVEBmjdvHoKCglC5cmW0bNkyw3ZFvZgrEVFRw7+4RAXk3r17WL16NQDgr7/+Qv/+/WFgYKDQprgWc42OjoZMJoNMJkN0dLSywyEiSoN9iIgKgFwux6hRo5CcnAxHR0d07NgRAGBpaYnAwEBERkYWq2KuRERFDRMiogKwdetW3LhxA3p6etJZIgBQUVGBubm50uIiIqIvmBARfQO5XJ7tMztv376VircuWrQIlSpVKshQiYgoG5gQEeWQr68v3NzcEBERIU0zMDCAvb19un1/JkyYgLCwMDRp0gRjx44tyFCJiCib2FmBKAd8fX3h6uqqkAwByLAG2YkTJ3Do0CGoqqpi27ZtUFPjbxAiosKICRFRNqXUIMtM6hpkERER0hmhKVOmZDgiNRERKV+hSIiio6MxduxYzJkzB87OzoiPj0/TJjw8HL169YKBgQEaNWqEGzduKCFSKskCAwPTnBn6WuoaZDNnzkRwcDCqVauGuXPnFkSIhZaqqio6duyIjh07snQHERVKhSIhGj16NNq1a4dFixahcePGmDlzZpo2y5YtQ8+ePeHu7o4qVaqgW7duHM+EClR2a4tFRkbi+vXr2LhxIwBgy5Yt0NbWzs/QCj0tLS2cOXMGZ86cgZaWlrLDISJKQ+kJUUhICA4dOgQHBwcAgIODAzZv3pzmy6ddu3bo06cPmjRpgr179yI0NBSPHz9WRshUQmW3tpiGhgaGDx8OABgyZAjs7OzyMywiIsoDSk+IPDw8YGRkJP1qNDY2hoaGBry8vBTapf5SMTAwgIGBASpXrpzheuPj4xEREaHwIMoNU1PTNCNLf83AwAD79u2Dr68vypUrh+XLlxdQdERElBtKT4iCg4NRpkwZhWn6+voICQnJcBk/Pz/Y2tqiYsWKGbZZsmQJDA0NpUeVKlXyLGYqmVRUVGBvb59pm2rVqmHJkiUAgHXr1qV5b5dU0dHR0NXVha6uLi91E1GhpPSESCaTpelTkJCQAHV19QyX2bBhQ5a/vGfOnInw8HDp8fr16zyJl0o2KysrODk5pVuDzNHREcuWLUNiYiK6dOmCXr16KSnKwikmJgYxMTHKDoOIKF1KHxTFxMQE4eHhCtOioqJgYmKSbvvz58/Dzs4OFhYWma5XU1MTmpqaeRYnUQorK6t0a5CtW7cO//77L/T09LBhwwbIZDJlh0pERNmk9DNEtra2CAoKQkJCAgBIl8psbGzStH306BECAwPRrVu3Ao2R6GspNcjq1asHc3NzvHjxQro7cvny5bxES0RUxCg9ITIxMYG9vT08PT0BABcuXMCYMWOgpaWFFStW4OnTpwCAFy9eYOvWrWjfvj0CAgJw7949bNq0SZmhEwH4MmDjkCFDEBsbCzs7O4wcOVLZIRERUQ4p/ZIZAGzevBkzZszAzZs38fnzZyxduhQAsH//fpibm8PAwEA6k7RmzRppuQMHDigrZCrmclK8dd26dbh27Rr09PSwfft2XiojIiqCZEIIoewgCkJERAQMDQ0RHh6e5a3TVLLlpHjr8+fPUb9+fcTGxmLTpk0YNWpUQYdbJERHR0NPTw/Alz6Curq6So6IiIqKgvr+VvolM6LCJCfFW3mpLPtUVFTQunVrtG7dOsMzbUREylQoLpkRFQbZLd5qaWkJFRUVrF+/HlevXuWlsmzQ1taGh4eHssMgIsoQEyKi/5OT4q1JSUmYMWMGgC93lZmbmxdAhERElF947pro/2S3eGt4eLjCpbIRI0bkc2RERJTfmBAR/Z/sFm89fPiwwqUy9onJWnR0NIyNjWFsbMzSHURUKPGSGdH/SSnemtlls+joaKxduxYAsGLFCl4qy4GPHz8qOwQiogzxpy3R/8mqeGtSUhJOnjyJ+Ph4dOrUiZfKiIiKESZERKlkVrz1zZs3ePr0KYyMjLBt2zbeVUZEVIzwkhnRV9Ir3hoQEIBff/0VALB161ZUqFBByVESEVFeYkJElI6U4q3Al7vKbG1tIYTAkCFD0L17d6XGRkREeY+XzIiyMGHCBLx69QoWFhZYvXq1ssMhIqJ8wDNEVOJlVsjV1dUVu3fvhoqKCvbs2ZPtW/NJkYqKCqytraX/ExEVNkyIqETLrJCrlpYWhg8fDgCYOXMmfvjhB2WFWeRpa2vj1q1byg6DiChDrHZPJVZKIdf0JCcn4/jx43jw4AG+//57eHh4QF1dvYAjJCIiVrsnykdZFXK9cuUKHjx4gFKlSmHfvn1MhoiIijkmRFQiZVbI9fnz57h+/ToAYMmSJTAzMyvI0IqlmJgYmJubw9zcHDExMcoOh4goDfYhohIpo0KukZGROHbsGACgadOm7DeUR4QQePXqlfR/IqLChmeIqERK724xuVyOY8eOITo6GuXLl0eHDh14VxkRUQnBhIhKpJRCrqldvXoVL1++hLq6OhwdHVG2bFmYmpoqKUIiIipITIioRPq6kOuzZ8/g7u4OAOjUqROMjY1hb2/PMXOIiEoI/rWnEiulkGtSUhKOHj0KAGjSpAlatWoFJycnWFlZKTlCIiIqKOxUTSWahYUFzp07h9jYWNStWxd///03atSowTNDREQlDBMiKtEmTJgAb29vlC1bFmfOnGGfoXwik8lQu3Zt6f9ERIUNEyIqVjKrS/a1HTt2YOvWrZDJZNi/fz+ToXyko6ODR48eKTsMIqIMMSGiYiOzumRf9we6ffs2xowZAwBYuHAh2rdvX6CxEhFR4cKOElQspNQl+3r06YiICLi6usLX11eaFhwcjG7duiE+Ph5dunTBzJkzCzpcIiIqZJgQUZGXVV0yAHBzc4NcLkdMTAy6deuGkJAQ1KlTB3v37mUH6gIQExODOnXqoE6dOizdQUSFEi+ZUZGXWV2yFBEREQgICMCMGTPg7e0NIyMjnDp1Kl8rJ9P/J4TA48ePpf8TERU2/GlMRV5Gdcm+tnz5chw6dAjq6uo4evQoLCws8jkyIiIqKniGiIq87NQbe/jwIQ4fPgwA+Ouvv9CyZcv8DouIiIoQniGiIi+9umSpvX79GsePHwcA/Prrrxg8eHABRUZEREUFEyIq8r6uS5bax48fsW/fPiQlJaFz585YtmxZAUdHRERFARMiKhZS6pKlPlMUGRmJffv2ITY2Fk2bNsWBAwegqqqqxCiJiKiwYh8iKjasrKxgaWmJwMBAvHnzBsOGDcPnz59RvXp1nDlzBrq6usoOscSSyWQwMzOT/k9EVNgwIaJiRUVFBRUrVsTw4cPx+PFjlCtXDm5ubjA2NlZ2aCWajo4OAgIClB0GEVGGeMmMipXExET07t0bly5dgq6uLs6cOYNq1aopOywiIirkeIaIipyMCrgmJydj0KBBOHHiBDQ1NXHy5ElYW1srO1wiIioCmBBRkZJRAdcOHTpg1apV2L9/vzTwop2dnRIjpdRiY2PRqlUrAMA///wDbW1tJUdERKSICREVGSkFXL8WFhaGwYMH49atW1BRUcG+ffvQsWNHJURIGZHL5bh9+7b0fyKiwoYJERUJGRVwlcvlOHPmDLy9vSGTybB9+3Y4OjoqIUIiIirKmBBRkZBeAVe5XI6TJ0/i7t27kMlk6NatG2xtbZUTIBERFWlMiKhI+LqAa3JyMo4fP44HDx5AJpOhZ8+eqFevXrYLvRIREaXGhIiKhNQFXBMSEnDkyBE8ffoUKioqcHR0RO3atdO0IyIiyi4mRFQkpBRwffv2Lfbt24egoCCoqamhV69esLS0BPDlbjNTU1MlR0pEREUREyIqElRUVFCvXj38/vvv+PjxI7S0tNCnTx+pHAQA2NvbQ0WFY40WVkZGRsoOgYgoQ0yIqEi4ffs2Bg0ahI8fP6JUqVLo27cvypUrB+DLmSF7e3tYWVkpOUrKiK6uLj58+KDsMIiIMsSEiAq9/fv3Y8iQIYiLi0Pt2rVx9uxZCCHSjFRNRET0rZgQUaEll8sxZ84cLFmyBADQsWNH7Nu3D4aGhkqOjIiIihv+rKZC6dOnT+jatauUDE2bNg0nT55kMlRExcbGwtbWFra2toiNjVV2OEREaTAhokJFLpfj8OHDqFu3Ls6cOQNNTU3s2bMHy5Ytg6qqqrLDo28kl8vh6ekJT0/PTEt3/Pvvv5DJZLC3ty/A6BRt3LgRFhYW0NLSQpMmTXD16tVcL/PPP/+gS5cuMDExgUwmw/Hjx9OsIykpCXPmzIGFhQW0tbVRtWpVLFiwQOH1WrJkCZo2bQp9fX2UK1cO3bt3x9OnT3O9z0TEhIgKkUePHqFr165wcnLC27dvUaZMGYwdOxZNmjRRdmhUQHbs2IE+ffrA3d0dgYGBBb79gwcPwtnZGbNnz8adO3fQsmVLODg4ZBpLdpaJjo5GgwYNsH79+gzXs2zZMmzevBnr16+Hr68v/vjjDyxfvhzr1q2T2nh6emLs2LG4ceMGLl68iKSkJHTo0AHR0dF58wIQlWAyIYRQdhAFISIiAoaGhggPD4eBgYGyw6GvXLx4ESNHjoS/vz8AoF69eujcuTM0NTUBAE5OTryLrAiLjo6Gnp4eACAqKgq6urrptqlQoQIuX76MuXPnonnz5vjf//5XoHF+9913aNy4MTZt2iRNs7KyQvfu3aXLt7ldRiaT4dixY+jevbvC9M6dO6N8+fLYvn27NO2nn36Cjo4O9uzZk+62P3z4gHLlysHT0xOtWrXKya4SFRkF9f3NM0SkVHK5HBs2bEDnzp3h7+8PNTU1dO3aFT179pSSIQBwc3NjlfRi7uDBg6hQoQJsbGzQr18/uLi44Ft/r/3+++/Q09PL9PH1Za2EhAR4e3ujQ4cOCtM7dOiAf//9N93tfMsyGWnRogUuX74MPz8/AMC9e/dw7do1dOzYMcNlwsPDAQBlypRJM8/V1RXW1tY4fvw4+vTpg1KlSmH16tX477//0KhRI1SoUAG3b9/OUYxExRnvMiOl8fPzw+jRo3HlyhUAX0aj7t69e7p/3CMiIhAYGAhzc/MCjpIKyvbt29GvXz8AQPfu3TFy5EhcvnwZ7dq1AwD8+eefWLVqFYyMjCCEQPPmzbFs2bJ0O9qPGjUKTk5OmW6vUqVKCs8/fvyI5ORklC9fXmF6+fLl8fbt23TX8S3LZGT69OkIDw9HrVq1oKqqiuTkZCxevBh9+vRJt70QApMnT0aLFi1Qt27dNPM7deqE4cOHw8fHBzt27MDx48cxcuRI7Nq1Cz4+Ppg0aRJWrlyJffv25ShOouKKCREVuKioKCxatAgrV65EYmIitLS0YGtrCxsbm0zHE2Lh1uLr6dOn+Pfff+Hi4gIA0NPTQ7du3bBjxw4pIXr48CFWrVqFXr16ITk5GdOmTYOzs7O0TGplypRJN7HODplMpvBcCJFmWl4s87WDBw9i79692LdvH+rUqYO7d+/C2dkZJiYmGDRoUJr248aNw/3793Ht2rV016erq4tSpUrBzs4O2trasLGxQWRkJHr06AEAqF+/Pvbu3ZujGImKM14yowIjl8uxb98+1KpVC8uWLUNiYiIcHBxw7tw5NGvWLMvBFVm4tWjT0dGBjo5OuvO2b9+Opk2bombNmtK0fv364ejRowgNDQXwJSGqU6cOAEBVVRWzZ8/G6dOn013ft1wyMzIygqqqapozO+/fv09zBig3y2Rk6tSpmDFjBn7++WfUq1cPAwYMwKRJk9LthzR+/HicPHkS7u7uqFy5cobrTJ2UfX2XpoqKCi9DE6XChIjynRACJ0+eRKNGjdCvXz8EBwejatWqOHnyJM6cOYNWrVpl2VGOhVuLNl1dXURHRyM6OjpNh+qkpCTs3r0bffv2VZj+448/Ql9fH3///TeEEHjx4gVq1KghzdfW1pb60Hxt1KhRuHv3bqYPa2trhWU0NDTQpEkTXLx4UWH6xYsX8f3336e7nW9ZJiMxMTFpfhSoqqoqJC1CCIwbNw5Hjx7FlStXYGFhkaNtEFHGeMmM8o0QAm5ubpg3bx68vLwAfElspk2bhl9//RVaWloAII074+rqmuG6WLi1+Dp9+jTevXuHunXr4uHDhwrzWrZsie3bt8PBwQGVKlWCurq6NO/ly5eoWrVquuv81ktmkydPxoABA2BtbY3mzZtjy5YtCAwMxKhRo6Q269evx7Fjx3D58uVsLxMVFYXnz59Lz/39/XH37l2UKVNGSvS7dOmCxYsXw9TUFHXq1MGdO3ewcuVKDBkyRFpu7Nix2LdvH06cOAF9fX3pzJShoSG0tbXT7I8QItOO6SXkJmOi7BElRHh4uAAgwsPDlR1KsRcbGyu2bt0qateuLQAIAEJHR0fMmDFDfPr0KcPlHj9+LFauXCnmzZsnPVauXCkeP35cgNFTQevcubP0PsnoMW/ePNG7d2+F5SZPnixmzpyZ5/Fs2LBBmJmZCQ0NDdG4cWPh6empMH/u3LnCzMwsR8u4u7unu1+DBg2S2kRERIiJEycKU1NToaWlJapWrSpmz54t4uPjpTYZvT4uLi5p9uP06dNCXV1djBs3TgQHB4v58+cLAGLLli3i9evXonv37sLIyEhcv349168ZUX4qqO9vjkNEeebp06dwcXGBi4sL3r9/D+BLv58RI0Zg6tSp2epTIZfLERgYyMKtxUxcXBx++uknAMCRI0eks4PZtXjxYsjlcvz2228QQmD//v1YuHAh/vvvP5QqVSofIiaiwqKgvr95yYxy5ePHjzh27Bh27typMO5KlSpVMHHiRAwbNixH9cdUVFR4a30xlJycjLNnz0r/z6lHjx7h2rVrOH78OIQQsLGxgbu7O5MhIsozTIgox169eoVTp07h6NGjCrWpVFRU0LFjRwwePBhdunRR6O9BlBscK4eI8hsTIsrSmzdvcPXqVVy+fBmXL1/GixcvFOY3atQIvXv3xsCBA1GxYsVsr5eXx4iIqLAoFAlRdHQ0pk2bhtKlSyMqKgrLli1TKNuQ4uTJk7h8+TLi4+Ph6OgoDdhGeSMpKQkBAQHw9fXFo0ePcOvWLXh5eSEoKEihnaqqKpo1a4aePXuiR48e33Trr6+vL9zc3BARESFNMzAwgL29PWuWERFRgSsUnaoHDhyIHj16oEePHti9ezfu3r2LlStXKrR58uQJBg4ciJs3b0IIAWtra5w6dSrN8PsZKemdqoUQiI6OxufPnxEcHIygoCDpERgYiCdPnsDPzw8JCQlpllVRUUHdunVhZ2eHtm3bZmvcoMz4+vpmeos9C7kWP9kp7kpElJ6C+v5WekIUEhKCatWqITQ0FFpaWvjw4QPMzMzw7t07hZGJR40ahXLlymHBggUAvozHUbp0aSxatChb20l5QXft2gUdHR1p/I3C9G9O2iYkJCAuLi7dR2xsLGJiYhAWFobPnz8jNDQUoaGhSExMzPJ10tLSgqWlJaysrNCkSRPY2NigcePG0pdZbsnlcqxZs0bhzNDXDAwMMHHiRF4+K0aYEBHRtyoxd5l5eHjAyMhIug3X2NgYGhoa8PLyQtu2baV2V65cwfTp06XnNWrUwJEjRzJcb3x8POLj46XnKSPaplcTqCRRU1NDxYoVYWJigsqVK8PExASVKlVCtWrVULNmTVSpUiXNEP9yuTzTBCYnXr16Jd2Sn5G4uDg8evQIZmZmebJNUr7o6Gjp/xEREd90pxkRlUwp3z/5ff5G6QlRcHBwmhFl9fX1ERISkmm79NqktmTJEsyfPz9vgy0GkpKS8Pr1a7x+/Ro3b95UdjgZWrp0qbJDoHxiYmKi7BCIqAj69OlTjoZxySmlJ0QymSzNIG0JCQlpbtn+ul16bVKbOXMmJk+eLD0PCwuDmZkZAgMD8/UFLWwiIiJQpUoVvH79ukT1neJ+c79LAu4397skCA8Ph6mp6TeV48kJpSdEJiYmaQo0RkVFpfkV+XW7yMjITH9pampqpnunmqGhYYl6I6UwMDDgfpcg3O+ShftdspTU/c7vfqVK77Vqa2uLoKAg6e6mlMtgNjY2Cu3atm0LPz8/6fnz58/Rpk2bgguUiIiIii2lJ0QmJiawt7eHp6cnAODChQsYM2YMtLS0sGLFCjx9+hTAl7vMLl26BOBLPxgvLy8MHz5caXETERFR8aH0S2YAsHnzZsyYMQM3b97E58+fpQ61+/fvh7m5OSwtLdGgQQMMHjwYU6ZMQUJCAlatWoUKFSpkexuampqYO3duupfRijPuN/e7JOB+c79LAu53/u630schIiIiIlI2pV8yIyIiIlI2JkRERERU4jEhIiIiohKPCRERERGVeIXiLrOcio6OxrRp01C6dGlERUVh2bJl6fY+P3nyJC5fvoz4+Hg4OjqiXbt2uZ6nTPm939HR0ahSpQpCQ0MBAEePHkWPHj3yf8eykBf7/fHjR/zxxx+Ijo7Ghg0bsr2cMuX3fhfX4x0eHo5hw4bh/PnzqFatGjZt2oRmzZpluZyy5fd+F9fjnZycjLFjx2Lfvn0wMzPD7t270ahRoyyXU7b83u/ierxTO3DgADZv3gwPD48cLZchUQQNGDBAHD16VAghxK5du8SkSZPStPH19RVNmzYVcrlcJCcni0aNGomgoKBczVO2/NxvIYRYu3atcHV1FRcvXhQXL14USUlJBbNjWcjtfgshhLe3t+jatasYNGhQjpZTpvzcbyGK7/GeOXOm2Ldvn7h9+7bo0qWLKFeunIiKispyOWXLz/0Wovge771794r//vtPvH//XnTs2FH88MMP2VpO2fJzv4Uovsc7RXBwsKhbt65o3bp1jpbLTJFLiIKDg4WWlpaIjY0VQgjx/v17oa2tLSIiIhTajRw5Uvz222/S8zFjxojZs2fnap4y5fd+JyUliQ4dOognT57k967kSF7sd4rffvstTWJQnI93ivT2uzgf78uXL0vTw8PDhbq6uvDy8spyOWXK7/0uzsc7JiZGmn7q1CnRqlWrbC2nTPm938X5eKcYP3682LBhg0JClNvjXeT6EHl4eMDIyEgq9GpsbAwNDQ14eXkptLty5QrMzMyk5zVq1JBGw/7WecqU3/t94cIFXL9+HbVq1cKPP/6I9+/f5/cuZUte7HeK9OrgFOfjnSK9/S7Ox9vOzk6anlLzqXLlylkup0z5vd/F+Xhra2tL0x8/fozff/89W8spU37vd3E+3gCwdetW9O/fHzo6OjlaLitFLiEKDg5OU/FWX19fqoGWUbvUbb51njLl9347ODggMjISnp6eCA4ORpcuXSCXy/Nrd7ItL/Y7J+svTsc7MyXlePv5+cHW1hYVK1bM0XIFLb/3u7gf70+fPmHBggXYsGGDQhHw4n68M9rv4ny8nz9/joiIiDT1TrNaLjuKXEIkk8mk7DJFQkIC1NXVM22Xus23zlOm/N7vlPmtWrWCu7s7nj9/jhs3buTHruRIXux3TtZfnI53drZR3I/3hg0bsHz58hwvV9Dye79Tli2ux1tfXx8dO3ZE06ZN0aNHD7x+/TpbyylLfu93yrLF7XgnJydj48aNcHZ2ztb6c3q8i1xCZGJiopAJA0BUVBRMTEwybRcZGSm1+dZ5ypTf+52asbExnJycFD5cypIX+52T9Ren451dxfV4nz9/HnZ2drCwsMjRcsqQ3/udWnE83hoaGrC2tsbBgwdRpUoV6cu/uB/vjPY7teJ0vP/9919s3rwZZcuWRalSpTBmzBhcu3YNpUqVynS57CpyCZGtrS2CgoKQkJAAANLpsK9Pn7Vt2xZ+fn7S8+fPn6NNmza5mqdM+b3fX1NTU0PDhg3zche+SV7sd2aK8/HOieJ2vB89eoTAwEB069YtR8spS37v99eK2/FOoaqqigYNGqBSpUo5Wq6g5fd+f624HO+mTZvi8ePHuHv3Lu7evYsFCxbA2toad+/ezXS5bMtJD/HC4qeffhIXLlwQQgjh4uIipkyZIoQQYvny5VKv+rt374qWLVsKIYRITEwUDRo0EG/evMnVPGXLz/0+fvy4ePz4sRBCiCdPnohZs2YV3I5lIbf7neJ///ufGDhwoMK04ny8U6S338X5eD9//lxMnDhR+Pv7C39/f3H37l2xcePGLJdTtvzc7+J8vB89eiQ+ffokhBAiNDRUdO3aVcjl8iyXU7b83O/ifLxTc3FxUbjLLLfHu0hWu//48SNmzJgBc3NzfP78GUuXLoWGhgaaNGmCmTNnwtHREQDg4uKCR48eISEhAT169FDIFL91njLl537Pnj0b69atQ6tWrdCmTRtMmjQp3buTlCEv9tvb2xuTJk3Chw8f8Ndff6FVq1bSvOJ8vDPa7+J6vN+8eQMbGxsEBQUprPPAgQPo3bt3hssVBvm538X1eAPAkCFDcPLkSXTr1g0VKlTAxIkTUa5cOWn9xfF4A5nvd3E+3qnt3LkTO3fuVBiYMTfHu0gmRERERER5qXCkjERERERKxISIiIiISjwmRERERFTiMSEiIiKiEo8JEREREZV4TIiIiIioxGNCRERERCUeEyIiIiIq8ZgQERERUYnHhIiIiEgJWCiicGFCVML9+++/6Nu3L2QyGUxNTdGpUyfUq1cPXbt2lSoI58a1a9fQt29f2NvbZzuePn36QCaToVy5cmjbti0cHR3RtGlT2Nvb4++//06zzK1bt1C2bNk0dZy+RVhYGFasWAFzc3MEBATken1F2YYNGzBx4kRUrVoVI0aMQGJiYoFuPyIiAmvXrkWlSpWgpqaGq1evKsyXy+XYt28fqlatinbt2uHixYsK7x8LCwvs2LEDnz59Snf9O3fuhKamJlq1agV7e3tUqlQJWlpasLe3R/v27WFsbAxbW1tcu3YNP/30E2QyGSwtLbF7925pHbt370atWrUgk8nQo0cPXLt2LVv7dvz4cezZs+fbX5xMFMX38Nq1ayGTyVC3bl1lh5Jjr1+/xp9//omkpKQcL3vhwgWcOXMmH6Kib5LTSrVU/Dx69EgAEFu3bhVCCBEVFSWaNWsm9PX1perD3+rjx4+iZcuWChWJs/Lw4UOFeIQQIjk5WaxcuVKoqKiIXr16SVWdhRAiMDBQDBo0SEREROQqViGEePv2rVi6dKkAIPz9/XO9vsLo/v37WbZxcXERgwYNEkIIcevWLVG+fHnx4sWLfI4sfc+fPxcAROXKlcXHjx/TzJ86dao4e/as9Dzl/bN69epM1+vi4iLOnz8vPR80aJCoVKmS9DwsLEz88ssvQgghfH19BQCxZs2aNOtZvXq1AJCjz4qDg4NUlTuvFcX3cIMGDYRMJhMAxI0bN5QdTrb5+/uLYcOGieTk5G9ex19//SV27NiRh1HRt+IZIoKOjo7Cc11dXcyePRuRkZFYtWpVrtZdtmxZmJub52gZXV3dNNNUVFQwadIkLF68GIcOHcKff/4pzatSpQp27twJfX39XMUKAOXLl0fTpk1zvZ7CKjw8HIsWLcqy3c6dO1GlShUAgLW1Nd6+fYuqVavmd3jpqlatGoyMjBASEoLBgwenma+jowNtbW3pecr7J733UWoWFhbo0KFDhvMNDQ0xaNAgAICWlpa0ra+lbEdTUzOLPfni1atXuHnzJq5evYqnT59ma5mcKGrv4du3b+PevXuYOnUqNDQ0sH37dmWHlC3Jycno168fZs+enatK8iNGjMC+ffvy5Iw85Q4TIkpXjRo1ACBPLkPl5o/F1yZPngxjY2OsXLkSycnJebbe1PIy3sIkISEB/fr1w7t377JsGxQUBFVV1QKIKnvq1KmD2bNn49SpU1izZk2erLN169ZZtrG1tc2TbaW2detW7N27F5qamti2bVuerx8oWu/h7du3Q1VVFZMmTULnzp1x4MABxMTEKDusLO3evRuGhoY5/sGXnlGjRmHq1Km5D4pypeh8aqhAvXz5EgBQq1YtAEBsbCzmzZuHMWPGoGHDhhg8eDAiIiIAAJcuXULfvn0xc+ZMNG7cGOvWrctwvZGRkbCxscEPP/yA7du357hfioaGBuzs7PDmzRv4+PggNDQUy5cvT9Nf4vjx45g4cSImTJgAAwMDrF+/HlevXsXPP/+MESNGwNXVFZUqVUKVKlWy/EWa1f4lJSVhzpw5GDduHDp27AgnJyeEhoZm+rp5e3tj8ODBcHBwwIMHD9CyZUvo6+tj9uzZSEhIwPjx46Gvr4/69evj1atX0rYyOw6p1+nt7Y1GjRrB0NBQ+tL9+++/8ezZM/j5+WHUqFE4ffp0mn319vbGsGHD8O7dO5w8eRLDhg2Dp6cnhgwZAgcHBxw8eBClS5fG0qVLAQDx8fGYNWsWpk6dih49eqBjx4549uxZmnhyso8ZmTdvHuzs7DBt2jT4+Phk2b4wSkpKwqNHj+Dg4ABHR0fs2rULCQkJCm2yOo5Azj5zS5YsgUwmg7m5OW7evAngS/+sXr16oXXr1nj//j2A9D8zQObvucyWy47Y2Fjs378fDg4OqFChAgYPHozIyEgcOnQo2+tQlvXr16Nt27bS89y8321tbXHp0iU8evSooHeDUlP2NTtSPn9/f4U+O+Hh4eK7774TpUuXFkFBQUIIIcaNGyf9//Pnz8LY2FgMHTpUxMTECB0dHakvxpYtW4SKiopCf55BgwZJfYj++ecfMWjQIJGQkJDteL42bdo0AUAcPnxYvHv3TixevFihv0RCQoIwNTWV2h8/flysW7dOPHjwQJiZmQlra2uxZs0acfv2bdGyZUshk8nE3bt3pfbu7u7S+rKzf0OGDBHz588XQnzpf6WjoyPGjRuX6evm5+cnrK2tRY0aNcS2bdvE27dvxYwZMwQAMX36dOHl5SX8/f1F5cqVxejRo6VtZbQ+Ib70pWrRooWwsLAQ69evF2/evBG//PKLMDAwEImJiWmORWbMzMzE3LlzhRBCPH78WFhbW4uqVauK7du3i1mzZglXV1chhBB9+vQRa9eulZYbO3asqFixoggPD//mfUxPSsxv374VFStWFNWrV5eOwdy5c4W7u7vUNqv3T0a+7kOUWmbr3Lp1a7b76xw5ckTs3r1bCCHEtWvXBADptUyR1XHMznsy9XtYCCEcHR1FvXr1FLYzevRo8fDhQyFExp8ZITJ/z2W2XHbs3r1bABBHjhwRQgiRlJQkKlSokG7/qpRj0KtXL4XpDRo0EMePHxdCCPHs2TOhqqoqYmNjpfnbtm0Turq6QldXV6irqws1NTXpeZs2bYS/v7+QyWQiKipKYb1NmzYVe/fuTTfukJAQAUCh71pu3+9lypQRf/zxRzZfOcoPagWfglFhdebMGbx69QrPnj3D999/jyNHjqBSpUoIDAzEiRMnULZsWalts2bNEBMTA3V1dfTt2xdNmjQBABgbG0MulyM0NDRNnx5PT0/s3LkT27Ztg5rat7/1ZDIZgC+3rJYrVw42NjYK8yMjIxEYGIg///wTkyZNQpcuXXD9+nXUrVsXZmZmqFChAiZMmADgy+WL2rVrY/v27Vi7dm2abWW1f/7+/nBxcUFgYCCAL/1J9uzZg8qVK2f6utWoUQOWlpYICgrC0KFDAQADBw7E0qVLYW9vL/UBadGihXTGJbP1AV/6UllYWEAmk2Hs2LEAAEdHR+zcuRPv37+HiYnJN73eVlZWsLKywoMHDzBkyBBp+sOHD7F//36F/lyzZ8/Gpk2bsG7dOsyePTvH+5iV8uXLY//+/Wjbti1GjRqV7l2H+cnFxSXNnWTZjR0Ajh49Kp2R/OGHH9CgQQNs27YNvXr1ktpkdRzLlSuX7c9cihkzZsDa2ho3btxAs2bNkJiYiFevXqFOnToAMv7MZPWey2i57Nq+fTuMjIzQuXNnAICqqioGDBiA5cuX49mzZ9KlewC4d+8ezM3NcenSJSQkJEBDQwOJiYnw9fVFgwYNAAD379+HpaWl1OcLAIYOHSq9/0aMGAE9PT2sXLlSmn/ixAmYmpqm6XP25MkTab1fu3//PgCgQoUK0rRv+UynVrlyZdy5cyebrxzlByZEJOnUqROGDRuWZvrDhw+ho6ODefPmpbvc1q1bcfXqVZw+fRrR0dEAvtwSndqnT5/QqVMnbNiwIVfJEAAp+Ujp5Pv1+sqUKQNnZ2dMmTIF27dvx/z586UvHJlMptAB19LSEpUrV8aLFy/S3Zaamlqm+3fr1i0IIRS+MHr27AkAOHv2bKav29dxp44rhaampnRZMavjAHzpO5K6/0jKH/mvL8vklIqKCgwNDRWmXblyBQAUplesWBGVKlXCrVu3AOR8H7OjdevWWLBgAWbPno127dple7m8MHjw4DSfkW3btuHff//NctmXL1/iwYMHGDlypDRNS0sLFy9eREBAgEJflMyOY1bvyfQ0adIEP/zwA1atWoWDBw/i7Nmz6NixozQ/o89MVu/hzD5rWXn+/Dn++ecfTJgwARoaGtL0wYMHY/ny5dixYweWLFkiTb937x5at26Nly9fwsPDAx06dICvry+0tLSk1+7evXsZJjHAl0Rm1KhRCtPu3buH2rVrK0wLDAxEfHy81GXgaylDOXydROXm/a6rqytdviTlYB8iylJ8fDwCAgLw+fNnhekfP34EAIwfPx6nT5/GkiVL4OjomO46ypYti+HDh2PcuHE5+kX9taSkJLi7u6NSpUqZ/uFbtWoVLly4ADU1NTg5OWHOnDkZti1Xrlymdwhltn8pv0T9/PwUpkdERGT5uuVUXq8vr7x580bheYUKFaCurp6v25w5cyYcHBwwbtw4+Pr65mjZDx8+5FNUmdu2bRsOHTqEnTt3So/Tp09DTU0NO3bsyNG6svOZ+5qzszOOHDmCwMBAHDhwAH369FGYn95nJjvvuZx81lLbsWMHhBD45ZdfFKZbWVnhu+++w65duxRunEhJdrp06YJTp04pTEvdpn79+uluTwiBhw8fppmfXkL0+PFj1KlTJ8MfbymJTm5/aKQmk8kK1Y0MJRETIspS7dq1ER8fj8WLFytM37ZtG65fv47169dj7ty5Wd7ZsmzZMtSoUQO9e/dGfHz8N8WyceNGvH37Fr/99luGfzw+f/6Mf/75B+3bt4ePjw8GDBigcDns67vT3rx5gx9++CHddWW1fymXHDZt2iRNE0Jgz549mb5u3yIv1pdyuTEvpFyqvHz5ssL0T58+wc7OLs+2kx6ZTIY9e/agbNmycHV1zfZyCQkJOH78eP4Flsl2nzx5gpo1aypMNzIygoODA3bs2JHtuyZz8plLrUePHqhUqRLmzZsHVVVVlClTRpqX0Wcmq/dcVp+1jCQnJ2PXrl1o1KgRGjZsmGb+4MGD8ebNG5w7d06alp2E6P79+xn+UHrx4gXi4+Olz2zq9W7atAlGRkbSo1evXpn+4KpUqRKAL8NY5JXIyEhpqAtSDiZEhNjYWADI8FZXS0tL9OjRAytXrkS/fv2wefNm9OzZE40aNZLuptq1axfu378PFxcXAMCDBw9w7949AF/++CUnJ0NDQwMHDhzA06dP4ezsnGE8GcWxc+dOTJkyBRMmTFC47JDyRZLyb0JCAhYvXgwhBNTU1NCzZ0+FL6JHjx5JQ+Z7eHhALpdL1/y/Xl9W+1etWjU4Ojpiy5YtGD9+PHbv3o2ffvoJDRo0yPR1S/26pEi55JH60ocQQmqT1foAIDExMd1LJynr0NXVxatXrxAWFoaLFy9mdAgQGxurcBzkcnmaJLZZs2ZwcHDA6tWrpcs2Pj4+UFVVlX7153QfM4ol5T2aWtmyZXHw4ME0Z6NS4k59J1RKLL/++ivq1auXo+2kXmd6lzpSXpe4uLgM92H37t1o3rx5uvO6dOmC4OBghcQus+OY3c9c6n+BL/1zxo0bBxcXlzRnhzL6zGT1nstoufj4eDRr1gyOjo7plqc4d+4cQkJCYGVlhePHj6d5pJx5TelvFRUVhZcvX6JBgwaoVasWNDQ0cO/ePYWEKCIiAgEBARkmMvfu3UPNmjUVzganrPfkyZO4e/eu9LCzs8s0IWrUqBH09fXT3C2Wm/d7QEAAvvvuuwy3SQVAKV25qdC4fv266N27twAgGjRoIA4cOKAwCnSK0NBQ0b9/f6GrqyvMzc2Fi4uLEOLLXSadO3cW+vr64ueffxaPHz8WZcqUEf379xcJCQnixIkTonLlysLAwEDs379fhIeHi4YNGwoAYsyYMcLPzy/DeCwsLESfPn3EsGHDRLt27YSjo6O4dOmSQvvAwECp/YQJE8THjx/FmzdvBADRtGlTMWvWLDFw4EBpFOHWrVuLJk2aCGdnZzFjxgzRrVs36U4bIYR4/fq1tL7x48eLgICATPdPiC8jGvfp00fo6OgIS0tLcfTo0Sxft0uXLolKlSpJr0tISIgYO3asACD69Okj/Pz8xIkTJ0SVKlWEvr6+dCdSRusTQoiLFy8KExMToa+vL/bv3y9evnwp+vXrJ+3L+/fvhY+PjzAxMRE2NjbizZs3aY7zu3fvxLx58wQAUa1aNXHgwAFx/PhxUalSJaGuri7+/PNPkZSUJLUPCwsTQ4YMEdbW1mL06NFi+PDh0h1J37qPqd2/f19MmDBBaGhoiA0bNoi3b9+mabNy5UrpLrPr16+LXr16CQBCT09PtG7dWrRt21a0aNFCVKhQQRgZGaUZVTgsLExs375dlClTRgAQCxYsUHhfXr9+Xfz8888CgGjcuLH4+++/pXkHDhwQjRs3FgBEz549xbVr19LE9/fff4syZcqIH374QVy9elVh3tOnT0WXLl0EAFG+fHmxffv2bB3HzN6TX7+Hg4ODpe2FhoYKCwsLhWMohMj0M5PZey6j5WJjY4WZmZmoUKGC8PHxSfOadO/eXQDI8qGmpibevn0rrl+/rnAH4OTJk8WCBQuEsbGx8PLyEkIIcfXqVVG2bNk020rxv//9T/z8888K065fvy60tbWluzBT1KxZU1y5ciXDdQkhRL9+/cTUqVOl57l5v798+VJoamqK9+/fZ7pNyl9MiKhEad26tVSSgqikefnypZg2bVqBbW/58uUiICAg1+vZuHGj6Nixo/Tc3d1dmJmZCVVVVRETEyOEEGLDhg2idevWIjY2VnrEx8dLy3Tv3l38/vvvadbbvHlzhWlhYWFCJpOJT58+ZRrTgwcPRP369XO7a1LsY8eOzZN10bfjJTMiohJi27ZtUjmS/BYUFITExESYmZnlel1fd5Zu0aIFIiIiUKNGDamD87179+Dp6QltbW3p8dNPP0nL3L9/P90O1dbW1grTvL29UalSJYU+VumpW7cuunfvnutBJJOSknDs2LE0/bSo4MmESOcCL1Ex1aJFC5iZmRX4+DVEyuLl5YWVK1eiZs2aCAgIwO7du/N9m1FRUbh8+TK6du2apx35Cxu5XI45c+ZgyJAhqF69eo6XF0Jg0aJFcHR0hJWVVT5ESDnBM0RUIiQnJ+Ovv/7CvXv3cOXKFZw8eVLZIREViA8fPuDs2bO4f/9+piU+8pKenh66detWrJMh4Mt4Ub///juePXuW4zJEAHD16lWMHDmSyVAhwTNEREREVOLxDBERERGVeEyIiIiIqMRjQkREREQlHhMiIiIiKvGYEBEREVGJx4SIiIiISjwmRERERFTiMSEiIiKiEo8JEREREZV4TIiIiIioxGNCRERERCXe/wP/v6w37gmXlwAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "x = np.linspace(0,0.05,500)\n",
    "plt.plot(d,pc,'o',label='Data from NLTH Analyses',color='grey')\n",
    "plt.plot(x,lognorm.cdf(x,params[0],scale=params[2]),label='Fitted lognormal distribution',color='k')\n",
    "\n",
    "plt.plot((0.0182,0.0182),(0,0.51),color='k',ls='--')\n",
    "plt.plot((0,0.0182),(0.51,0.51),color='k',ls='--')\n",
    "\n",
    "plt.text(0.0185,0.02,r'$\\Delta_{D}=0.0182$ m')\n",
    "plt.text(0.001,0.53,r'$POE = 0.49$')\n",
    "\n",
    "plt.xlim(0,0.04)\n",
    "plt.ylim(0,1)\n",
    "plt.xlabel(r'Peak Displacement from NLTH Analyses, $\\Delta_{NLTH}$ (m)',fontsize=12)\n",
    "plt.ylabel(r'CDF',fontsize=12)\n",
    "plt.legend()\n",
    "\n",
    "plt.savefig('C:/Users/rober/Documents/ROSE/PostDoc/Tanks/Paper/Figure15.tiff',dpi=300)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "14cf8ff8-26dc-488d-96cf-1e7d9cd1bc68",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "params = lognorm.fit(d,floc=0.0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "c38d847d-bc41-47a4-8fa6-fbddb7a21e1f",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.42791742607650884, 0.0, 0.018010762842673302)"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "params"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "e422847d-8b28-49d3-b844-c3aee7a51620",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.018010762842673302"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "lognorm.ppf(0.5,params[0],scale=params[2])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 62,
   "id": "59272be3-f545-42e5-b86e-e0709b1bacf9",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x149bfc79310>]"
      ]
     },
     "execution_count": 62,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAGdCAYAAACyzRGfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABLC0lEQVR4nO3deVxVdf7H8dcFgYsK1wUNKEfNMRuzQlLSxspMJ6zUaRMXzH3BJW0q06YyS1MnRzNLUdtcwWXMNc3MLUcFTakxS9OfpQKu2AVE1nt+fzgwoqBclnu58H4+Hvcx3ev3nvvhDHXffleTYRgGIiIiIg7i5uwCREREpHJR+BARERGHUvgQERERh1L4EBEREYdS+BARERGHUvgQERERh1L4EBEREYdS+BARERGHquLsAq5ls9lISEjAx8cHk8nk7HJERESkCAzDICUlhcDAQNzcbty3Ue7CR0JCAvXq1XN2GSIiIlIMJ0+e5Lbbbrthm3IXPnx8fIArxfv6+jq5GhERESmK5ORk6tWrl/c9fiPlLnzkDrX4+voqfIiIiLiYokyZ0IRTERERcSiFDxEREXEohQ8RERFxKIUPERERcSiFDxEREXEohQ8RERFxKIUPERERcSiFDxEREXEohQ8RERFxqGKFjw0bNhASEsKvv/5aaJs1a9YwcuRIhgwZwubNm4tbn4iIiFQwdm+vfvbsWbKzs9m7d2+hbX7++WcmTJhATEwMhmHQokUL1q5dy6233lqiYkVERMT12d3zUbduXZ544okbtnn//fcJDQ3FZDLh5uZG69atmT17drGLFBERkYqjWMMubm43ftuWLVuoX79+3vPGjRuzffv2AttmZGSQnJyc7yEiIiKlLykpiaeffppvvvnGqXWUyYTT+Ph4atWqlffcx8eHhISEAttOmjQJi8WS96hXr15ZlCQiIlKp7d69m+bNm/PFF1/Qv39/srKynFZLmYQPk8mE2WzOe56ZmYmHh0eBbceOHYvVas17nDx5sixKEhERqZRsNhvvvfceDz30ECdOnOCPf/wjX3zxRaHfy45g94TToggMDMRqteY9T0lJITAwsMC2Xl5eeHl5lUUZIiIildr58+fp3bs3X375JQDdunVjzpw5+Pr6OrWuMgkfjz76KEeOHMl7fvToUR555JGy+CgREREpwLfffkv37t2Jj4/HbDbzwQcfMGDAAEwmk7NLK96wi2EY+f4XYOrUqRw+fBgg394e2dnZxMbGMnDgwJLWKiIiIjdhs9mYOHEibdu2JT4+niZNmhATE8PAgQPLRfCAYvR8pKamsnDhQgDmz5/P8OHD8fPzIyoqigYNGtCkSRPuvfde+vbty8svv0xmZibTp0/H39+/1IsXERGR/zlz5gy9evXi66+/BqBXr17MmjWL6tWrO7my/EzG1d0X5UBycjIWiwWr1er0MSkRERFXsXXrVnr06MHp06fx9vbmo48+ok+fPg7r7bDn+1tnu4iIiLiwnJwcxo8fT/v27Tl9+jRNmzZl37599O3bt9wMs1yrTCacioiISNlLTEykZ8+ebN26FYB+/foxc+ZMqlat6uTKbkzhQ0RExAV9/fXXhIeHc/bsWapVq0ZkZCTh4eE3fE+OzSD2eBJnU9Kp62MmpGEt3N0c3zui8CEiIuJCsrOzeeutt3j33XcxDIN77rmHpUuXcuedd97wfRsPJjJ+7SESrel5rwVYzIzr1JTQZgFlXXY+mvMhIiLiIk6dOkW7du2YOHEihmEwePBg9uzZU6TgEbFof77gAXDamk7Eov1sPJhYlmVfR+FDRETEBWzYsIGgoCC+/fZbfHx8iI6OJjIyEm9v7xu+L8dmMH7tIQpa2pr72vi1h8ixOW7xq8KHiIhIOZaVlcWrr77K448/zoULFwgODmb//v2EhYUV6f2xx5Ou6/G4mgEkWtOJPZ5UShXfnOZ8iIiIlFMnTpygW7du7N69G4Dhw4czdepUu85EO5tSePAoTrvSoPAhIiJSDq1Zs4Y+ffpw8eJFLBYLn3zyCc8884zd16nrY755IzvalQYNu4iIiJQjmZmZ/O1vf6NLly5cvHiRli1bcuDAgWIFD4CQhrUIsJgpbEGtiSurXkIa1ip2zfZS+BARESknjh8/Tps2bZg+fToAL774Ijt37qRhw4bFvqa7m4lxnZoCXBdAcp+P69TUoft9KHyIiIiUAytXrqR58+bs3buXmjVrsnr1aqZNm4anp2eJrx3aLIDZ4cH4W/IPrfhbzMwOD3b4Ph+a8yEiIuJEGRkZvPzyy3z44YcAtG7dmqioKOrXr1+qnxPaLIAOTf21w6mIiEhldvToUcLCwti/fz8Ao0ePZsKECXh4eBTYvqTbo7u7mWjdqHap1F4SCh8iIiJOsHTpUgYOHEhKSgq1a9dmwYIFPP7444W2L0/bo5eU5nyIiIg40OXLlxkyZAjdunUjJSWFBx98kLi4uJsGj/K0PXpJKXyIiIg4yOHDh2nVqhVz5szBZDLRe+jf+HrzN9x2222Fvqc8bo9eUgofIiIiDrBo0SKCmgfzww8/4Fa1BnWee5ttPu1o+88dN+y5KI/bo5eUwoeIiEgZSktLo3///vTq1Yv0y2l4/eEeAvp+gHfD5sDNh07K4/boJaUJpyIiImXk0KFDPPfccxw6dAhMJiwPdMfyQBgmN/e8NgZXNvsav/YQHZr6X7d6pTxuj15S6vkQEREpZYZh8Nlnn9GiRQsOHTpE7Tp1uSVsIjXa9MgXPPLaU/jQSXncHr2kFD5ERERKUWpqKr1796Zfv35cvnyZDh068M+orzDXv+em7y1o6KQ8bo9eUgofIiIidsixGew+doHVcfHsPnYh3yqTH374gZYtW7Jw4ULc3NyYOHEiGzdu5I4G9Yp07cKGTsrb9uglpTkfIiIiRVTYRl9vPvknTu5ex8iRI0lPT+fWW28lKiqKBx98EPjf0Mlpa3qBS2ZNXAkSNxo6KU/bo5eUwoeIiEgR5G70dW14SDibxHNh3bn00w4AOnbsyIIFC/Dz88trkzt0ErFoPybIdw17hk7Ky/boJaVhFxERkZsobKOvzDPHSJg/8krwcHNj8pQprFu3Ll/wyFXRhk5KQj0fIiIiN3HtRl+GYZB6YD1JWz6GnGzcfetQp/NoHnqmL25uhf+9viINnZSEwoeIiMhNXL0KxZaeyoWNM0k7/G8AvP94P7UfH4W7t0+RNvqqKEMnJaHwISIilUpxjqXPXYWSkXiE86unkG09A25VqNm2Lz4tOmMymfK1kxtT+BARkUqjuMfSt2xQEw6u5/SGeWDLxt1yC3W6vIpXwB1A0VaryP9owqmIiFQKxT2WPikpiWefeZrf1s8GWzZV73iAwD4z8gUPcL2NvpxJPR8iIuJy7B06udmx9IWdrbJnzx7CwsI4ceIEnp6eDHh5HN9VC+F0ckZeG/8i9JxIfgofIiLiUoozdGLPsfStG9XGZrPxz3/+k9dee43s7GwaNWrEsmXLCA4OLtacEclP4UNERFxGYRt95Q6dFLZfhj3H0p8/f54+ffqwfv16AMLCwpg7dy6+vr6AVquUBs35EBERl3CzoRO4MnRy9VkruYq6CiXxcBxBQUGsX78eLy8v5syZQ1RUVF7wkNKh8CEiIi7BnqGTa93sWHoMG7YDKxneswvx8fE0adKE2NhYBg0alLeMVkqPwoeIiLgEe4ZOrnWjY+ltl37nzLJxnNz0KTk5OYSHh7Nv3z7uueeekpYshVD4EBERl1DUoRN7jqVPP/EDp+ePJP3XA3h7e/Ppp5+yYMECqlevXio1S8E04VRERFxCaR5Lv/voOT6aPoVlS9/HZrPRtGlTli1bxl133VVm9cv/qOdDRERcwo2GTuzZ6Ovc2TOMi+hOdOQ0bDYbffv2JTY2VsHDgRQ+RETEZZT0WPrNmzdz7733smXLFqpVq8aCBQv49NNPqVatWlmWLdfQsIuIiDhcSTbqKs6x9NnZ2bz11lu8++67GIbB3XffzbJly7jzzjtL60cSOyh8iIiIQxX3cLer2bPRV3x8PD169GDHjh0ADB48mOnTp+Pt7W1/8VIqNOwiIiIOU9zD3Yr9eRs3EhQUxI4dO/Dx8SEqKorIyEgFDydT+BAREYcoyQ6l9srKymLMmDF07NiR8+fP07x5c7777ju6detW4mtLySl8iIiIQ5Rkh1J7nDhxgrZt2zJlyhQAhg0bxq5du2jcuHGJriulR3M+RETEIUqyQ2lRrV27lt69e3Px4kUsFguffPIJzzzzTLGvJ2VDPR8iIuIQJd2h9EYyMzN56aWX6Ny5MxcvXqRly5bs379fwaOcUvgQERGHuNnhbiaurHq50Q6lBTl+/DgPPvgg06ZNA2DUqFHs3LmT22+/vWQFS5lR+BAREYcorR1Kr7Zy5UqaN29ObGwsNWrUYNWqVUyfPh1PT8/SKVrKhMKHiIg4TEl3KM2VkZHBiBEjeOaZZ7BarbRq1Yq4uDi6dOlSFmVLKdOEUxERcaji7FB6taNHjxIWFsb+/fsBeOWVV5g4cSIeHh5lWbaUIoUPERFxOHt2KL3asmXLGDBgACkpKdSuXZv58+fzxBNPlEGFUpY07CIiIuXe5cuXiYiIICwsjJSUFNq0aUNcXJyCh4tS+BARkXLt8OHDtGrVisjISEwmE6+99hpbt27ltttuc3ZpUkwadhERkXJr8eLFDB48mEuXLlGnTh0WLVrEX/7yF2eXJSWkng8RESl30tLSGDBgAOHh4Vy6dIm2bdvy/fffK3hUEHb3fFy6dInRo0dTs2ZNUlNTmTJlCl5eXte1e/fdd7HZbLi5uZGWlsY777yDyVT0tdsiIlK+5diMYq9YuZFDhw7RtWtXfvzxR0wmE2+++SZvvPEG7u7upVC1lAd2h4+IiAieeuopnnrqKRYsWMDYsWPzdpXLtW7dOo4fP868efMA6NevH8uWLSMsLKx0qhYREafaeDCR8WsP5TsoLsBiZlynpkXeq6Mgn3/+OcOGDSMtLQ1/f38WL15Mu3btSqNkKUfsGnZJSEhg+fLldOzYEYCOHTsSGRlJSkpKvnY//fRTvte8vb2xWq2lUK6IiDjbxoOJRCzaf90Jtaet6UQs2s/Gg4l2XzM1NZXevXvTt29f0tLSaN++PXFxcQoeFZRd4WPbtm34+flhNl/Zma5OnTp4enoSGxubr13nzp1Zs2YNK1as4OLFi5w/f55evXoVeM2MjAySk5PzPUREpHzKsRmMX3sIo4A/y31t/NpD5NgKalGw//znP7Rs2ZIFCxbg5ubGhAkT+Oqrr7jllltKpWYpf+wKH/Hx8dSqlf/AHx8fHxISEvK91qRJE5YsWcLzzz9Pjx49+Oyzz/D29i7wmpMmTcJiseQ96tWrZ+ePICIixZFjM9h97AKr4+LZfexCkQJD7PGk63o8rmYAidZ0Yo8n3fRahmEwb948QkJC+PnnnwkMDGTr1q38/e9/x81N6yEqMrv+3zWZTHm9HrkyMzML3NI2PT2dNWvWcOzYMbp37052dnaB1xw7dixWqzXvcfLkSXtKEhGRYth4MJE2U7bQfd4eRkbH0X3eHtpM2XLTIZOzKYUHD3vapaSk0LNnTwYNGkR6ejodO3YkLi6Ohx56qMg/g7guu8JHYGDgdXM3UlNTCQwMzPfatm3b2LdvH+3bt2fbtm0cOHDgukmpuby8vPD19c33EBGRslOSORt1fcyF/llR2x04cIDg4GCioqJwd3dnypQprFu3jjp16hTtBxCXZ1f4aNu2LadOnSIzMxMgb7glJCQkX7sVK1bQuHFj4EpgmTBhAjt27CiNekVEpARKOmcjpGEtAixmCltQa+LKqpeQhrWu+zPDMJg1axatW7fm6NGj1KtXjx07djB69GgNs1Qydvd8hIaGsn37dgA2bdrE0KFDMZvNTJ06lcOHDwMQFBTEgQMH8t5nMpmuCygiIuJ4JZ2z4e5mYlynpgDXBZDc5+M6Nb1uvw+r1UrXrl0ZNmwYGRkZdOrUibi4OB544IFi/iTiyuze5yMyMpIxY8YQExNDUlISkydPBiAqKooGDRrQpEkT+vXrx4kTJ5g8eTJ+fn6cPHmS1157rdSLFxER+5TGnI3QZgHMDg++bp8P/0L2+di3bx9du3bl+PHjeHh4MGXKFEaNGqWNJysxk2EYRV8P5QDJyclYLBasVqvmf4iIlLLdxy7Qfd6em7aLGtjqpkfe32yHU8Mw+OCDD3jllVfIysqiQYMGLF26VD3hFZQ93986WE5EpBLJnbNx2ppe4LwPE1d6MAqas3EtdzdToQHl4sWL9OvXj1WrVgHw9NNP88knn1CjRo1i1y4Vh2b4iIhUIsWds2GPPXv20Lx5c1atWoWnpyczZ85kxYoVCh6SR+FDRMSFFWejsNw5G/6W/Mth/S1mZocHF/tsFpvNxtSpU3nwwQf57bffaNSoEbt27WL48OGa3yH5aNhFRMRFleRwt9BmAXRo6l9qp9JeuHCB3r17s379egDCwsKYO3eu5u5JgTThVETEBeVuFHbtf8Bzo0NJejDstXPnTrp3786pU6fw8vJixowZDBo0SL0dlYw9398adhERcTFlcbhbcdhsNiZNmpS3AeUdd9xBTEwMgwcPVvCQG1L4EBFxMaV5uFtxnT17lscff5zXXnuNnJwcwsPD+e6777j33nvL7DOl4tCcDxERF1Nah7sV1/bt2+nevTuJiYl4e3vz4Ycf0rdvX/V2SJGp50NExMWUxuFuxZGTk8Pbb79Nu3btSExM5E9/+hN79+6lX79+Ch5iF/V8iIi4mNLcKKyoTp8+Tc+ePdmyZQsAffv2ZebMmVSrVq3UPkMqD/V8iIg4UXH26XDERmFX27x5M0FBQWzZsoWqVauyYMECPv30UwUPKTb1fIiIOElJ9+mw53C34sjOzmb8+PFMnDgRwzC4++67WbZsGXfeeWeJry2Vm/b5EBFxgtLap+Nmh7sVV3x8PD169GDHjh0ADBw4kBkzZuDt7V3ia0vFpIPlRETKsZvt02Hiyj4dHZr63zRI3Ohwt+LauHEjvXr14vz581SvXp25c+fSvXv3Uv0Mqdw050NExMHKwz4dBcnKymLs2LF07NiR8+fPExQUxP79+xU8pNSp50NExMGcvU9HQU6ePEm3bt3YtWsXAMOGDWPq1KmYzaW7XFcEFD5ERBzOWft0FGbdunX07t2bpKQkfH19+eSTT3j22Wcd8tlSOWnYRUTEwXL36ShsNoeJK6teSnOfjoJkZmby0ksv0alTJ5KSkmjRogUHDhxQ8JAyp/AhIuJgjt6noyC//vorDz74INOmTQNg1KhR7Ny5k9tvv73MPlMkl8KHiIgT5O7T4W/JP7TibzEXeZltcX3xxRc0b96c2NhYatSowapVq5g+fTpeXl5l9pkiV9OcDxERJwltFkCHpv5lsk9HQTIyMnjllVeYOXMmAK1atSI6Opr69euXyeeJFEbhQ0TEicpin46CHDt2jLCwML777jsAXn75Zd599108PDzK/LNFrqXwISJSwS1fvpwBAwaQnJxM7dq1mT9/Pk888YSzy5JKTHM+REQqqPT0dIYOHUrXrl1JTk6mTZs2xMXFKXiI0yl8iIhUQEeOHKFVq1bMnj0bgLFjx7J161Zuu+02J1cmomEXEZEKZ8mSJQwePJjU1FTq1KnDwoULeeyxx5xdlkge9XyIiFQQaWlpDBw4kJ49e5Kamkrbtm2Ji4tT8JByR+FDRKQC+Omnn7j//vv5+OOPMZlMvPnmm2zevJnAwEBnlyZyHQ27iIi4uPnz5zN06FDS0tK45ZZbWLJkCe3atXN2WSKFUs+HiIiLunTpEr1796ZPnz6kpaXRvn17vv/+ewUPKfcUPkREXNB//vMfWrRowYIFC3Bzc+Odd95h48aN3HLLLc4uTeSmNOwiIuJCDMPgk08+YcSIEaSnpxMYGMiSJUt4+OGHnV2aSJEpfIiIuIiUlBSGDBnCkiVLAAgNDWXBggXUqVPHyZWJ2EfDLiIiLiAuLo777ruPJUuW4O7uzuTJk1m/fr2Ch7gk9XyIiJRjhmEQGRnJiy++SEZGBvXq1SM6OpoHHnjA2aWJFJvCh4hIOWW1Whk4cCDLly8HoFOnTnz22WfUrl32p+CKlCUNu4iIlEP79u0jODiY5cuXU6VKFf75z3+yevVqBQ+pENTzISJSjhiGwcyZM3n55ZfJysqifv36LF26lPvvv9/ZpYmUGoUPEZFy4uLFi/Tr149Vq1YB8NRTT/HJJ59Qs2ZN5xYmUso07CIiUg7ExMTQvHlzVq1ahaenJx988AH/+te/FDykQlL4EBFxIsMw+Oc//0mbNm347bffuP3229m1axcjRozAZDI5uzyRMqFhFxERJ7lw4QJ9+vRh3bp1AHTt2pW5c+disVicXJlI2VLPh4iIE/z73/8mKCiIdevW4eXlxezZs4mOjlbwkEpB4UNExIFsNhuTJ0/m4Ycf5tSpU9xxxx3ExMQwZMgQDbNIpaFhFxERBzl37hzPP/88GzduBKBnz57Mnj0bHx8fJ1cm4lgKHyIiDrB9+3Z69OhBQkIC3t7ezJw5k379+qm3QyolDbuIiJShnJwc3nnnHdq1a0dCQgJ/+tOfiI2NpX///goeUmmp50NEpIycPn2a8PBwvvnmGwD69OnDhx9+SLVq1ZxcmYhzKXyIiJSBb775hp49e3LmzBmqVq3K7Nmzef75551dlki5oGEXEZFSlJOTw5tvvkmHDh04c+YMzZo1Y9++fQoeIldRz4eISClJSEigR48ebN++HYCBAwcyY8YMvL29nVyZSPmi8CEiUgq++uorwsPDOX/+PNWrV2fOnDn06NHD2WWJlEsadhERKYHs7GxeHTOG0NBQzp8/T+M/3U3s3n0KHiI3oPAhIlJMJ0+e5N6QB/jHlCkAVG/+BBmPj6ffF6fYeDDRydWJlF8KHyIixbB+/XruuvteDh3Yi8mzKn5dxlD7LxGYqnhy2ppOxKL9CiAihVD4EBGxQ1ZWFi+//DJPPvkkKdaLePr/kYA+M6h2Z5u8NsZ//3f82kPk2IyCLyRSidk94fTSpUuMHj2amjVrkpqaypQpU/Dy8iqwbVpaGnPmzKF27drccccdtGrVqsQFi4g4y6+//kq3bt2IiYkBwOe+ztRs2xdTFY/r2hpAojWd2ONJtG5U28GVipRvdvd8RERE0L59eyZMmEBwcDBjx44tsF1SUhLPPfccf/3rX3n++ecVPETEpa1atYrmzZsTExNDjRo1GDPtY2q1H1Rg8Lja2ZR0B1Uo4jrsCh8JCQksX76cjh07AtCxY0ciIyNJSUm5rm23bt145ZVXaNiwYelUKiLiBBkZGYwcOZKnnnqK33//nfvvv58DBw7QufNfi/T+uj7msi1QxAXZFT62bduGn58fZvOVf5nq1KmDp6cnsbGx+dqtW7eOX375hdjYWB5//HHGjh1LVlZWgdfMyMggOTk530NEpDw4duwYf/7zn/nggw8AeOmll9ixYwcNGjQgpGEtAixmCjsazgQEWMyENKzlsHpFXIVd4SM+Pp5atfL/i+Tj40NCQkK+1xYvXsz999/PqFGjWLhwIYsXL+aNN94o8JqTJk3CYrHkPerVq2fnjyAiUvqWL19OcHAw3333HbVq1WLt2rVMnToVT09PANzdTIzr1BTgugCS+3xcp6a4u+nkWpFr2RU+TCZTXq9HrszMTDw88o95/vjjj7Rp0wZPT09q167NwIEDWbBgQYHXHDt2LFarNe9x8uRJO38EEZHSk56eztChQ+natSvJycn8+c9/Ji4ujieffPK6tqHNApgdHoy/Jf9/F/0tZmaHBxPaLMBRZYu4FLtWuwQGBmK1WvO9lpqaSmBgYL7XsrOzycnJyXt+zz33kJSUVOA1vby8Cl0tIyLiSL/88gtdu3YlLi4OuPKXo/Hjx1/3F6yrhTYLoENTf2KPJ3E2JZ26PleGWtTjIVI4u8JH27ZtGTRoEJmZmXh6euYNt4SEhORrd8899/DLL7/870OqVOHOO+8shXJFRMpGVFQUgwYNIjU1lTp16rBw4UIee+yxIr3X3c2k5bQidrBr2CUwMJDQ0NC8Exs3bdrE0KFDMZvNTJ06lcOHDwMwcuRI1q9fT0ZGBgC7du1ixIgRpVy6iEjJXb58mYEDB9KjRw9SU1N5+OGHiYuLK3LwEBH72b3JWGRkJGPGjCEmJoakpCQmT54MXPlbQ4MGDWjSpAmtW7fmnXfeYcSIETRp0gSz2Uz//v1LvXgRkZL46aef6Nq1KwcPHsRkMvHGG2/wxhtvUKWKDvwWKUsmwzDK1d6/ycnJWCwWrFYrvr6+zi5HRCqoBQsWEBERQVpaGrfccguLFy/m0UcfdXZZIi7Lnu9vne0iIpXKpUuX6Nu3L7179yYtLY1HH32UuLg4BQ8RB1L4EJFK4+DBg7Rs2ZLPP/8cNzc33n77bb766iv8/f2dXZpIpaKBTRGp8AzD4NNPP2X48OGkp6cTGBjIkiVLePjhh51dmkilpPAhIhVaSkoKERERLF68GIDHHnuMhQsXUqdOHSdXJlJ5adhFRCqs77//nhYtWrB48WLc3d2ZNGkSX375pYKHiJOp50NEKhzDMJgzZw6jRo0iIyOD2267jejoaP785z87uzQRQeFDRCoYq9XKoEGDWLZsGQBPPvkkn3/+ObVrF74DaY7N0PboIg6k8CEiFcZ3331HWFgYx44do0qVKkyZMoUXX3wRk6nwILHxYCLj1x4i0Zqe91qAxcy4Tk11MJxIGdGcDxFxeYZhMHPmTB544AGOHTtG/fr12blzJ3/7299uGjwiFu3PFzwATlvTiVi0n40HE8u6dJFKSeFDRFzaxYsXeeaZZ3jhhRfIzMzkr3/9KwcOHOD++++/4ftybAbj1x6ioC2ec18bv/YQObZytQm0SIWg8CEiLismJobg4GC++OILPDw8mDFjBitXrqRmzZo3fW/s8aTrejyuZgCJ1nRijyeVYsUiAgofIuKCDMNg2rRptGnThl9//ZXbb7+dXbt28cILL9xwmOVqZ1MKDx7FaSciRacJpyLiUi5cuECfPn1Yt24dAM899xzz5s3DYrHYdZ26PuZSbSciRaeeDxFxGbt27aJ58+asW7cOLy8vZs2axdKlS+0OHgAhDWsRYDFTWD+JiSurXkIa1ipRzSJyPYUPESn3bDYbU6ZM4aGHHuLkyZM0btyYPXv2EBERUeRhlmu5u5kY16kpwHUBJPf5uE5Ntd+HSBlQ+BCRcu3cuXM88cQTjBkzhpycHHr06MF3331HUFBQia8d2iyA2eHB+FvyD634W8zMDg/WPh8iZURzPkSk3NqxYwfdu3cnISEBs9nMhx9+SL9+/Yrd21GQ0GYBdGjqrx1ORRxI4UNEyp2cnBwmTZrEuHHjsNls3HnnnSxfvpxmzZqVyee5u5lo3ajw7ddFpHQpfIhIuXLmzBnCw8PZvHkzAL179+ajjz6iWrVqTq5MREqLwoeIlBtbtmyhR48enDlzhqpVqzJr1ix69+5dpPfqcDgR16HwISJOl5OTw9tvv80777yDYRg0a9aMpUuX0rRp0yK9X4fDibgWrXYREadKSEigffv2vP322xiGwYABA4iJibEreOhwOBHXovAhIk7z1VdfERQUxLZt26hevTqLFy9m3rx5VK1atUjv1+FwIq5J4UNEHC47O5vXXnuN0NBQzp07x7333st3331Hjx497LqODocTcU2a8yEiDnXq1Cm6d+/Ozp07AYiIiGDatGmYzfafoaLD4URck8KHiDjM+vXr6d27NxcuXMDX15d58+bRtWvXYl9Ph8OJuCYNu4hImcvKyuKVV17hySef5MKFC9x3333s37+/RMEDdDiciKtS+BCRMvXbb7/x0EMPMXXqVABeeOEF/v3vf9OoUaMSX1uHw4m4JoUPESkzq1evJigoiD179lCjRg1WrlzJjBkz8PLyKrXP0OFwIq5Hcz5EpNRlZmYyevRoZsyYAUBISAhLly6lQYMGZfJ5OhxOxLUofIhIqfq///s/wsLC2LdvHwAvvfQS7777Lp6enmX6uTocTsR1KHyISKlZsWIF/fv3Jzk5mVq1avH555/TqVOnIr1XZ7OIVB4KHyJSYunp6bz00kvMmjULgAceeIDo6Gjq1atXpPfrbBaRykUTTkWkRH755Rdat26dFzzGjBnDtm3b7AoeOptFpHJR+BCRYouKiiI4OJi4uDj8/PzYsGEDkyZNwsPDo0jv19ksIpWTwoeI2O3y5csMGjSIHj16kJqaykMPPURcXByhoaF2XUdns4hUTprzISJ2+fnnn+natSv/+c9/MJlMvP7667z55puY3NzZfeyCXRNGdTaLSOWk8CEiRbZgwQIiIiJIS0vjlltuYdGiRbRv377YE0Z1NotI5aRhFxG5qUuXLtG3b1969+5NWloa7dq1Iy4uLi94FHfCqM5mEamcFD5EJE+OzWD3sQusjotn97EL5NgMfvzxR0JCQvj8889xc3Pj7bffZtOmTfj7+5d4wqjOZhGpnDTsIiLA9XttGIZBlaPbiP/yIzLS0wkICGDJkiW0bds27z32TBgtbPfR3LNZrh228dc+HyIVlsKHSAVS3F1Cc4dOcvsnbBlpJG2axaVD2wAIfuBhNnyxjLp16+Z7X2lNGNXZLCKVi8KHSAVR3Emf1w6dZJ79P86tnkJ2UjyY3KjxUC88OvSktl+d695bmhNGdTaLSOWhOR8iFUBJJn3mDp0YhkFK3AYSF7xEdlI87j5+3NJjEpZWz3E6JbPAvTY0YVREikPhQ8TFlXTS59mUdGwZaZxf8w+SvvoIcrLwbtSSgL4fYL7trnztrqUJoyJSHAofIi6upLuEXvztMImfjyTt52/BzZ2aj/SjzjNv4O7tm69dYUMnuRNG/S35/9zfYmZ2eLAmjIrIdTTnQ8TFFXfSp2EYfPTRR7z00ktkZ2bi7luXOp1H43XrnfnambgSJG40dKIJoyJiD4UPERdXnEmfv//+O/3792flypUAtG4Xyqm7++Burp5v+MaeoRNNGBWRotKwi4iLs3fSZ2xsLM2bN2flypV4eHgwY8YM/r35S+YOeEhDJyLiEOr5EHFxuZM+IxbtxwSF9ly4mWD69Om8+uqrZGVlcfvtt7N06VJatGgBaOhERBxH4UOkArjZLqEhgV506dKFtWvXAvDss8/y8ccfY7FY8l1HQyci4ggKHyIVRGE9FzF7dhP0eDdOnjyJl5cX06dPZ8iQIZhM6tEQEedQ+BCpQK7uubDZbEyd+h6vvfYaOTk5NG7cmGXLlhEUFOTcIkWk0lP4EKmAzp07R+/evdmwYQMA3bt3Z86cOfj4+Di5MhERhQ+RCmfHjh10796dhIQEzGYzM2fOpH///hpmEZFyQ0ttRSoIm83GxIkTeeSRR0hISODOO+8kNjaWAQMGKHiISLlid/i4dOkSw4YN4/XXX2fUqFFkZGTcsH10dDRt27Ytbn0iUgRnzpwhNDSU119/HZvNxvPPP8/evXu5++67nV2aiMh17A4fERERtG/fngkTJhAcHMzYsWMLbZuQkMDEiRNLVKCI3NiWLVsICgri66+/pmrVqnz22WfMnz+f6tWrO7s0EZEC2RU+EhISWL58OR07dgSgY8eOREZGkpKSUmD7yZMnExERUfIqReQ6OTk5vPXWW7Rv357Tp09z1113sXfvXvr06ePs0kREbsiu8LFt2zb8/Pwwm69swVynTh08PT2JjY29ru28efMIDw+natWqN7xmRkYGycnJ+R4icmMJCQm0b9+e8ePHYxgG/fv3JzY2lqZNmzq7NBGRm7IrfMTHx1OrVv6TLX18fEhISMj32tGjR0lOTiYkJOSm15w0aRIWiyXvUa9ePXtKEql0Nm3aRFBQENu2baNatWosWrSIjz/++KZBX0SkvLArfJhMprxej1yZmZl4eHjkPc/JyWHWrFmMGjWqSNccO3YsVqs173Hy5El7ShKpcHJsBruPXWB1XDy7j10gx3bltJbs7Gz+/ve/Exoayrlz57j33nvZv38/PXv2dHLFIiL2sWufj8DAQKxWa77XUlNTCQwMzHu+a9cuIiMj+fTTT4Er4SQzM5MaNWrw+++/X3dNLy8vvLy8ilG6SMWz8WDideezBFjMDG1Zg7njR7Fz504AhgwZwrRp0/D29nZWqSIixWYyDMO4ebMrEhISuOOOO0hKSsLT05OEhAQaN27MhQsX8npE0tPTOX36dN57VqxYwYoVK4iOjqZBgwY3/Yzk5GQsFgtWqxVfX1/7fyIRF7XxYCIRi/Zz7b+Ql4/t5fz66dguJ+Pj48PHH39M165dnVKjiEhh7Pn+trvnIzQ0lO3bt9OhQwc2bdrE0KFDMZvNTJ06lU6dOtGkSZN8ISN3gmpRgodIZZVjMxi/9lC+4GHkZPP7jgUkx64EoFpgY/Zt+5I7Gv/ROUWKiJQSu7dXj4yMZMyYMcTExJCUlMTkyZMBiIqKokGDBjRp0qTUixSp6GKPJ+Ubasm2nuXcmilkJhwGwOe+TtRs248LbjWdVaKISKmxa9jFETTsIpXR6rh4RkbHAZD2yx4ufPk+tvRUTF7V8Os4kqpNHgBgRrcgugTd6sRKRUQKVmbDLiJSNur6mDFysri47XNS9q0GwDPgDvw6j8ajhn++diIirk7hQ6QcqIOVC9FjuHTqv8MsLf9KzYd7Y3K/sozdBPhbzIQ0rHWDq4iIuAaFDxEn+9e//kX//v25ZLXiZq6O3xMv4v3H+/P+PPc82nGdmuLuptNpRcT12X2wnIiUjvT0dIYPH86zzz6L1WrlgQce4LM1W7j9vofztfO3mJkdHkxoswAnVSoiUrrU8yHiBL/88gthYWEcOHAAgFdffZV33nkHDw8Pej5qEHs8ibMp6dT1uTLUoh4PEalIFD5EHCw6OppBgwaRkpKCn58fCxcuJDQ0NO/P3d1MtG5U24kVioiULQ27iDjI5cuXGTx4MN27dyclJYWHHnqIuLi4fMFDRKQyUPgQcYCff/6Z+++/n7lz52IymXj99df55ptvuPVW7dkhIpWPhl1EytjChQuJiIjg0qVL1K1bl8WLF9O+fXtnlyUi4jTq+RApI5cuXaJfv348//zzXLp0iXbt2hEXF6fgISKVnsKHSBn48ccfCQkJ4bPPPsPNzY3x48ezadMmAgK0XFZERMMuIqXIMAw+/fQzhg0fTkb6Zfzq3kJ0VBSPtnvE2aWJiJQbCh8ipSQ1NZUu3fuwZd2/ADA3aI75yZd4PcZGVt1EbRImIvJfGnYRKUCOzWD3sQusjotn97EL5NhufPjzDz/8wJ1333sleJjcqPHQ89TtOh73ajU4bU0nYtF+Nh5MdFD1IiLlm3o+RK6x8WAi49ceItGanvdagMXMuE5Nr+u9MAyDuXPnMnLkSDIyMnCvXhu/LqMx33bX/9pw5XyW8WsP0aGpv3YrFZFKTz0fIlfZeDCRiEX78wUPoMDei+TkZLp3786QIUPIyMjA+/YWBPT9IF/wyGUAidZ0Yo8nlfWPICJS7qnnQ+S/cmwG49ceoqABlmt7L76PO0BYWBhHjx6lSpUqhI8YwxaPEEymG+f5synpN/xzEZHKQD0fIv8Vezzpuh6PqxlAwu+XGT3+H7Ru3ZqjR4/yhz/8gW+//ZZBw0bdNHgA1PUxl2LFIiKuST0fIv91s14JW3oqFzZ8wLQjuwDo0qULn376KbVq1SLHZhBgMXPaml5gz4kJ8LdcOaFWRKSyU8+HyH/dqFciI/EIiZ+PJO3ILqp4ePD+++/zxRdfUKvWlTDh7mZiXKemwJWgcbXc5+M6NdVkUxERFD5E8oQ0rEWAxZwvPBiGQfLe1ZxeNJps6xm8avrz7bc7GTlyJCZT/iAR2iyA2eHB+Fvyhxh/i5nZ4cHa50NE5L807CLyX7m9FxGL9mMCsi+ncOHL97l8NAaAqnc8wPzPPqHV/XcWeo3QZgF0aOpP7PEkzqakU9fnylCLejxERP7HZBjGjXdPcrDk5GQsFgtWqxVfX19nlyOV0MaDibz80Qp+XvIOOcnnwL0KDZ6IYNaEMXS8O9DZ5YmIlEv2fH+r50PkKjabjR++XMjP814jJyeHgHoNePfDT+n1ZFv1XoiIlBKFD5H/On/+PM8//zwbNmwAoFu3bsyZM0c9cCIipUwTTkWAb7/9lqCgIDZs2IDZbGbu3LksWbJEwUNEpAwofEilZrPZmDhxIm3btiU+Pp4mTZoQExPDwIEDr1vNIiIipUPDLlJpnTlzhl69evH1118D0KtXL2bNmkX16tWdXJmISMWm8CGV0tatW+nRowenT5/G29ubWbNm0adPH2eXJSJSKWjYRSqVnJwcxo8fT/v27Tl9+jR33XUX+/btU/AQEXEg9XxIpZGYmEjPnj3ZunUrAP369WPmzJlUrVrVyZWJiFQuCh9SKXz99deEh4dz9uxZqlWrRmRkJOHh4c4uS0SkUlL4kAopx2YQezyJxIuprPnsAxbMno5hGNxzzz0sW7aMJk2aOLtEEZFKS+FDKpyNBxMZv/YQJ0+e4vza98g49SMAjz/XixXz5+Dt7e3kCkVEKjeFD6lQNh5MJGLRftKO7eX8+unYLidj8vTGL3QEh25/iO3Hfie0mcKHiIgzKXxIhZFjMxi36geStn5KcuxKADxvaYRfl1fxqHnlQLjxaw/Roam/zmkREXEihQ+pMNb8+we+nz2KjISfAfC5rxM12/bDVMUDAANItKYTezyJ1o1qO7FSEZHKTeFDKoQ1a9bwfK/nyUi2YvKqhl/HkVRt8kCBbc+mpDu4OhERuZo2GROXlpmZyd/+9je6dOlCarIVz4DGBPSZUWjwAKjrY3ZghSIici31fEi5lbtc9mxKOnV9zIQ0rJVvrsbx48cJCwtj7969AIwcNYrdtR7j7KUcjAKuZwL8LVeuIyIizqPwIeVS7nLZROv/hkgCLGbGdWpKaLMAVq5cSb9+/bBardSsWZPPP/+czp075612MUG+AJIbWcZ1aqrJpiIiTmYyDKOgvyQ6TXJyMhaLBavViq+vr7PLESfIDRDX/mKaACM7i3sS1rA26jMAWrduTXR0NH/4wx/yvf9GwUVEREqfPd/f6vmQciXHZjB+7aECh00yLyZwfvUUfjtzDIDRo0czYcIEPDw88rULbRZAh6b+NxyyERER51H4kHIl9nhSvh6LXJd+2sGFjTMxMi/j5u3LP2bO4aX+3Qq9jrubSctpRUTKKYUPKVeuXQZry8rg4pZ5pMZtBMDrtrvw6/QKf7zvQWeUJyIipUDhQ8qVq5fBZl04xbnVk8k69ytgwtK6K5Y2PTC5uWu5rIiIC1P4kHIlpGEtAixmju7awIWvPsLISsetag38nnwJ74bNtVxWRKQCUPiQciUj/TK+ez/m/LpoALz+cA9+nV6mSvVaWi4rIlJBKHxIuXHo0CGee+45Dh06hMlk4tZ2vXALfgaTmztwpcdDy2VFRFyfwoc4nWEYfP755wwbNozLly/j7+/PkiVLeOjhtlouKyJSASl8iFOlpqYydOhQFi5cCECHDh1YtGgRdevWBdByWRGRCkgHy4nT/PDDD7Rs2ZKFCxfi5ubGxIkT2bhxY17wEBGRikk9H+JwhmEwb948Ro4cSXp6OrfeeitRUVE8+KD27hARqQwUPsShkpOTGTx4MNHRV1azdOzYkQULFuDn5+fkykRExFE07CIOc+DAAe677z6io6OpUqUK//jHP1i3bp2Ch4hIJaOeDylzhmEwa9Ys/va3v5GZmckf/vAHoqOjad26tbNLExERJ7C75+PSpUsMGzaM119/nVGjRpGRkXFdG6vVynPPPYevry/Nmzdnz549pVKsuJ7ff/+drl27Mnz4cDIzM+ncuTMHDhxQ8BARqcTsDh8RERG0b9+eCRMmEBwczNixY69rM2XKFJ5++mm2bt1KvXr16NKlC5cuXSqVgsV17N27l+DgYFasWIGHhwfTp09n1apV1KqlrdFFRCozk2EYRlEbJyQk0KhRIy5evIjZbObcuXPUr1+fM2fO4OPjk9duy5YttGvXDrgywdDPz49///vftGzZ8qafkZycjMViwWq14uvrW4wfSZzNMAxmzJjB6NGjycrKomHDhixdurRI//+LiIhrsuf7266ej23btuHn54fZfOVE0Tp16uDp6UlsbGy+drnBA8DX1xdfX19uu+22Aq+ZkZFBcnJyvoe4rqSkJJ566ilefPFFsrKyeOaZZ9i/f7+Ch4iI5LErfMTHx1/XZe7j40NCQkKh7zly5Aht27YlIKDg8zgmTZqExWLJe9SrV8+ekqQc2bNnD82bN2f16tV4enry4Ycfsnz5cmrUqOHs0kREpByxK3yYTKa8Xo9cmZmZeHh4FPqejz76iPfee6/QPx87dixWqzXvcfLkSXtKknLAZrPx3nvv8eCDD3LixAkaNWrE7t27GTZsGCaTzmIREZH87FpqGxgYiNVqzfdaamoqgYGBBbb/6quvaNeuHQ0bNiz0ml5eXnh5edlThpQj58+fp0+fPqxfvx6AsLAw5s6dq/k6IiJSKLt6Ptq2bcupU6fIzMwEyBtuCQkJua7tjz/+yIkTJ+jSpUsplCnl0bfffktQUBDr16/Hy8uLOXPmEBUVpeAhIiI3ZFf4CAwMJDQ0lO3btwOwadMmhg4ditlsZurUqRw+fBiAY8eOMW/ePDp06MCvv/7K999/z+zZs0u/enEKm83Gu+++yyOPPEJ8fDxNmjQhNjaWQYMG5RtmybEZ7D52gdVx8ew+doEcW5EXVomISAVm9w6nkZGRjBkzhpiYGJKSkpg8eTIAUVFRNGjQAF9f37wekhkzZuS9L/csD3FtZ8+epVevXmzatAmAXr16MWvWLKpXr56v3caDiYxfe4hEa3reawEWM+M6NSW0WcGTj0VEpHKwa58PR9A+H+XXtm3b6NGjB4mJiXh7e/PRRx/Rp0+f6yaVbjyYSMSi/Vz7i5XbanZ4sAKIiEgFU2b7fEjllJOTw/jx43n00UdJTEykadOm7N27l759+14XPHJsBuPXHroueAB5r41fe0hDMCIilZjCh9zQ6dOn+ctf/sJbb72FzWajX79+7N27l7vuuqvA9rHHk/INtVzLABKt6cQeTyqjikVEpLzTqbZSqM2bN9OzZ0/Onj1LtWrVmD17Nr169brhe86mFB48itNOREQqHvV8yHWys7N5/fXX+ctf/sLZs2e5++672bdv302DB0BdH/NN29jTTkREKh6FD8knPj6eRx99lIkTJ2IYBoMHDyYmJoY777yzSO8PaViLAIuZwvY1NXFl1UtIQ51sKyJSWSl8SJ4NGzYQFBTEjh078PHxISoqisjISLy9vYt8DXc3E+M6NQW4LoDkPh/XqSnubtp2XUSkslL4ELKysnj11Vd5/PHHOX/+PM2bN+e7776jW7duxbpeaLMAZocH42/JP7TibzFrma2IiGjCaWV34sQJunfvzq5duwAYPnw477333nUHCNortFkAHZr6E3s8ibMp6dT1uTLUoh4PERFR+KjE1q5dS+/evbl48SIWi4VPPvmEZ555Jl+bHJtR7ADh7maidaPaZVG6iIi4MIWPSigzM5OxY8cybdo0AFq2bEl0dDS33357vnbaIl1ERMqC5nxUMsePH+fBBx/MCx6jRo1i586dBQaPiEX7r9sw7LQ1nYhF+9l4MNFhNYuISMWi8FGJrFy5kubNmxMbG0vNmjVZvXo106dPx9PTM187bZEuIiJlSeGjEsjIyGDEiBE888wzWK1WWrVqxYEDB+jcuXOB7bVFuoiIlCWFjwru8JFfuPe+ED788EMAXn7lFXbs2EH9+vULfY+2SBcRkbKk8FGBvTZ1DnfdE8ThH3/AzduXus+OY2fNx/jm8Pkbvk9bpIuISFlS+KiALl++zBNhvZn0yhByMtLwuq0pAX0+wLtRyyJNGNUW6SIiUpYUPiqYw4cP06pVK75ctgAw4ds6jFu6T6KKrx9QtAmj2iJdRETKksJHBbJ48WLuu+8+fvjhB9yqWqjb9W1qPtQLk5t7vnZFmTCqLdJFRKSsaJOxCiAtLY0XXniBTz75BIC7Wz5AUsgQqlS/8bDIzSaMaot0EREpCwofLu7QoUN07dqVH3/8EZPJxJtvvkn7HkMJ/3TvTd9blAmj2iJdRERKm8KHC/v8888ZNmwYaWlp+Pv7s3jxYtq1a0eOzSDAYua0Nb3AjcJMXBk+0YRRERFxBs35cEGpqan07t2bvn37kpaWRocOHYiLi6Ndu3aAJoyKiEj5pvDhYv7zn//QsmVLFixYgJubGxMmTGDjxo3ccsst+dppwqiIiJRXGnZxEYZh8PHHH/PCCy+Qnp5OYGAgUVFRPPTQQ4W+RxNGRUSkPFL4cAEpKSkMHjyYqKgoADp27Mj8+fOpU6fOTd+rCaMiIlLeaNilnDtw4ADBwcFERUXh7u7OlClTWLduXZGCh4iISHmkno9yyjAMZs+ezYsvvkhmZib16tUjOjqaBx54wNmliYiIlIjCRzlktVoZMGAAK1asAKBz58589tln1KqlpbEiIuL6FD7KmX379tG1a1eOHz+Oh4cHkyZPpnXn5/n2xGXqXrygCaMiIuLyFD7KCcMw+OCDD3jllVfIysqiQYMGjHr3I6J/9WLmxzF57QIsZsZ1aqqlsiIi4rI04bQcuHjxIk8//TSjRo0iKyuLp59+mvcWb+D97w0SrfnPXzltTSdi0X42Hkx0UrUiIiIlo/DhZHv27KF58+asWrUKT09PZs6cydJly/nn9vgCt0bPfW382kPk2ApqISIiUr4pfDiJzWZj6tSpPPjgg/z22280atSI3bt3M3z4cPb+evG6Ho+rGUCiNZ3Y40mOK1hERKSUaM6HE1y4cIHevXuzfv16AMLCwpg7dy6+vr7AzY+6z1XUdiIiIuWJwocD5NiMvC3OEw/HMfHlIZw6dQovLy9mzJjBoEGDMJn+t4KlKEfd29NORESkPFH4KGMbDyYyfu0hEn5PI3nPCn7/dhEYNm5tcDvrV63k3nvvve49IQ1rEWAxc9qaXuC8DxNXDogLaah9P0RExPVozkcZ2ngwkYhF+zmVcJqzy8bx+44FYNiodtcjuD81hUT3ugW+z93NxLhOTYErQeNquc/HdWqq/T5ERMQlKXyUkRybwfi1h7h84gcSP3+B9F8PYKriRe2OL1D7ib/h5ul9wxUroc0CmB0ejL8l/9CKv8XM7PBg7fMhIiIuS8MuZWT30XP89OVnWP8dBYYNj9p/wK/Lq3jWqQ/kX7FS2Kmzoc0C6NDUP2++SF0fs3Y4FRERl6fwUQZOnz7N0F7PYo39NwDV7m5PrfZDcPO8foLozVasuLuZCg0nIiIirkjho5Rt3ryZ8PBwzpw5g8nDTK2/DKV6s3aFtteKFRERqWwUPorg6qWyhQ19ZGdnM378eCZOnIhhGNx99924PfoiVq+6WrEiIiJyFYWPm8hdKnv1jqPXHu4WHx9Pjx492LFjBwCDBg3i/fffZ/ux34lYtB8T5AsgWrEiIiKVmckwjHJ1QEhycjIWiwWr1Zq346ez5C6VvfYG5caF2eHBcOp7evXqxfnz56levTrz5s2jW7du+a5xs/AiIiLi6uz5/q404aMoQyfXtm8zZUvhZ6zkZJO1N5qE7dEABAUFsWzZMho3blzizxYREXE19nx/V4phl+L0PsQeTyo0eGQnn+P8mn+QEf8TAMOGDWPq1KmYzQVPHtWKFRERkf+p8JuM5Q6dXBskTlvTiVi0n40HEwt8X2FLYNOOxpD42QtkxP+EybMqo/8RyYcfflho8BAREZH8KnT4yN1ltKBxpdzXCttl9NolsEZOFklbPubcv97Blp6Cp39jAvp+wF+ffrb0CxcREanAKvSwy42GTuDGu4xefbhblvUM51ZPITPxCAA+LbpQq20fAmr5aKmsiIiInSp0z8fNdg+9Ubvcw90uHdlF4mcvkJl4BDevatR5+nVqPzoQk7uHlsqKiIgUQ4Xu+Sjq7qEFtcvIyODLuZM498VMADwDm1Cn86tUsdTFX0tlRUREiq1Ch4+rh07s2WX02LFjhIWF8d133wHw0ssv06X/SySl52iprIiISAlV6PCRO3Rizy6jy5cvZ8CAASQnJ1O7dm3mz5/PE0884ciyRUREKrQKPecDrhxLPzs8GH9L/qEVf4uZ2eHBeUMn6enpDB06lK5du5KcnEybNm2Ii4tT8BARESllFbrnI1doswA6NPUvdJfRI0eO0LVrV77//nsAxo4dy9tvv02VKpXi9oiIiDhUpfl2LWyX0SVLljB48GBSU1OpU6cOCxcu5LHHHnNChSIiIpWD3eHj0qVLjB49mpo1a5KamsqUKVPw8vK6rt2aNWv45ptvyMjI4Nlnn6V9+/alUnBpSUtL44UXXuCTTz4BoG3btixevJjAwEAnVyYiIlKx2R0+IiIieOqpp3jqqadYsGABY8eOZdq0afna/Pzzz0yYMIGYmBgMw6BFixasXbuWW2+9tdQKL4mffvqJrl27cvDgQUwmE2+++SZvvPEG7u7uzi5NRESkwrNrwmlCQgLLly+nY8eOAHTs2JHIyEhSUlLytXv//fcJDQ3FZDLh5uZG69atmT17dulVXQLz58+nRYsWHDx4EH9/fzZv3sxbb72l4CEiIuIgdoWPbdu24efnl3eIWp06dfD09CQ2NjZfuy1btlC/fv28540bN2b79u0FXjMjI4Pk5OR8j7Jw6dIlevfuTZ8+fUhLS6N9+/bExcXRrl27Mvk8ERERKZhd4SM+Pp5atfJvyOXj40NCQsIN2xXUJtekSZOwWCx5j3r16tlTUpHNmTOHBQsW4ObmxoQJE9i4cSO33HJLmXyWiIiIFM6u8GEyma47Oj4zMxMPD48btiuoTa6xY8ditVrzHidPnrSnpCJ74YUXCAsLY+vWrfz973/XMIuIiIiT2DXhNDAwEKvVmu+11NTU61aIXNsuJSWl0FUkXl5eBa6WKW1VqlQhOjq6zD9HREREbsyuno+2bdty6tQpMjMzAfKGUkJCQvK1e/TRRzly5Eje86NHj/LII4+UtFYRERGpAOwKH4GBgYSGhuZNHt20aRNDhw7FbDYzdepUDh8+DMCQIUPYvHkzANnZ2cTGxjJw4MBSLl1ERERckd37fERGRjJmzBhiYmJISkpi8uTJAERFRdGgQQOaNGnCvffeS9++fXn55ZfJzMxk+vTp+Pv7l3rxIiIi4npMhmEUdNq80yQnJ2OxWLBarfj6+jq7HBERESkCe76/K/yptiIiIlK+KHyIiIiIQyl8iIiIiEMpfIiIiIhDKXyIiIiIQyl8iIiIiEMpfIiIiIhDKXyIiIiIQyl8iIiIiEPZvb16WcvdcDU5OdnJlYiIiEhR5X5vF2Xj9HIXPlJSUgCoV6+ekysRERERe6WkpGCxWG7Yptyd7WKz2UhISMDHxweTyVSq105OTqZevXqcPHlS58aUId1nx9B9dgzdZ8fQfXacsrrXhmGQkpJCYGAgbm43ntVR7no+3NzcuO2228r0M3x9ffXL7QC6z46h++wYus+OofvsOGVxr2/W45FLE05FRETEoRQ+RERExKEqVfjw8vJi3LhxeHl5ObuUCk332TF0nx1D99kxdJ8dpzzc63I34VREREQqtkrV8yEiIiLOp/AhIiIiDqXwISIiIg6l8CEiIiIOVe42GSsNly5dYvTo0dSsWZPU1FSmTJlS4KzeNWvW8M0335CRkcGzzz5L+/btnVCt6yrKfbZarQwYMICvvvqKRo0aMXv2bFq1auWkil1TUX+fc0VHRxMZGcm2bdscV2QFYM99TktLY86cOdSuXZs77rhDv9N2KOp9fvfdd7HZbLi5uZGWlsY777xT6rteVwYbNmxg3LhxLFu2jAYNGhTYxinfhUYF1KtXL2PlypWGYRjG/PnzjRdffPG6Nj/99JPRsmVLw2azGTk5OUbz5s2NU6dOObpUl1aU+zx27FhjyZIlxr59+4xOnToZdevWNVJTUx1dqksryn3OFR8fbzRr1sx4+OGHHVRdxVHU+3zhwgXj8ccfN/7v//7PkeVVGEW5z2vXrjUGDBiQ97xv375GdHS0w2qsKM6cOWOsWbPGAIzjx48X2MZZ34UVLnzEx8cbZrPZuHz5smEYhnH27FnD29vbSE5Oztdu8ODBxhtvvJH3fOjQocbf//53h9bqyop6n7/55pu8f7ZarYaHh4cRGxvr0FpdWVHvc64RI0YYH330kcKHney5zx06dDC2bt3q4AorhqLe53/84x9GWFhY3vOhQ4cac+bMcWitFUVOTs4Nw4ezvgsr3JyPbdu24efnh9lsBqBOnTp4enoSGxubr92WLVuoX79+3vPGjRuzfft2h9bqyop6n9u1a5f3z7nnCJT12T0VSVHvM8C8efMIDw+natWqji7T5RX1Pq9bt45ffvmF2NhYHn/8ccaOHUtWVpYzSnZJRb3PnTt3Zs2aNaxYsYKLFy9y/vx5evXq5YySXd7NDnhz1ndhhZvzER8fT61atfK95uPjQ0JCwg3bFdRGClfU+3y1I0eO0LZtWwICAsq6vAqjqPf56NGjJCcnExISwqFDhxxZYoVQ1Pu8ePFi7r//fkaNGkX//v1p3rw5hmEwefJkR5brsop6n5s0acKSJUvo0aMHDz/8MP/617/w9vZ2ZKmVhrO+Cytcz4fJZMpL1bkyMzPx8PC4YbuC2kjhinqfr/bRRx/x3nvvlXVpFUpR7nNOTg6zZs1i1KhRDq6u4ijq7/OPP/5ImzZt8PT0pHbt2gwcOJAFCxY4slSXZs9/N9LT01mzZg3Hjh2je/fuZGdnO6rMSsVZ34UVLnwEBgZitVrzvZaamkpgYOAN26WkpFzXRgpX1Puc66uvvqJdu3Y0bNjQEeVVGEW5z7t27SIyMpLatWtTo0YNhg4dys6dO6lRo4aDq3VdRf19zs7OJicnJ+/5PffcQ1JSkkNqrAiKep+3bdvGvn37aN++Pdu2bePAgQNMmzbNkaVWGs76Lqxw4aNt27acOnWKzMxMgLzuo5CQkHztHn30UY4cOZL3/OjRozzyyCOOK9TFFfU+w5W/LZ44cYIuXbo4tMaKoCj3uWXLlhw6dIi4uDji4uJ4++23adGiBXFxcc4o2SUV9ff5nnvu4Zdffsl7XqVKFe68807HFeriinqfV6xYQePGjYErX44TJkxgx44dji22knDWd2GFCx+BgYGEhobmTZjZtGkTQ4cOxWw2M3XqVA4fPgzAkCFD2Lx5M3DlbzOxsbEMHDjQaXW7mqLe52PHjjFv3jw6dOjAr7/+yvfff8/s2bOdWbpLKcp9NpvNNGjQIO+RO6GvsDX9cr2i/j6PHDmS9evXk5GRAVzpdRoxYoTT6nY1Rb3PQUFBHDhwIO99JpOpwL/YyM0Z/z071rjqDNny8F1Y4SacAkRGRjJmzBhiYmJISkrKmwwWFRVFgwYNaNKkCffeey99+/bl5ZdfJjMzk+nTp+Pv7+/kyl3Lze6zr69v3t90ZsyYkfe+6OhoZ5Xskory+ywlV5T73Lp1a9555x1GjBhBkyZNMJvN9O/f38mVu5ai3Od+/fpx4sQJJk+ejJ+fHydPnuS1115zcuWuJzU1lYULFwIwf/58hg8fjp+fX7n4LjQZV8chERERkTJW4YZdREREpHxT+BARERGHUvgQERERh1L4EBEREYdS+BARERGHUvgQERERh1L4EBEREYdS+BARERGHUvgQERERh1L4EBEREYdS+BARERGHUvgQERERh/p/90yF7qyMRqsAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "x1 = np.linspace(0,1,1000)\n",
    "plt.plot(pc,lognorm.cdf(d,params[0],scale=params[2]),'o')\n",
    "plt.plot(x1,x1,color='k')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "id": "17bd70c2-a343-4a68-9aa3-568a6f46449e",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "lilfor = np.zeros(len(d))\n",
    "\n",
    "for i in range(len(d)):\n",
    "    lilfor[i] = abs(pc[i]-lognorm.cdf(d[i],params[0],scale=params[2]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 64,
   "id": "5d71d795-bd7f-4d8f-926e-e8b4248fb2e2",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.07744635678031964"
      ]
     },
     "execution_count": 64,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "max(lilfor)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 71,
   "id": "324e4884-47b3-406f-a9b5-b0ce3a66176c",
   "metadata": {},
   "outputs": [],
   "source": [
    "w = np.zeros(30)\n",
    "\n",
    "for i in range(len(w)):\n",
    "    uplift = np.loadtxt('./Results/UpliftI_GM'+str(i+1)+'.out')\n",
    "    w[i] = max(abs(uplift[:,1]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "id": "1c350663-5b31-4ce4-83d0-4fd838252568",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "w = np.sort(w)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 73,
   "id": "43ed4870-0687-470f-b5bc-aa1c5883cc08",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.06138457666666666"
      ]
     },
     "execution_count": 73,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.mean(w)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 74,
   "id": "88ba9f1e-53ef-4ce4-a431-75836602709f",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.057103"
      ]
     },
     "execution_count": 74,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.median(w)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 78,
   "id": "be84969b-05c1-47a8-8f16-4c754282a7ef",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "pcw = np.linspace(0,1,len(w))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 82,
   "id": "ffde18a5-a862-46fb-b1e5-2672fe1a4e2f",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x149bfc15690>]"
      ]
     },
     "execution_count": 82,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAGdCAYAAACyzRGfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAAA+uElEQVR4nO3deXhU5f3+8fdMtmFJAiFAiESCiChFKiBRrEvYNGARF6i1on6tooDaorVKUIsICgjVtr8iWFpXFBS0LWBVFAy4EtmUTTZRIAlrcLJA1jm/Pw4ZGDJJZpLJZJb7dV1z5cyZZ4bnFDq5Pec5n4/FMAwDERERET+xNvUEREREJLwofIiIiIhfKXyIiIiIXyl8iIiIiF8pfIiIiIhfKXyIiIiIXyl8iIiIiF8pfIiIiIhfRTb1BM7kcDjIzc0lNjYWi8XS1NMRERERDxiGQWFhIcnJyVittZ/bCLjwkZubS0pKSlNPQ0REROph3759dOzYsdYxARc+YmNjAXPycXFxTTwbERER8URBQQEpKSnO3+O1CbjwUXWpJS4uTuFDREQkyHiyZEILTkVERMSvFD5ERETErxQ+RERExK8UPkRERMSvFD5ERETErxQ+RERExK8UPkRERMSvFD5ERETErwKuyJiIiEioqHQYZO/J51BhCe1ibaR1TiDCqr5l9Qof77//PpMmTeLtt98mNTXV7ZglS5awYsUKSktLGTFiBIMGDWrIPEVERILKB5vzmLx0K3n2Eue+DvE2Jg3rTkaPDk04s6bndfg4dOgQFRUVfP311zWO+e6775g6dSpr1qzBMAwuvvhili5dyllnndWgyYqIiASDDzbnMXb+eowz9h+wlzB2/nrmjOod1gHE6zUf7dq149prr611zF/+8hcyMjKwWCxYrVb69evHnDlz6j1JERGRYFHpMJi8dGu14AE4901eupVKh7sR4aFeC06t1trftnLlSjp16uR83rVrV1atWuV2bGlpKQUFBS4PERGRYJW9J9/lUsuZDCDPXkL2nnz/TSrANMrdLjk5OSQkJDifx8bGkpub63bstGnTiI+Pdz5SUlIaY0oiIiJ+caiw5uBRn3GhqFHCh8ViwWazOZ+XlZURFRXldmxmZiZ2u9352LdvX2NMSURExC/axdrqHuTFuFDUKLfaJicnY7fbnc8LCwtJTk52OzYmJoaYmJjGmIaIiIjfpXVOoEO8jQP2ErfrPixAUrx52224apQzHwMHDmTHjh3O57t27aJ///6N8UeJiIgElAirhUnDugNm0Dhd1fNJw7qHdb2PeoUPwzBcfgLMmjWL7du3AzBmzBg+/vhjACoqKsjOzmb06NENnauIiEhQyOjRgTmjepMU73ppJSneFva32UI9LrsUFRXx+uuvA/Dqq69y//33k5iYyIIFC0hNTaVbt278/Oc/58477+Thhx+mrKyM559/nqSkJJ9PXkREJFBl9OjA4O5JqnDqhsU4/fRFACgoKCA+Ph673U5cXFxTT0dEREQ84M3vb/V2ERGRRqX+JnImhQ8REWk06m8SgBwOqCiB6OZNNgWFDxERaRTqb+Jn5SVQmAuFB6DoIBQdOvnz9O1D5uP8a+Hm15tsqgofIiLic3X1N7Fg9jcZ3D1Jl2DqYhhw/CgU5EJhnpufeWboOHHM888sOtR48/WAwoeIiPicN/1N+nVp47+JBaLSIrDvg5/2gX3vyZ/7wJ5z6kxGZZlnnxVpg9gOEJsELdtBy/an/Txtu0Xbxj2muqbZpH+6iIiEJPU3OckwzDMSP/14KlQ4f+41f3p6xqJFWzNYxCWf9jMJYpMhroO5r1lrsAT+mSSFDxER8bmw6m/icJiXP47tgfzvIf/kz2N7zO1SD7q12+Ih/mxolQLxKSd/djwVLFomQWR04x+Lnyh8iIiIz4VcfxPDgIIcOLIDju42Q0VV2Dj2g3n3SG1atj8tVKRAq7Ndn9vCq66VwoeIiPhcVX+TsfPXYwGXABLQ/U0qSs1wcWQHHNkJR7af3N4F5cU1v88SYQaKhHMgobP5s3XVz04Q1cx/xxAEFD5ERKRRVPU3ObPOR1Ig1PkoOw6Hv4NDW82fR3bC4e3m2gzD4f491siT4aKLa8hI6GyevYiI8u8xBDGFDxERaTRN3t/E4TAvjxzaCge3mI9DW83LJTWFjJg4SDzv5KMrtO1mbrdOVcDwEYUPERFpVBFWi39upy07Dgc3Q+5GOPCtGTIObYPy4+7HN28D7X8GbS+AtudB4smQ0bJdUNwxEswUPkREQlDI91MpLYIDmyBvI+R9YwaOI9vdn82ItJlnL9r9zAwb7bub2woZTUbhQ0QkxIRcP5XyE2bA2L/2VNg4shPc3UfTsj10uAg69DwZNHqY6zKsEX6etNRG4UNEJIQEfT8VwzAXfe77GvZ/DfuzzTMcjorqY2OTocPPIfmik4Hj52ZNDAl4Ch8iIiEiKPuplB2HnHUng8ZaM2wUH64+rkVb6NgXknufChwt2/l9uuIbCh8iIiEiKPqplBTAvmz48XP48QszeDjKXcdYo8zLJh37nnq0OlvrM0KIwoeISIgIyH4qx/Nh75dm0Pjxc3O9xpmLQlsmwdmXnAwaaWbwUFGukKbwISISIgKin0r5CTNsfJ8Fuz8xb3k9U6tOkHo5dLrMfLTurLMaYUbhQ0QkRDRJPxVHpRkwdn9iBo69X0FlqeuYxG4ng8YvoFM/s2GahDWFDxGREOG3fiqFB2Hnctj1EexZXb0lfGwydOkP5/SHzldCbPuG/XkSchQ+RERCSKP0U3E4IG8D7FgOOz4wa22cLiYOUq+Ac9LN0NHmXF1GkVopfIiIhBif9FMpKYDdK80zHDs/guJDrq8n94auV8O5A83tCP06Ec/pX4uISJCrqZS617fTFh2C75bBtqXm5ZTTC3tFx5pnNc67Bs4drEsp0iAKHyIiQazBpdR/2gvbTgaOvV/islKkzbnQ9RozcJzdDyKjfX8AEpYUPkREglS9S6kf3Q1b/wNbl1Rfv5HcG7pfB+cPg8RzG2nmEu4UPkREgpDXpdQL8mDLu7BpMeSuP220xbwN9oLr4PxroVWKfw5AwprCh4hIEPKklHqx/Qg/LJ9Dl4Pvw55PcV5SsUTAOVdB9+HQbah6pIjfKXyIiAShmkqkR1POIOs6ro/4nHTrRqK/qjz1YsolcOFI6H49tGzrn4mKuKHwISIShM4skd7d8gMjI1ZxfcTntLYUOfcXt+pGiz6/hh43QetO/p6miFsKHyIiQSitcwLd4srpV7ySkRFZ/Mz6o/O1PCOBdysv54tm/Xntd3dCQyuaiviYwoeISDBxOOD7lUSsf533K97DGmW2oy81IvnIcTGLKq/iM8eFOLAyZ3jvhpdSF2kECh8iIsHgeD5smA9rX4JjewCwAvZW3fln0WW8XtSXn4gFvKzzIdIEFD5ERAKVYUDOOvj6n7D53VPdYmPi4ee/hl6jiO/Qk/EOg8saUkpdxM8UPkREAk35Cdi0yAwded+c2p/UE/reDReOgOgWzt31KqUu0oQUPkREfKSmHiseKzpkBo6v/wnHj5r7ImKgx41m6Dirj7rFSkhQ+BAR8YEG9Vg5uBW+mg3fvg2VZea++LMh7W7odRs0T2jEmYv4n8KHiEgD1avHimHA7hXw5WyzdX2Vsy6Gy+43e6uoTb2EKP3LFhFpAK97rDgqYdsS+PTPcGCTOdBihQuGQb/7ISXNj7MXaRoKHyIiDeBJj5U8ewlf7z7IpcUr4dPn4OhO88WoFtDnDrjkXmid6pf5igQChQ8RkQaoqcdKlRjKGBmxip7vPgwncs2dtni4ZIz50HoOCUMKHyIiDXBmj5UqMZRxa8QKxkQupZ3lJzgBtGhrXlq5+Ldgi/PrPEUCicKHiEgDpHVOoEO8jQP2Egwgigp+FZHF/ZH/oYMlH4ADJNIu4xGsfW6HqGZNOl+RQGBt6gmIiASzCKuFScO6Y8XBTdbVrIj+A09HvUQHSz45Rhsyy+/mmxtXYr30XgUPkZN05kNEpCEcDjL4km/aTqFl4fcAHDbi+XvF9XzSYigTR/6ca9RjRcSFwoeISH19nwXLn4AD39ISMJq1Zu8F97Ap+VdktG7Nn9RjRcQthQ8REW8d+g4+egJ2LgfAiI5l//m/5duUW0lISGSIQodIrRQ+REQ8VXgQsqbB+lfBcIA1kh/P+TX37h3Ed9nRkL0L2KWW9iJ10IJTEZG6lJ+AVTPh//WGdS+bweP8X7J68FLSNw/lu4Jol+FVZdU/2JzXRBMWCWwKHyIiNTEM2LYMZqfBJ1OhrAiSe8Od71P5q/k8mnWixrLqYJZVr3S4GyES3nTZRUTEnSM74f1HzeZvAHFnweCn4Gc3gtVK9u6jHpVVz96TT78ubfwzZ5EgofAhInK60kJYPRO+fAEc5RARDZc9AFf8AaJbOIfVVVbd23Ei4UThQ0QEzEssW/4NH06EwpNrNbpeAxnToE2XasNrKqte33Ei4UThQ0Tkp73w3sOw80PzeevOkDEdumXU+JYzy6qfyQIkxdtI66zGcSJn0oJTEQlflRXw5WyYfYkZPCKi4aoJMO6rWoMHnCqrDmbQOF3V80nDuqveh4gbCh8iEp5yN8I/B5iXWcqPw9mXwZjPoH8mRHl2qSSjRwfmjOpNUrzr+KR4G3NG9VadD5Ea6LKLiISX8hPwydPmGQ/DAbZ4GDwFet0GVu//eyyjRwcGd08ie08+hwpLaBdrXmrRGQ+RmnkdPoqLi3nkkUdo3bo1RUVFzJgxg5iYmGrjnnnmGRwOB1arlePHjzNlyhQsFv2fUUSa0P618J+xcGSH+bzHTXDNNIht36CPjbBadDutiBe8jvljx45l0KBBTJ06ld69e5OZmVltzLJly9izZw+PP/44EydOJDc3l7ffftsnExYR8VpFKXz8JPxrsBk8WibBLW/BiJcgtj2VDoMvdx/lvxtz+HL3URUGE2lkXp35yM3NZdGiRfzjH/8AYMiQIYwZM4bJkycTGxvrHLdt2zYKCwudz5s1a4bdbvfRlEVEvJC7Af49Fg5vM5/3vNm8k6W5eRfKB5vzmLx0q0vBMPVmEWlcXp35yMrKIjExEZvNXFzVtm1boqOjyc7Odhl33XXXsWTJEhYvXsyxY8c4cuQIt912m9vPLC0tpaCgwOUhItJgleWw8mmYN9AMHi3aws3z4cZ/uASPsfPXV6tUqt4sIo3Lq/CRk5NDQoLrPeuxsbHk5ua67OvWrRtvvvkmt99+O7/5zW94+eWXadasmdvPnDZtGvHx8c5HSkqKl4cgInKGo7vhpWtg9bNgVMLPboBxa+CCYc4hlQ6DyUu3qjeLSBPwKnxYLBbnWY8qZWVlREVFVRtbUlLCkiVL2L17N7fccgsVFRVuPzMzMxO73e587Nu3z5spiYicYhiw8U148UrIWWfeyTLiJRj5CrRwXRCavSff494sIuJbXoWP5OTkams3ioqKSE5OdtmXlZXF2rVrGTRoEFlZWWzYsIHnnnvO7WfGxMQQFxfn8hAR8dqJn+Cdu8y7WcqKoNMvYMzn5h0tbqg3i0jT8Sp8pKens3//fsrKygCcl1vS0tJcxi1evJiuXbsCZmCZOnUqq1ev9sV8RUSq+/FLmHsFbH4HLBEw4HG4Yym0qvkyrnqziDQdr898ZGRksGrVKgCWL1/OuHHjsNlszJo1i+3btwNw0UUXsWHDBuf7LBZLtYAiItJgDgesmgmvDAX7XmidCr/9EK78I1gjan1rVW+WmqoPWTDvelFvFhHf87rI2Ny5c5kwYQJr1qwhPz+f6dOnA7BgwQJSU1Pp1q0bv/3tb9m7dy/Tp08nMTGRffv2MXHiRJ9PXkTCWPFR+Pc9sOtj83nPm2HoLLB5dum2qjfL2PnrsYDLwlP1ZhFpXBbDMAJqKXdBQQHx8fHY7Xat/xAR9/Z9DYv+Dwr2Q6QNrv0z9BpVr49SnQ8R3/Dm97d6u4hI8DAMWPMiLH8cHOWQ0AV+9Rok9aj3R6o3i4j/KXyISHAoKYAlD8DW/5jPuw+H6/7u8WWW2qg3i4h/KXyISOA7shMW3AJHd4I1Eq5+Gi65F85oVlnpMHQGQyQIKHyISGDb+REsvgtK7RB3llkwLKX63XNauyESPLzuaisi4heGAZ/9Bd4YaQaPlEvhnqwag4d6tIgED4UPEQk85Sfg3dHw8STAgN53mEXDWrarNlQ9WkSCj8KHiAQW+36zKdymReb6jmv/DMP+CpHRboerR4tI8NGaDxEJHPvXwoJfQ/FhaN7GvI029fJa36IeLSLBR+FDRALD1iXmpZaKEmh/IdzyJrQ6u863qUeLSPDRZRcRaVqGAV/8P3j7djN4dL0GfvuBR8ED1KNFJBgpfIhI06msgPf+YFYsxYC+o+HXb0JMS48/oqpHC1AtgKhHi0hgUvgQkaZRWgQLb4G1/wIscM0zMHQmRHh/NTijRwfmjOpNUrzrpZWkeBtzRvVWnQ+RAKM1HyLifwV58Oav4MC3ENkMbpoHFwxr0EeqR4tI8FD4EBH/OrobXr8eftoLLdrCLW9Bxz71+ih35dTVo0Uk8Cl8iIj/5H0D828yb6VNOAdGvQsJnev1USqnLhK8tOZDRPxjz6fw8rVm8EjqCb9d3qDgoXLqIsFL4UNEGt+2ZeYZj7JCSL0C/u89aNm2Xh+lcuoiwU/hQ0Qa1/rX4e3boLIUzv8l3LoYbHH1/jiVUxcJfgofItJ4Pv8bLLkfDAf0ug1GvgpRDas0qnLqIsFPC05FpHGsehY+edrc/sV4GPQkWBp+26vKqYsEP535EBHfMgxYOfVU8Bj4Jxg82SfBA1ROXSQUKHyIiO8YBnz0J1g903x+9dNwxR98+keonLpI8FP4EBHfMAz4IBO++Jv5fMhMuOz+RvmjVE5dJLhpzYeINJzDAf97+GSfFuCXf4GL72zUP1Ll1EWCl8KHiDSMwwFLfwcbXgcsMPzv0GuUX/7oCKtF5dRFgpDCh4h4zdlTpeA4l26ZSvtdC8Fihevnws9vburpiUiAU/gQEa+c6qlygsmRr9A+8iMqsbIpbSYXKXiIiAe04FREPHaqp8oJHot8gzsiP8JhWHik7B5uWNVBPVVExCMKHyLikVM9VQz+GPkWoyP/B0Bmxd2847gSUE8VEfGMwoeIeKSqp8rvIv7NfZFLAHii/P94q7I/oJ4qIuI5hQ8R8cihwhLuiVjKQ1GLAZhSPorXK692O05EpDZacCoiHulx4L8Mj1oAwLPlN/OvyqFux6mniojUReFDROq2dQnnrHkMgDkVw3ihcni1IRbMCqPqqSIiddFlFxGp3fdZ8M5dWAwH+zqP5NmKX6uniog0iMKHiNQsZx0svBUqy+CC60i57UXmjOqjnioi0iC67CIi7h3eDvNHQFkRnJMON/0TrBHqqSIiDabwIRIGnOXQPQ0L9hx4/QY4kQ9n9YGb34DIGOfL6qkiIg2h8CES4k6VQz91C2yHeBuThnV3f5mkxA5vjISCHEg8D25dDDEt/ThjEQl1WvMhEsJOlUN3rb1xwF7C2Pnrq5dDryiDt0bBoS3Qsj2Megea6+4VEfEthQ+REHWqHHp1VftcyqEbBvz3PtizGqJbwq2LoNXZ/pquiIQRhQ+REFVVDr0m1cqhr3gKNr0N1kj41WvQ4ef+maiIhB2FD5EQ5WmZ80OFJfD1P+Gz58wdw/4G5w5sxJmJSLjTglOREOVpmfPz7F9A1h/NJ+kTodetjTgrERGd+RAJWWmdE+gQb6tWjbSKBbg89gDnfz4eDAf0ug2uesSPMxSRcKXwIRKiIqwWJg3rDuC2HHoiduZFzcJSVgSdr4RfPg8WFQoTkcan8CESwjJ6dGDOqN7VyqGfHWdlefKLNDueCwldYOSrEBHVRLMUkXCjNR8iIa5aOfSWMVzyTSbWTRvBFg+/eUu1PETErxQ+RMKASzn01bPMW2otEeYZj8SuTTs5EQk7Ch8iQcLr/izubP0vrJxibg+dCV36+36iIiJ1UPgQCQJe92dxJ3cjvHuvuX3JGOh7l+8nKiLiAS04FQlwXvdncafoMCz8DVScgC4D4eqnG2m2IiJ1U/gQCWBe92dx+yHlsOj/zC61bbrCyJchQic9RaTpKHyIBDCv+7O489Gf4MfPzGZxv37DvMNFRKQJKXyIBDCv+rO4881b8NUL5vYNc6FtNx/NTESk/hQ+RAKYp/1Z3I7L+waW/s7cvuJhuGCYD2cmIlJ/Ch8iAcyT/iwd4s3bbl0UH4WFo6CiBM4dDP0nNvZURUQ85vWqs+LiYh555BFat25NUVERM2bMICYmxu3Y48eP8+KLL9KmTRvOO+88Lr300gZPWCScVPVnGTt/PRZwWXhaFUgmDevuWu+jsgLe+S3Y90LrznDTPLBG+HHWIiK18/rMx9ixYxk0aBBTp06ld+/eZGZmuh2Xn5/PyJEjuf7667n99tsVPETqqab+LEnxNuaM6l29zsfKp+D7LIhqbi4wbdbaf5MVEfGAxTCMWu7Rc5Wbm0uXLl04duwYNpuNw4cP06lTJw4ePEhsbKzL2KuvvpqJEyeSnp7u1YQKCgqIj4/HbrcTFxfn1XtFQplHFU63LYW3RpnbI16GHjf6f6IiEpa8+f3t1ZmPrKwsEhMTsdnM/wJr27Yt0dHRZGdnu4xbtmwZO3fuJDs7m6FDh5KZmUl5ebnbzywtLaWgoMDlISLVVfVnGX7RWfTr0qZ68MjfA/+5z9zud7+Ch4gELK/WfOTk5JCQ4LqwLTY2ltzcXJd9b7zxBpdccgnjx4/nrrvuolevXhiGwfTp06t95rRp05g8eXI9pi4SHjw641FRahYSK7VDxzQY9GRTTFVExCNenfmwWCzOsx5VysrKiIqKctm3ZcsWLr/8cqKjo2nTpg2jR4/mtddec/uZmZmZ2O1252Pfvn1eHoJI6Ppgcx6Xz1jJLfO+4vcLN3LLvK+4fMbK6iXVP3wM8jZCs4STFUyj3H6eiEgg8Cp8JCcnY7fbXfYVFRWRnJzssq+iooLKykrn8549e5Kf774CY0xMDHFxcS4PEfGip8vmd+Dreeb2jf+A+I5+nqmIiHe8Ch/p6ens37+fsrIyAOfllrS0NJdxPXv2ZOfOnc7nkZGRnH/++Q2dq0jY8Liny+GdsOT35o7LH4Kug/01RRGRevP6zEdGRgarVq0CYPny5YwbNw6bzcasWbPYvn07AL///e957733KC0tBeCLL77ggQce8PHURUKXJz1d8u0FlLw5CsoKodMvoP9j/pugiEgDeF1kbO7cuUyYMIE1a9aQn5/vXES6YMECUlNT6datG/369WPKlCk88MADdOvWDZvNxl133eXzyYuEKk96ukyKfJUWx76DFm3hpn+pU62IBA2v6nz4g+p8iMCXu49yy7yvanz9Ouvn/C16NgYWLLf/B85J99vcRETcabQ6HyLiH7X1dEmxHOTpqJcAMK58RMFDRIKOwodIAKrq6QK4BJBIKvhb1GxiLSfIb9MH61WPNM0ERUQaQOFDJEC56+ny+8h36WXdRXlULAm3vap1HiISlPTNJRLAMnp0YHD3JLL35FO551N+8dl/AYga/jdoldLEsxMRqR+FD5EAdXpZ9eToEi79diIWDOg1Sn1bRCSoKXyIBKAPNucxeenWk7U+DOZE/QVLRA7FLVNpkTGjqacnItIgWvMhEmDOLKv+64hPGBLxNWVGBL8+ejcf7Cxs4hmKiDSMwodIADmzrHoXSw6TIs2mjDMrbmazcY5ZVt0RUOV5RES8ovAhEkBOL6seRQV/jZpNM0sZqysv5J+VQzGAPHsJ2XvcN2oUEQkGCh8iAeT0suq/i3yXHtYfyDda8ofyMRin/d/Vk/LrIiKBSuFDJIC0izVrevSy7GRchHlb7WPld3GY1m7HiYgEI4UPkQCS1jmBznHwXNQcIiwG71ZezvuOS5yvW4AO8TbSOic03SRFRBpI4UMkgERYLbySsozO1gPkGQk8WX6H87WqMuuThnUnwuqu64uISHBQ+BAJJLs+ptPuNwGYFv0ABbRwvpQUb2POqN5k9OjQVLMTEfEJFRkTCRQnjsF/7ze30+7l+YwHueVkhdN2sealFp3xEJFQoPAhEijeexgK86DNuTDoSSKsFvp1adPUsxIR8TmFD5FAsPkd2LwYLBFwwz+ojGxG9u6jOushIiFJ4UOkqRUehPf+YG5f+TAf/JTM5NdXOouNgXmHy6Rh3bXeQ0RCghacijQlw4D3HjLXe3T4OR+2GeXS16XKAXsJY+ev54PNeU00URER31H4EGlKm9+B75aBNYrK62bz5Hs7cde1pWqf+rqISChQ+BBpKkWH4X9/NLevfJjs48nVznicTn1dRCRUKHyINJX/PQwn8qF9D7j8IY/7taivi4gEO4UPkaaw9b+w9T/m3S3DZ0NktMf9WtTXRUSCncKHiL8VHz11d8vlD0LyRYDZ16VDvI2abqhVXxcRCRUKHyL+9sGjUHwY2l4AVz3i3B1htTBpWHeAagFEfV1EJJQofIj403f/g02LwGI9ebklxuXljB4dmDOqN0nxrpdW1NdFREKJioyJ+MuJY7DsQXO73/3QsY/bYRk9OjC4exLZ6usiIiFKZz5E/GX5E1B0wOzd0n9ijcMqHYaCh4iENJ35EPGHPZ/ChtfN7ev+DlHN3A77YHMek5duVWl1EQlpOvMh0tjKS2DZeHO7z53QqZ/bYR9szlNpdREJCwofIo3t0z/D0V3Qsj0MetLtkEqHweSlW1VaXUTCgsKHSGM6tA0+e97cHvIsNGvldlj2nnyVVheRsKHwIdJYHA5Y+ntwlEO3odB9eI1DVVpdRMKJwodIY1n3MuxbA9EtYehMsNR8x4pKq4tIOFH4EGkMBXnw8ZPm9sA/QXzHWoertLqIhBOFD5HG8P4fobQAzuoDfe+uc7hKq4tIOFH4EPG1bctg21KwRsKwv4E1wqO3qbS6iIQLFRkT8aWSAvjfH83ty34HST28ertKq4tIOFD4EPGllVOhMBdad3bpWOuNCKuFfl3a+HhiIiKBQ+FDxFdyN8LX88ztXz5fYwn106mPi4iEI4UPEV9wOOC9h8BwQI8R0KV/nW9RHxcRCVdacCriC+tfhZx1EB0L1zxd53D1cRGRcKbwIdJQxUdO1fQY8DjEJtU6XH1cRCTcKXyINNRHk6DkJ0i60KOaHurjIiLhTuFDpCF+/BI2zje3r30eIupeRqU+LiIS7hQ+ROqrstxcZArQ+w5I6evR29THRUTCncKHSH2teREObYVmCTDoSY/fpj4uIhLuFD5E6sOeA1nTzO3BT0Fzz4OC+riISLhT+BCpjw8nQlkRpFwCF93q9dvVx0VEwpmKjIl4a9fHsPU/YImAa58Da/0yvPq4iEi4UvgQ8UZFKfzvZM+WS8bU2jjOk9Lp6uMiIuFI4UPEG1/Ohvzd0DIJ0ifUOEyl00VEaqY1HyKeKsiF1bPM7cFPgS3O7TCVThcRqZ3Ch4inPvoTlBebi0x7/srtEJVOFxGpm8KHiCd+/AI2LQIsMORZsLhfFKrS6SIidVP4EKmLo/LUItM+/wfJF9U4VKXTRUTq5nX4KC4u5r777uPxxx9n/PjxlJaW1jp+4cKFpKen13d+Ik1v3ctwcBPYWsGAJ2odqtLpIiJ18zp8jB07lkGDBjF16lR69+5NZmZmjWNzc3N5+umnGzRBkSZ1PB9WTjW3BzwOLWq/LVal00VE6uZV+MjNzWXRokUMGTIEgCFDhjB37lwKCwvdjp8+fTpjx45t+CxFmsrKKXDiGLTvAX3urHO4SqeLiNTNq/CRlZVFYmIiNpt5yrht27ZER0eTnZ1dbey8efMYNWoUzZs3r/UzS0tLKSgocHmIBIS8b2Dty+b2kGchwrOyOCqdLiJSO6+KjOXk5JCQ4Hq6ODY2ltzcXJd9u3btoqCggLS0NLZu3VrrZ06bNo3Jkyd7Mw2RxmcYJxeZGtDjJkj9hVdvV+l0EZGaeRU+LBaL86xHlbKyMqKiopzPKysreeGFF5g5c6ZHn5mZmclDDz3kfF5QUEBKSoo30xLxvU2LYN9XENUcBk+p10eodLqIiHtehY/k5GTsdrvLvqKiIpKTk53Pv/jiC+bOnctLL70EmOGkrKyMVq1a8dNPP1X7zJiYGGJiYuoxdZFGUloIy0/e1XLlwxB/lsvLnvRsERGRmnkVPtLT07nnnnsoKysjOjraebklLS3NOaZv374ul1oWL17M4sWLWbhwoY+mLNLIVs+CogPQujP0u9/lJfVsERFpOK8WnCYnJ5ORkcGqVasAWL58OePGjcNmszFr1iy2b9+OzWYjNTXV+ahaoJqamtoY8xfxraO7zeZxABnTIfLUWTn1bBER8Q2v63zMnTuXt956i6lTp/Ltt98663gsWLCATZs2+XyCIn710Z/AUQ7nDoJuGc7d6tkiIuI7FsMwAurbsqCggPj4eOx2O3Fx7ruGijSK71fBa9eBJQLGfQltuzlf+nL3UW6Z91WdH7Fg9KVaZCoiYcmb39/q7SICZv+WDyea233vdgkeoJ4tIiK+pPAhArD+NTi42ezfkj6h2svq2SIi4jsKHyIl9lP9W9IzoXn1vivq2SIi4jsKHyKrZ8HxI5B4HvS9y+0Q9WwREfEdhQ8Jb0d3w1dzzO1rnoGIqBqHqmeLiIhveFVkTCTknH5rbdfBdQ5XzxYRkYZT+JDw9f0q+G6ZeWvt1U97/Db1bBERaRiFDwlPLrfW3gXtzq/zLerpIiLiGwofEp42vH7y1tp48w6XOqini4iI72jBqYSfEjusmGJu13Br7enU00VExLcUPiT8VN1a26arWc20FurpIiLiewofEl68uLUWIHtPfrUzHqczgDx7Cdl78n04SRGR0KbwIeGl6tbaLgM9urVWPV1ERHxP4UPCx57Vp26tveZpsNR9p4p6uoiI+J7Ch4QHRyV8cPLW2ot/C+0u8Oht6ukiIuJ7Ch8SHja8Dgc3mbfW9p/o8dvU00VExPcUPiT0lRTU2bW2NurpIiLiWyoyJqHvs+eg+DAkdIGL3XetrYt6uoiI+I7Ch4S2Yz/Cly+Y21dPhcjoOt9SUxl19XQREfENhQ8JbSsmQ2UpdL4Sug2pc7jKqIuIND6t+ZDQtS8bNr8DWMyutXXcWqsy6iIi/qHwIaHJMOCDkw3jet0KHXrWOlxl1EVE/EfhQ0LT5ncgZy1EtYABT9Q5XGXURUT8R+FDQk/5Cfj4SXP78gchNqnOt6iMuoiI/yh8SOj56gWw74O4s6DffR69RWXURUT8R+FDQkvRIfj0OXN74CSIbu7R21RGXUTEfxQ+JLR88jSUFUFyb7hwpMdvUxl1ERH/UfiQ0HFwC6x/zdy+5hmwevfPW2XURUT8Q0XGJDQYBnw4EQwHdB8OnfrV62NURl1EpPEpfEho2PkRfJ8FEdEwaHKDPkpl1EVEGpfChwS/ynJY/pi5fckYSOhc89Aa+raIiIj/KHxI8Fv3ChzZAc3bwBV/qHGY+raIiAQGLTiV4HbiJ/jkGXM7PROatXI7TH1bREQCh8KHBLdPZ8GJfEjsBn3udDtEfVtERAKLwocEr/zv4au55vY1T0OE+6uI6tsiIhJYFD4keH00CRzl0GUAnDuoxmHq2yIiElgUPiQ4/fgFbFsCFitcPRUsNd+xor4tIiKBReFDgo/DYRYUA+h9O7T/Wa3D1bdFRCSwKHxI8Nn0NuRugOhY6P9YncPVt0VEJLAofEhwKTsOH5+sYHrFQ9CynUdvU98WEZHAoSJjEly+/DsU5kL82XDpOK/eqr4tIiKBQeFDgkdBHnz2vLk9+EmIqr5AtK7y6erbIiLS9BQ+JHisnArlx6FjX/jZjdVeVvl0EZHgoDUfEhzyvoGNb5jb10yrdmutyqeLiAQPhQ8JfIYBHz4GGNDjJkjp6/KyyqeLiAQXhQ8JfN8tgx8+hYgYGPRktZdVPl1EJLgofEhgKy+B5Y+b25c9AK3OrjZE5dNFRIKLwocEtq9egGM/QMskuPxBt0NUPl1EJLgofEjgKjwAn/7Z3B48GWJauh2m8ukiIsFF4UMC14qnoKwIzuoDF/6qxmEqny4iElwUPiQw5aw7dWttxgyw1v5PVeXTRUSCh4qMSeAxDHh/grnd8+Zqt9bWROXTRUSCg8KHBJ5Ni2F/NkQ1d3trbW1UPl1EJPApfEhgKSuGjyeZ21c8BHHJQN09W0REJHgofEhg+fyvUJBjdq3tdz+gni0iIqFGC04lcPy0zwwfAFdPgahm6tkiIhKCvA4fxcXF3HfffTz++OOMHz+e0tLSamPsdjsjR44kLi6OXr168dVXX/lkshLiPvoTVJRAp19A9+Hq2SIiEqK8Dh9jx45l0KBBTJ06ld69e5OZmVltzIwZM7jxxhv55JNPSElJYfjw4RQXF/tkwhKifvwCtrwLWCBjOlgs6tkiIhKivAofubm5LFq0iCFDhgAwZMgQ5s6dS2Fhocu4QYMGccstt9CnTx/mz5/PsWPH2Lp1q+9mLaHFUQnvP2pu97kDOvQE1LNFRCRUeRU+srKySExMxGYzCzm1bduW6OhosrOzXcYNGDDAuR0XF0dcXBwdO3Z0+5mlpaUUFBS4PCTMrH8NDnwLMXHQ/3HnbvVsEREJTV6Fj5ycHBISXPtjxMbGkpubW+N7duzYQXp6Oh06uL8rYdq0acTHxzsfKSkp3kxJgt3xfFgx2dxOz4SWbZ0vqWeLiEho8ip8WCwW51mPKmVlZURFRdX4ntmzZzNz5swaX8/MzMRutzsf+/bt82ZKEuxWPAUnjkG77pB2j8tL6tkiIhKavAofycnJ2O12l31FRUUkJye7Hf/hhx8yYMAAOnfuXONnxsTEOC/NVD0kTOSsh3WvmNtDZ0FE9bIz6tkiIhJ6vCoylp6ezj333ENZWRnR0dHOyy1paWnVxm7ZsoW9e/cyevRo38xUQovDAf/7I2DAhSMh9Rc1DlXPFhGR0OL1mY+MjAxWrVoFwPLlyxk3bhw2m41Zs2axfft2AHbv3s28efMYPHgwP/zwA9988w1z5szx/ewleG2cDzlrIToWBk+pc3hVz5bhF51Fvy5tFDxERIKY1+XV586dy4QJE1izZg35+flMnz4dgAULFpCamkpcXBzp6ens37+fv/71r873LVy40HezloDidd+V4/nw8ZPmdvoEiNOlExGRcGIxDCOgykMWFBQQHx+P3W7X+o8gUK++K8segrX/grYXwJhPIaLmBcsiIhIcvPn9rd4uUm/16ruSuxHWvmRuD52p4CEiEoYUPqRe6tV3xeGA/z1sjugxAjpf4YeZiohIoFH4kHqpV9+Vb96E/V9DdEu4emrjT1JERAKSwofUi9d9V4qPwvInzO2rHtUiUxGRMKbwIfXidd+V5Y/BiXxo3wMuHduIMxMRkUCn8CH14lXfle9XwTcLzL3D/qpFpiIiYU7hQ+rF474rlaWw7EFzR9+7oePFfpujiIgEJoUPqTeP+q58+mfI3w0tk2DgE000UxERCSReVzgVOV2tfVcOb4fPnjcHDn0WbPFNO1kREQkIOvMhjcPhgKXjwVEO52XABdc19YxERCRA6MyHNEhN5dX/0WMrF+79AqKam5VMLWoEJyIiJoUPqbeq8upnVjmtsB8gZe00c+Vp/8eg1dlNMT0REQlQuuwi9VJzeXWDqVEv0cpSzHeWc6hMu7cJZiciIoFM4UPqpaby6sOsX3JNxFrKjQjGl9xD9o8FTTA7EREJZAofUi/uyqsnUMCTUa8CMLtyON8ZZ3tchl1ERMKHwofUi7vy6k9FvUIbSyHbHGczu+L6GseJiEh4U/iQejmzvPo11mx+GfEVFYaVh8vvpYLIU+XVRURETqPwIfVyenn11hQyNeolAOZWDmOr0Rk4WV7dqltsRUTElW61lXrL6NGBObf2IuLfd9PWUcAOx1n8reJGkuJtTBrW3SyvLiIicgaFD2mQDOMzcHyOYYkgr/9zvHp231Pl1UVERNxQ+JD6+2kvvPcHACzpE7jqqoymnY+IiAQFhQ9xqnQY7hvEueOohH+PhdIC6NgXLn/Iv5MVEZGgpfAhQM09Wmpcu/Hl3+HHzyCqBdz4D4jQPyUREfGM7nYRZ4+WMyuWHrCXMHb+ej7YnOf6hgObYMUUc3vIdEg4x08zFRGRUKDwEeZq7tGCc9/kpVupdJx8VnYc3hkNjnI4/5fQ6zZ/TVVEREKEwkeYq6lHSxUDyLOXkL0n39zxwaNweBu0aAfD/goW3dUiIiLeUfgIc572XjlUWALfLoL1rwEWuGketEhs3MmJiEhIUvgIc572XjnbyIVl480nVz0C56Q32pxERCS0KXyEuTN7tJzJAnSKs3LRV+OhrAhSr4CrHvXjDEVEJNQofIS503u0nBlAqp7PP+sdLAc3Q/NEuHEeWCP8OkcREQktCh9i9mgZ1ZukeNdLMEnxNpZetpOUPYsAi1nPI079WkREpGFUGUoAM4AM7p7kWuE0cicRr041Bwz8E5w7sGknKSIiIUHhI4y5K6fer0sb88WCPPjH7WY9j+7D4fIHm3ayIiISMhQ+wlSt5dTPT4C3b4eig9CuOwx/QfU8RETEZxQ+wlBVOfUzq5qa5dTX8fl5C0nenw22eLh5PsS0bJJ5iohIaNKC0zBTVzn1ByMXk7x3KYY1Eka+Am26+HmGIiIS6hQ+wkxt5dRHRKzid5H/BuD7tCnQZYA/pyYiImFC4SPM1FROvZ91C9Mi/wnA3yuGszlpuD+nJSIiYUThI8y4K6fe07Kbf0Q9R5SlkiWV/fhzxUiPy66LiIh4S+EjzJxZTv08yz5ejZ5BrOUEX1Z254/l95IU35y0zglNOk8REQldCh9h5vRy6p0tebwR/QytLUVscJzL6PI/UEY0k4Z1J8KqW2tFRKRx6FbbMJTRowMv39Cebu//jrbY2eY4mzvKHiE2vjWzhnUno4dKqIuISONR+AhHR3eT/vkdwBFOxHXmx8te5cW2HUnrnKAzHiIi0ugUPsLNwS3w2vVQfAjadKXZ7f8lI/6spp6ViIiEEYWPEOGuT0u1sxj718H8G6HkJ2h/Idz2b2jZtknmKyIi4UvhIwTU2qelav3GjuWw+E4oK4KOfeHWRdCsdRPNWEREwpnudglyVX1azqxaavZpWc8Hm/Mgex4suNkMHp2vNM94KHiIiEgT0ZmPIFZXn5ZoKih690FwvG/uvGgU/PJ5iIz25zRFRERcKHwEsdr6tCRxlBei/0pvxy5zx8BJcPmDYNHdLCIi0rQUPoJYTX1a+ls3MCtqLm0shdiN5my7dCaXXjHKz7MTERFxT2s+gtiZ/VdacIJnIufxcvRM2lgK2eRI5ZdlT2OcN6SJZigiIlKdznwEsao+LQfsJVxl3cCUyFdIsR7GYVj4V+UQZlX8ioT4OPVpERGRgKLwEcQirBamD4ijdNnTXB2xDoB9jrY8XD6GbOMCAPVpERGRgKPwEayKDsOnf+aqtf+CiDIqiOBfFRn8reJGimlWvc6HiIhIgFD4CDY/7YM1c2Hty1BebO7rfCWWa2bQs7g9z9RW4VRERCQAKHwEA0clld+v5tin/yRh7/tYjUpzf3JvGPgnOCedCIuFfk07SxEREY94HT6Ki4t55JFHaN26NUVFRcyYMYOYmJhq45YsWcKKFSsoLS1lxIgRDBo0yCcTDhuOSshZB9+9x4n1b9HsRB6JJ1/6rPJnvBtzPVf3G0VGl+QmnaaIiIi3LIZhuCuQWaPbb7+dG264gRtuuIHXXnuNjRs38txzz7mM+e6777j99ttZs2YNhmFw8cUXs3TpUs46q+7uqQUFBcTHx2O324mLi/PuaIKZoxIObYW9X5mPPaug+LDzZbvRnKWV/XizciBbjVSqLqjMGdVb6zpERKTJefP726vwkZubS5cuXTh27Bg2m43Dhw/TqVMnDh48SGxsrHPcmDFjaNeuHU899RQA9913H61bt2bq1Kk+nXzQqSiFooNQeMB85H8PR3bA4e/g8Haz98ppjJg4PirvyX9KerPC0ZtSXMuiW4CkeBufPTpA6ztERKRJefP726vLLllZWSQmJmKzmcWt2rZtS3R0NNnZ2QwcONA5buXKlTz66KPO5127duWdd95x+5mlpaWUlpa6TL5RHN4Ba/8FhgEYXvzEy/EGVJabQaL8OJQVm9ulhXDiWO1zjG5pdpw9+1I4ux9rKs7jnpfW1zjcAPLsJWTvyadflzYN/p9IRETEH7wKHzk5OSQkuBasio2NJTc3t9Zx7sZUmTZtGpMnT/ZmGvVj32veJdLUIqKhZRLEJkGrFEjsBm3PM38mngcRp/5KDm7M8egjayqzLiIiEoi8Ch8Wi8V51qNKWVkZUVFRtY5zN6ZKZmYmDz30kPN5QUEBKSkp3kzLM607wxV/ACwnm6vV9hPPxpkHW/01ayTExEJ0i5OPlubPlu3NVvYeNnc7s3x6Q8eJiIgEAq/CR3JyMna73WVfUVERycnJtY4rLCysNqZKTEyM27tlfK5NF/O21CByevl0dwtzqtZ8qHy6iIgEE68ay6Wnp7N//37KysoAnJdS0tLSXMYNHDiQHTt2OJ/v2rWL/v37N3SuYSfCamHSsO6A8zyLU9VzlU8XEZFg41X4SE5OJiMjg1WrVgGwfPlyxo0bh81mY9asWWzfvh0w73b5+OOPAaioqCA7O5vRo0f7eOrhIaNHB+aM6k1SvOullaR4m26zFRGRoOR1nY8jR44wYcIEUlNTyc/PZ/r06URHR9OnTx8yMzMZMWIEAC+//DJbtmyhrKyMG264weMzHyF9q20DVDoMsvfkc0jl00VEJAA1Wp0Pf1D4EBERCT7e/P726rKLiIiISEMpfIiIiIhfKXyIiIiIXyl8iIiIiF8pfIiIiIhfKXyIiIiIXyl8iIiIiF8pfIiIiIhfKXyIiIiIX3nV1dYfqgquFhQUNPFMRERExFNVv7c9KZwecOGjsLAQgJSUlCaeiYiIiHirsLCQ+Pj4WscEXG8Xh8NBbm4usbGxWCy+bZxWUFBASkoK+/btC4u+MTre0KbjDW063tAXasdsGAaFhYUkJydjtda+qiPgznxYrVY6duzYqH9GXFxcSPxFe0rHG9p0vKFNxxv6QumY6zrjUUULTkVERMSvFD5ERETEr8IqfMTExDBp0iRiYmKaeip+oeMNbTre0KbjDX3heMxVAm7BqYiIiIS2sDrzISIiIk1P4UNERET8SuFDRERE/ErhQ0RERPwq4IqM1UdxcTGPPPIIrVu3pqioiBkzZrhdPbxkyRJWrFhBaWkpI0aMYNCgQR69Fmgaerx2u527776bDz/8kC5dujBnzhwuvfRSfx+Gx3zx91tl4cKFzJ07l6ysLD/MvH58dbzHjx/nxRdfpE2bNpx33nkh/3f8zDPP4HA4sFqtHD9+nClTpvi8SrKveHq8R44c4dlnn6W4uJjZs2e7vBaK31k1HW+ofmfV9vdbJRi+s+rFCAG33Xab8e677xqGYRivvvqq8eCDD1Ybs23bNqNv376Gw+EwKisrjV69ehn79++v87VA1NDjzczMNN58801j7dq1xrBhw4x27doZRUVFfj0GbzT0eKvk5OQYPXr0MK666ip/TLvefHG8R48eNYYOHWp8//33fpt3QzT0mJcuXWrcfffdzrF33nmnsXDhQv9Mvh48OV7DMIx169YZ1113nXHHHXe47A/F7yzDqPl4Q/E7yzBqPt4qwfKdVR9BHz5ycnIMm81mnDhxwjAMwzh06JDRrFkzo6CgwGXcvffeazzxxBPO5+PGjTMee+yxOl8LNL443hUrVjj32+12IyoqysjOzvbD7L3ni+Ot8sADDxizZ88O6P8j++p4Bw8ebHzyySd+mXND+eKYn332WePmm292ee3FF1/0w+y95+nxVnniiSeq/XIKxe+sKu6ONxS/s6q4O94qwfCdVV9Bv+YjKyuLxMREbDYbAG3btiU6Oprs7GyXcStXrqRTp07O5127dmXVqlV1vhZofHG8AwYMcO6v6inQ2P106ssXxwswb948Ro0aRfPmzf0z8XryxfEuW7aMnTt3kp2dzdChQ8nMzKS8vNx/B+ElXxzzddddx5IlS1i8eDHHjh3jyJEj3Hbbbf47CC94erxV3DXoCsXvrCrujjcUv7Oq1NSALVi+s+or6Nd85OTkkJCQ4LIvNjaW3NzcWsedPqa21wKNL473dDt27CA9PZ0OHTo0zoQbyBfHu2vXLgoKCkhLS2Pr1q2NP+kG8MXxvvHGG1xyySWMHz+eu+66i169emEYBtOnT2/8A6gHXxxzt27dePPNN/nNb37DVVddxTvvvEOzZs0af/L14OnxevMZofCd5alQ+c6qTTB9Z9VX0J/5sFgszoRZpaysjKioqFrHnT6mttcCjS+O93SzZ89m5syZjTNZH2jo8VZWVvLCCy8wfvx4f0y3wXzx97tlyxYuv/xyoqOjadOmDaNHj+a1115r/MnXk6/+TZeUlLBkyRJ2797NLbfcQkVFReNOvJ48PV5vPiMUvrM8FSrfWTUJtu+s+gr68JGcnIzdbnfZV1RURHJycq3jCgsLnWNqey3Q+OJ4q3z44YcMGDCAzp07N96EG6ihx/vFF18wd+5c2rRpQ6tWrRg3bhyfffYZrVq18sf0veaLv9+KigoqKyudr/Xs2ZP8/PxGnHXD+OKYs7KyWLt2LYMGDSIrK4sNGzbw3HPPNf7k68HT4/XmM0LhO8sTofSdVZNg+86qr6APH+np6ezfv5+ysjIA56mttLQ0l3EDBw5kx44dzue7du2if//+db4WaHxxvGD+1/HevXsZPny4H2Zdfw093r59+7J161Y2btzIxo0beeqpp7j44ovZuHGj347BG774++3Zsyc7d+50vhYZGcn555/f2FOvN18c8+LFi+natStgfvlPnTqV1atX+2P6XvP0eGsTit9ZdQm176yaBNt3Vn0FffhITk4mIyPDudhq+fLljBs3DpvNxqxZs9i+fTsAY8aM4eOPPwbM/zLMzs5m9OjRdb4WaHxxvLt372bevHkMHjyYH374gW+++YY5c+Y0zQHVoaHHa7PZSE1NdT6qFoKlpqY21SHVyhd/v7///e957733KC0tBcz/knrggQea4Gg844tjvuiii9iwYYPzMy0Wi9e/3PzF0+OtYph3JbrsC8XvrCrujjcUv7OqnHm8wfadVV8h0dX2yJEjTJgwgdTUVPLz85k+fTrR0dH06dOHzMxMRowYAcDLL7/Mli1bKCsr44YbbnD5L4XaXgs0DTnevLw80tLS2L9/v8tnLly4kJtvvrkpDqdOvvj7rfLKK6/wyiuvBHTBHl8c7/z581m9ejXdunXj+PHjPPHEE011OB5p6DE7HA6efPJJmjdvTmJiIvv27WPixIkB26rc0+Ndt24dDz74IIcPH+bFF1/kyiuvdH5GKH5nuTveUP7Oqu3vt0owfGfVR0iEDxEREQkeQX/ZRURERIKLwoeIiIj4lcKHiIiI+JXCh4iIiPiVwoeIiIj4lcKHiIiI+JXCh4iIiPiVwoeIiIj4lcKHiIiI+JXCh4iIiPiVwoeIiIj4lcKHiIiI+NX/B8oZeqnAoBjtAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "x = np.linspace(0,0.15,500)\n",
    "plt.plot(w,pcw,'o')\n",
    "plt.plot(x,lognorm.cdf(x,paramsw[0],scale=paramsw[2]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "id": "62deed4d-84c0-4f27-a922-1a12ade12eac",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "paramsw = lognorm.fit(w,floc=0.0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 81,
   "id": "b2dd71ad-5a58-4d25-9ac9-f843ec2e945f",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.5392485735008187, 0.0, 0.0537324734468408)"
      ]
     },
     "execution_count": 81,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "paramsw"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "b5f5faee-fe4f-40bf-867a-ec8373f77f39",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
